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Preface 


Solitons were first discovered by a Scottish engineer, J. Scott Russell, in 1834 while 
riding his horse by a water channel when a boat suddenly stopped. A hump of water 
rolled off the prow of the boat and moved rapidly down the channel for several miles, 
preserving its shape and speed. The observation was surprising because the hump 
did not rise and fall, or spread and die out, as ordinary water waves do. 

In the 150 years or so since the discovery of Scott Russell, solitons have been 
discovered in numerous systems besides hydrodynamics. Probably the most im- 
portant application of these is in the context of optics where they can propagate 
in optical fibers without distortion: they are being studied for high-data-rate (tera- 
bits) communication. Particle physicists have realized that solitons may also exist 
in their models of fundamental particles, and cosmologists have realized that such 
humps of energy may be propagating in the far reaches of outer space. There is even 
speculation that all the fundamental particles (electrons, quarks etc.) may be viewed 
as solitons owing to their quantum properties, leading to a “dual” description of 
fundamental matter. 

In this book I describe the simplest kinds of solitons, called “kinks” in one 
spatial dimension and “domain walls” in three dimensions. These are also humps 
of energy as in Scott Russell’s solitons. However, they also have a topological 
basis that is absent in hydrodynamical solitons. This leads to several differences 
e.g. water solitons cannot stand still and have to propagate with a certain velocity, 
while domain walls can propagate with any velocity. Another important point in this 
regard is that strict solitons, such as those encountered in hydrodynamics, preserve 
their identity after scattering. The kinks and domain walls discussed in this book 
do not necessarily have this property, and can dissipate their energy on collision, 
and even annihilate altogether. 

Why focus on kinks and domain walls? Because they are known to exist in many 
laboratory systems and may exist in other exotic settings such as the early uni- 
verse. They provide a simple setting for discussing non-linear and non-perturbative 


xi 
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physics. They can give an insight into the dynamics of phase transitions. Lessons 
learned from the study of kinks and domain walls may also be applied to other 
more complicated topological defects. Domain walls are good pedagogy as one 
can introduce novel field theoretic, cosmological, and quantum issues without ex- 
traneous complexities that occur with their higher co-dimension defects (strings 
and monopoles). 

The chapters of this book can be approximately categorized under four different 
headings. The first two chapters discuss solitons as classical solutions, the next 
three describe their microscopic classical and quantum properties, followed by 
another three chapters that discuss macroscopic properties and applications. The 
very last chapter discusses two real-world systems with kinks and, very briefly, 
Scott Russell’s soliton. The book should be accessible to a theoretically inclined 
graduate student, and a large part of the book should also be accessible to an 
advanced undergraduate. At the end of every chapter, I have listed a few “open 
questions” to inspire the reader to take the subject further. Some of these questions 
are intentionally open-ended so as to promote greater exploration. Needless to say, 
there are no known answers to most of the open questions (that is why they are 
“open”) and the solutions to some will be fit to print. 

Every time I think about research in this area, I feel very fortunate for having 
unwittingly chosen it, for my journey on the “soliton train” has weaved through 
a vast landscape of physical phenomena, each with its own flavor, idiosyncrasies, 
and wonder. I hope that this book, as it starts out in classical solitons, then moves 
on to quantum effects, phase transitions, gravitation, and cosmology, and a bit of 
condensed matter physics, has captured some of that wonder for the reader. 

This is not the first book on solitons and hopefully not the last one either. In 
this book I have presented a rather personal perspective of the subject, with some 
effort to completeness but focusing on topics that have intrigued me. Throughout, I 
have included some material that is not found in the published literature. Prominent 
among these is Section 4.5, where it is shown that the leading quantum correction 
to the kink mass is negative. The discussion of Section 6.5, with its emphasis on 
a bifurcation of correlation scales, also expresses a new viewpoint. I had partic- 
ular difficulty deciding whether to include or omit discussion of domain walls in 
supersymmetric theories. On the one hand, many beautiful results can be derived 
for supersymmetric domain walls. On the other, the high degree of symmetry is 
certainly not realized (or is broken) in the real world. Also, non-supersymmetric 
domain walls are less constrained by symmetries and hence have richer possibilities. 
In the end, I decided not to include a discussion of supersymmetric walls, noting 
the excellent review by David Tong (see below). Some other must-read references 
are: 
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. Rajaraman, R. (1982). Solitons and Instantons. Amsterdam: North-Holland. 

. Rebbi, C. and Soliani, G. (1984). Solitons and Particles. Singapore: World Scientific. 

. Coleman, S. (1985). Aspects of Symmetry. Cambridge: Cambridge University Press. 

. Vilenkin, A. and Shellard, E. P. S. (1994). Cosmic Strings and Other Topological Defects. 

Cambridge: Cambridge University Press. 

5. Hindmarsh, M. B. and Kibble, T. W. B. (1995). Cosmic Strings. Rep. Prog. Phys., 58, 
477-562. 

6. Arodz, H., Dziarmaga, J. and Zurek, W. H., eds. (2003). Patterns of Symmetry Breaking. 
Dordrecht: Kluwer Academic Publishers. 

7. Volovik, G. E. (2003). The Universe in a Helium Droplet. Oxford: Oxford University 
Press. 

8. Manton, N. and Sutcliffe, P. (2004). Topological Solitons. Cambridge: Cambridge Uni- 
versity Press. 

9. Tong, D. (2005). TASI Lectures on Solitons, [hep-th/0509216]. 
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1 
Classical kinks 


Kink solutions are special cases of “non-dissipative” solutions, for which the energy 
density at a given point does not vanish with time in the long time limit. On the 
contrary, a dissipative solution is one whose energy density at any given location 
tends to zero if we wait long enough [35], 


lim; Max, {Too(t, x)} = 0, dissipative solution (1.1) 


where Too(t, x) is the time-time component of the energy-momentum tensor, or 
the energy density, and is assumed to satisfy Too > 0. Dissipationless solutions are 
special because they survive indefinitely in the system. 

In this book we are interested in solutions that do not dissipate. In fact, for 
the most part, the solutions we discuss are static, though in a few cases we also 
discuss field configurations that dissipate. However, in these cases the dissipation 
is very slow and hence it is possible to treat the dissipation as a small perturbation. 
In addition to being dissipationless, kinks are also characterized by a topological 
charge. Just like electric charge, topological charge is conserved and this leads to 
important quantum properties. 

In this chapter, we begin by studying kinks as classical solutions in certain field 
theories, and devise methods to find such solutions. The simplest field theories that 
have kink solutions are first described to gain intuition. These field theories are also 
realized in laboratory systems as we discuss in Chapter 9. The simple examples 
set the stage for the topological classification of kinks and similar objects in higher 
dimensions (Section 1.10), and are valuable signposts in our discussion of the more 
complicated systems of Chapter 2. 
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Figure 1.1 Shape of the ¢* potential. 


1.1 Z kink 


The prototypical kink is the so-called “Z3 kink.” It is based on a field theory with 
a single real scalar field, 6, in 1 + 1 dimensions. The action is 


1 
S= J d’x 5.0) = vo 
= J Px 5 oy —2@?- a (1.2) 
7 ee ae 


where u = 0, 1, and à and 7 are parameters. The Lagrangian is invariant under the 
transformation @ — —@ and hence possesses a “reflectional” Z) symmetry. The 
potential for @ (see Fig. 1.1) is 

4 


À 2 À À 
Vib) =F" -= a e+ 504+ 7 


(1.3) 


where m? = ànņ?. The potential has two minima: ¢ = +n, that are related by the 

reflectional symmetry. The “vacuum manifold,’ labeled by the classical field con- 

figurations with lowest energy, has two-fold degeneracy since V(¢) = V(—@). 
The equations of motion can be derived from the action 


db — p + AG? -np =0 (1.4) 


where 0, = 0/ot and similarly for əy. A solution is (t,x) = +n, and another 
is o(t, x) = —n. These have vanishing energy density and are called the “trivial 
vacua.” The action describing excitations (sometimes called “mesons”’) about one 
of the trivial vacua can be derived by setting, for example, 6 = n + yY, where yw is 
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the excitation field. Then 
1 m? X n 
S= | &x| -0y - — E 3_ Eyt 1.5 
J x z. uY) 5 wv zim av (1.5) 
where 


my =V2m (1.6) 


is the mass of the meson. 

Next consider the situation in which different parts of space are in different 
vacua. For example, ¢(t, —oo) = —nand g(t, +oo) = +n. In this case, the function 
o(t, x) has to go from —7n to +n as x goes from —oo to +00. By continuity of the field 
there must be at least one point in space, xo, such that f(t, x9) = 0. Since V(O) Æ 0, 
there is potential energy in the vicinity of xo, and the energy of this state is not 
zero. The solution of the classical equation of motion that interpolates between the 
different boundary conditions related by Z2 transformations is called the “Z2 kink.” 

We might wonder why the Z3 kink cannot evolve into the trivial vacuum? For 
this to happen, the boundary condition at, say, x = +00 would have to change in a 
continuous way from +7 to —n. However, a small deviation of the field at infinity 
from one of the two vacua costs an infinite amount of potential energy. This is 
because as ¢ is changed, the field in an infinite region of space lies at a non-zero 
value of the potential (see Fig. 1.1). Hence, there is an infinite energy barrier to 
changing the boundary condition. ! 

A way to characterize the Z2 kink is to notice the presence of a conserved current 


1 
jE = — e ad (1.7) 

2n 
where u, v = 0, 1 and e“” is the antisymmetric symbol in two dimensions (€°! = 
1). By the antisymmetry of €’”, it is clear that j is conserved: 0, j“ = 0. Hence 


1 
ge J drj? = zl = +00) — $lx = -20)] (1.8) 


is a conserved charge in the model. For the trivial vacua Q = 0, and for the kink con- 
figuration described above Q = 1. So the kink configuration cannot relax into the 
vacuum — it is in a sector that carries a different value of the conserved “topological 
charge.” 

To obtain the field configuration with boundary conditions (boo) = +n, we 
solve the equations of motion in Eq. (1.4). We set time derivatives to vanish since 


1 In Chapter 2 we will come across an example where the vacuum manifold is a continuum and correspondingly 
there is a continuum of boundary conditions that can be chosen as opposed to the discrete choice in the Z2 case. 
This will lead to some new considerations. 
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we are looking for static solutions. Then, the kink solution is 


d(x) =n ianh( Snr) (1.9) 


In fact, one can Lorentz boost this solution to get 


A(t, x) = ntanh( nx) (1.10) 


x — vt 
V1 —v? 
(Recall that we are working in units in which the speed of light is unity i.e. c = 1.) 


The solution in Eq. (1.10) represents a kink moving at velocity v. 
Another class of solutions is obtained by translating the solution in Eq. (1.9) 


Ok(x3a) =n tanh ne — o) (1.12) 


It is easily checked that translations do not change the energy of the kink. This is 
often stated as saying that the kink has a zero energy fluctuation mode (or simply 
a “zero mode”). To explain this statement, we need to consider small fluctuations 
of the field about the kink solution, similar to Eq. (1.5). We now have 


$ = x(x) + Wt, x) (1.13) 


where yx denotes the kink solution. The fluctuation field, w, obeys the linearized 
equation 


where 


X 


(1.11) 


dy — Iy +A- n )y = 0 (1.14) 
To find the fluctuation eigenmodes we set 
wae f(x) (1.15) 
where f(x) obeys 
-0f HAGA- 1°) = of (1.16) 


We will discuss all the solutions to this equation in Chapter 4. Here we focus on the 
translation mode. Since translations cost zero energy, there has to be an eigenmode 
with w = 0. This can be obtained by directly solving Eq. (1.16) or by noting that 
for small a, the solution in Eq. (1.12) can be Taylor expanded as 

dok 


Ox (x3 a) = ġk(x; a = 0) + a — (1.17) 
dx |a—0 
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Figure 1.2 The curve ranging from —1 to +1 as x goes from —oo to +00 shows 
the Z2 kink profile for à = 2 and n = 1. The energy density of the kink has also 
been plotted on the same graph for convenience, and to show that all the energy is 
localized in the narrow region where the field has a gradient. 


Comparing Eqs. (1.17) and (1.13), the zero mode solution is 


= 2 Bscch? E x (1.18) 
Lyn ND 2" l 


The solution in Eq. (1.9) can be used to calculate the energy density of the kink 


_ dbx 
fox) = Fe 


1 1 
a= AOS + 5 (eb) + V(¢x) 
= 0+ Vk) + V(dx) 


hn À 
= a soent ( ne) (1.19) 


where the second line is written to explicitly show that (3yo)? = 2V ($). The kink 
profile and the energy density are shown in Fig. 1.2. The total energy is 


2/2 m? 
E= faen r (1.20) 


As is apparent from the solution and also the energy density profile, the half- 
width of the kink is, 
o O S 


w = = (1.21) 
àn m My 


On the x > 0 side of the kink we have ¢ ~ +7 while on the x < 0 side we have 
gd ~ —n. At the center of the kink, ¢ = 0, and hence the Z) symmetry is restored 
in the core of the kink. Therefore the interior of the kink is a relic of the symmetric 
phase of the system. 


6 Classical kinks 


1.2 Rescaling 


It is convenient to rescale variables in the action in Eq. (1.2) as follows 
o=-, yh = VA nx” (1.22) 
Then the rescaled action is 
S=7 J dy 5a. — o — n| (1.23) 


where derivatives are now with respect to y”. The overall multiplicative factor, n?, 
does not enter the classical equations of motion. Hence the classical A* action is 
free of parameters.” 


1.3 Derrick’s argument 


In the context of rescaling, we now give Derrick’s result that there can be no static, 
finite energy solutions in scalar field theories in more than one spatial dimension 
[45]. Consider the general action in n spatial dimensions 


S= ae È X Ouo- von] (1.24) 


where the potential is assumed to satisfy V (°) > 0. The index on ° means that the 
model can contain an arbitrary number of scalar fields. Let a purported static, finite 
energy solution to the equations of motion be ¢5(x“) and consider the rescaled field 
configuration 


or") = polax") (1.25) 


where œ > 0 is the rescaling parameter. Then the energy of the rescaled field con- 
figuration is 


E[o%] = fos E + v(2)| (1.26) 


where the sum over a is implicit and V denotes the spatial gradient. Now define 
y” = ax" and this gives 


—n+2 
E[®5] = fes É 5 (VEO) ta” v (680) | (1.27) 


2 In quantum theory, however, the value of the action enters the path integral evaluation of the transition amplitudes 
and this will depend on 7. So the properties of the quantized kink also depend on the value of n (see Chapter 4). 
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Since the kinetic terms are non-negative, we find that with n > 2 and a > 1 this 
gives 


E[S] < Elo] (1.28) 


and hence ¢ cannot be an extremum of the energy. Only ifn = 1 can ġ§ be a static, 
finite energy solution. 

In more than one spatial dimension, Derrick’s argument allows for static solutions 
of infinite energy. The next section describes one such static solution in three spatial 
dimensions. 


1.4 Domain walls 


When kink solutions are placed in more than one spatial dimension, they become 
extended planar structures called “domain walls.” The field configuration for a Z2 
domain wall in the yz-plane in three spatial dimensions is 


plt, x,y,z) =n ianh( Snr) (1.29) 


The energy density of the wall is concentrated over all the yz-plane and is given 
by Eq. (1.19). The new aspects of domain walls are that they can be curved and 
deformations can propagate along them. These will be discussed in detail in Chap- 
ter 7. 

Another feature of the planar domain wall is that it is invariant under boosts in 
the plane parallel to the wall. This is simply because the solution is independent of 
t, y and z and any transformations of these coordinates do not affect the solution. 


1.5 Bogomolnyi method for Z, kink 


Rather than directly solve the equations of motion, as was done in Section 1.1, we 
can also obtain the kink solution by the clever method discovered by Bogomolnyi 
[20]. The method is to obtain a first-order differential equation by manipulating the 
energy functional into a “whole square” form 


1 2,1 2 
E= fa Fo + z0) + vo) 


2 
1 1 
z fa [300 + z9 F /2V($)) + VIO. 


(+00) 


1 1 
= J dx [300 +50 VVO j + f do’ /2V(6) 


(oo) 
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Then, for fixed values of @ at too, the energy is minimized if 


dæ =0 (1.30) 
and 
Ord = V2V(d) = 0. (1.31) 


Further, the minimum value of the energy is 


$(+00) 
Emin =f ap VVO). (1.32) 
ġ(=00) 

The energy can only be minimized provided a solution to Eq. (1.31) exists with 
the correct boundary conditions. This relates the choice of sign in Eq. (1.31) to 
the boundary conditions and to the sign in Eq. (1.32). In our case, for the Z2 
kink boundary conditions (¢(++00) > ¢(—0oo)), we take the — sign in Eq. (1.31). 
Inserting 


À 
VV) = [Fo — ¢’) (1.33) 


in Eq. (1.31) we get the kink solution in Eq. (1.9). 
The energy of the kink follows from Eq. (1.32) 


E= = —_ 1.34 
3 3 2x ( ) 
where m = VAn is the mass scale in the model (see Eq. (1.3)). 
1.6 Z, antikink 
In an identical manner, we can construct antikink solutions that have Q = —1. The 


boundary conditions necessary to get Q = —1 are ġ(+00) = #7 (see Eq. (1.8)). 
In the Bogomolnyi method, antikinks are obtained by taking the opposite choice of 
signs to the ones in the previous section 


1 1 
E= J dx Foo + 5 (a+ VVO) — VVD (1.35) 


This leads to the antikink solution 


k = —n ianh( Snr) (1.36) 
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1.7 Many kinks 


The kink solution is well-localized and so it should be possible to write down field 
configurations with many kinks. However, a peculiarity of the Z kink system is 
that a kink must necessarily be followed by an antikink since the asymptotic fields 
are restricted to lie in the vacuum: @ = +n. It is not possible to have neighboring 
Z> kinks or a system with topological charge |Q| > 1. 

There is a simple scheme, called the “product ansatz,” to write down approxi- 
mate multi-kink field configurations, i.e. alternating kinks and antikinks. Suppose 
we have kinks at locations x = k; and antikinks at x =/;, where i, j label the 
various kinks and antikinks. The locations are assumed to be consistent with the 
requirement that kinks and antikinks alternate: ...1; < k; < li+ı ... Then an ap- 
proximate field configuration that describes N kinks and N’ antikinks is given 
by the product of the solutions for the individual objects with a normalization 
factor 


1 N N’ 
O) = NNT [ae -ko | [i-e - 1) (1.37) 
i=1 j= 


j=l 


where œx is the kink solution. Note that |N — N’| < 1 since kinks and antikinks 
must alternate. 

The product ansatz is a good approximation as long as the kinks are separated 
by distances that are much larger than their widths. In that case, in the vicinity of a 
particular kink, say at x = k;, only the factor d(x — ki) is non-trivial. All the other 
factors in Eq. (1.37) multiply together to give +1. 

Another scheme to write down approximate multi-kink solutions is “additive” 
[109]. If ġ; denotes the ith kink (or antikink) in a sequence of N kinks and antikinks, 
we have 


N 
(x) =J gi EMN- Dn, N2 = N (mod 2) (1.38) 


i=l 


where the sign is + if the leftmost object is a kink and — if it is an antikink. 

Neither the product or the additive ansatz yields a multi-kink solution to the 
equations of motion. Instead they give field configurations that resemble several 
widely spaced kinks that have been patched together in a smooth way. If the multi- 
kink configuration given by either of the ansätze is evolved using the equation of 
motion, the kinks will start moving due to forces exerted by the other kinks. In the 
next section we discuss the inter-kink forces. 
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-a-R —a +R 


Figure 1.3 A widely separated kink-antikink. 


1.8 Inter-kink force 


Consider a kink at x = —a and an antikink at x = +a where the separation 2a is 
much larger than the kink width (see Fig. 1.3). We would like to evaluate the force 
on the kink owing to the antikink [109]. 

The energy-momentum tensor for the action Eq. (1.2) with a general potential 


V($)is 
1 
Tav = 3 Pdh — Suv E - vo (1.39) 


where g,,, is the metric tensor that we take to be the flat metric, that is, guv = 
diag(1, —1). The force exerted on a kink is given by Newton’s second law, by 
the rate of change of its momentum. The momentum of a kink can be found by 
integrating the kink’s momentum density, T” = —Tp;, in a large region around the 
kink. If the kink is located at x = —a, let us choose to look at the momentum, P, 
of the field in the region (—a — R, —a + R) 


—a+R 

P= -f dx 0,60. (1.40) 
—a—R 

After using the field equation of motion (for a general potential) and on performing 

the integration, the force on the field in this region is 


dP 1 x 1 > —a+R 
F=—= [-30) = Odor vo) ray) 


7 dt —a—R 


To proceed further we need to know the field ¢ in the interval (—a — R, —a + R). 
This may be obtained using the additive ansatz given in Eq. (1.38) which we take 
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as an initial condition 


P(t = 0, x) = ba) + pr) — pkl) (1.42) 


In addition, we assume that the kinks are initially at rest 


jol 50 (1.43) 
t=0 


The expression for the force is further simplified by using the Bogomolnyi equation 
(Eq. (1.31)) which is satisfied by both ¢, and 6 


(ag) = 2V(¢) (1.44) 
This gives 


—a+R 


F= EI + V (k + $k — Px(O)) — V (Qx) — vē] (1.45) 


—a—R 


The terms involving the potential can be expanded since the field is nearly in the 
vacuum at x = —a + R. Let us define 


dj =o(-—a+ R), pk = fka R) 
Apk = ox(—a £ R) — pk(+00) 
Ady = $x(~a R) — $x(—00) (1.46) 


(Note that the argument in the very last term is —oo, independent of the signs in 
the other terms. This is because both x = —a + R lie to the left of the antikink.) 
Also define 


mi, = V"(px(00)) = V"(Gx(00)) (1.47) 


Then the force is 
F = (Ospi bf — OG, PL) + M (Ap Ad! — Ad, Ady) (148) 


Let us illustrate this formula for the Z2 kink, where 


Px(x) = n tanh(o(x + a)) 
k(x) = —n tanh(o(x — a)) (1.49) 


with o = ./A/2 n. Inserting these expressions in Eq. (1.48) and retaining only the 
leading order behavior gives 


4m4 
F= a oo (1.50) 
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where / = 2a is the kink separation. The force is attractive since it is acting on the 
kink at x = —a and points toward the antikink at x = +a. 

The result for the force could have been guessed from other considerations. The 
kinks are interacting by the exchange of massive scalars of mass my. As described 
in many quantum field theory texts [119] the force mediated by scalar interactions 
is the Yukawa force which goes like exp(—m,,/). The dimensionful prefactor of 
the force can be deduced on dimensional grounds while the numerical coefficient 
requires more detailed analysis. 


1.9 Sine-Gordon kink 


The sine-Gordon model is a scalar field theory in 1 + 1 space-time dimensions 
given by the Lagrangian 


1 2 Q 
L= 59) -g (1 — cos(B¢)) (1.51) 


The model is invariant under @ — ¢ + 27n where n is any integer and thus pos- 
sesses Z symmetry. The vacua are given by ¢ = 27rn/B and are labeled by the 
integer n. 

As in the Z3 case, the classical kink solutions can be found directly from the 
second-order equations of motion or by using the Bogomolnyi method (see Sec- 
tion 1.5). The kinks are solutions that interpolate between neighboring vacua. The 
unit charge kink solution is 


dk = S an! (e~) + o(—0oo) (1.52) 


where the inverse tangent is taken to lie in the interval (—2r/2, +7 /2). The antikink 
with @(—0oo) = 27/6 and ¢(+00) = 0 is obtained from Eq. (1.52) by replacing x 
by —x. 


oy = tan! (e77) + (+00) (1.53) 


The width of the kink follows directly from these solutions and is ~1/,/a. 
The energy of the kink also follows from Bogomolnyi’s method (Eq. (1.32)) 


8 Ja 
Exc = B 


Defining my = ./a — the mass of excitations of the true vacuum — and Ji = fap 
we get 


(1.54) 


(1.55) 
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While the Z% and sine-Gordon kinks are similar as classical solutions, there 
are some notable differences. For example, it is possible to have consecutive sine- 
Gordon kinks whereas in the Z3 case, kinks can only neighbor antikinks. In addition, 
the sine-Gordon system allows non-dissipative classical bound states of kink and 
antikink — the so-called “breather” solutions — while no such solutions are known 
in the Z2 case (though see Section 3.1). The sine-Gordon kink is also much more 
amenable to a quantum analysis as we discuss in Chapter 4. 

We can use the additive ansatz described in Section 1.7 to construct field configu- 
rations for many kinks. Specializing to a kink-kink pair (¢(—co) = 0 to (+00) = 
4x /ß) and a kink-antikink pair (@(—0oo) = 0 and back to ¢(+00) = 0), we have 


Ox(t, x) = SEn (ev@-9)) + tan! (evae—)] (1.56) 


A(t, xX) = Sien (evee-)) + tan! (ever)| — = (1.57) 


with b > a. 

The additive ansatz described above gives approximate solutions to the equations 
of motion for widely separated kinks (b — a >> 1/,/a). A one-parameter family 
of exact, non-dissipative, breather solutions composed of a kink and an antikink is 


dt. x)= tan! || (1.58) 
where n = Va —@/q and the tan~! function is taken to lie in the range 
(—/2, +7/2). The frequency of oscillation, œw, is the parameter that labels the 
different breathers of the one-parameter family. 

To see the breather as a bound state of a kink and an antikink, note that 
p(t, too) = 0. Also, if n > 1, then ¢@(¢,0) + 27/6 during the time when 
n sin(@t) >> 1. Hence the breather splits up into a kink and an antikink for part 
of the oscillation period. For the remainder of the oscillation period, the kink and 
an antikink overlap and a clear separation cannot be made. 

The constant energy of the breather is evaluated by substituting the solution at 
t = 0 (for convenience) in the sine-Gordon Hamiltonian with the result 


16 —— | f 
Ey = B2 6 Aa w2 = 2EG 1 -= 24 (1.59) 
a 


As expected, when w — 0, the breather energy is twice the kink energy. 
As in Section 1.8 we can find the force on a kink owing to an antikink: from 
Eq. (1.48) the leading order behavior of the force is 
20m? a 


where / is the kink separation. 
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1.10 Topology: zo 


The kinks in the Zz and sine-Gordon models can be viewed as arising purely 
for topological reasons, as we now explain. A very important advantage of the 
topological viewpoint is that it is generalizable to a wide variety of models and 
can be used to classify a large set of solutions. When applied to field theories in 
higher spatial dimensions, topological considerations are convenient in order to 
demonstrate the existence of solutions such as strings and monopoles. 

Consider a field theory for a set of fields denoted by ® that is invariant under 
transformations belonging to a group G. This means that the Hamiltonian of the 
theory is invariant under G: 


HIS] = H[®*] (1.61) 


where g € G and ® represents ® after it has been transformed by the action of 
g. The group G is a symmetry of the system, if Eq. (1.61) holds for every g € G 
and for every possible ®. Now, let the Hamiltonian be minimized when ® = ®o. 
Then, from Eq. (1.61), it is also minimized with ® = i for any g € G, and the 
manifold of lowest energy states — “vacuum manifold” — is labeled by the set of 
field configurations $. However, there will exist a subgroup (sometimes trivial), 
H of G, whose elements do not move ®o: 


P! = By (1.62) 


Hence, a group element gh € G acting on ®o has the same result as g acting on 
o (because h acts first and does not move ®o). So, while the configuration 25 
has the same energy as ®ọ for any g, not all the fs are distinct from each other. 
The distinct és are labeled by the set of elements {gh : h € H} = gH. The set 
of elements {gH : g € G} are said to form a “coset space” and the set is denoted 
by G/H; each element of the space is a coset (more precisely a “left coset” since 
g multiplies H from the left). Therefore the vacuum manifold is isomorphic to the 
coset space G/H. 

We have so far connected the symmetries of the model to the vacuum manifold. 
Now we discuss the tools for describing the topology of the vacuum manifold. This 
will lead to a description of the topology of the vacuum manifold directly in terms 
of the symmetries of the model. 

The topology of a manifold, M, is classified by the homotopy groups, ,,(M; xo), 
n = 0, 1,2,...The idea is to consider maps from n-spheres to M, with the 
image of an n-sphere in M containing one common base point, xo (see Fig. 1.4). 
If two maps can be continuously deformed into each other, they are considered 
to be topologically equivalent. In this way, the set of maps is divided into equiva- 
lence classes of maps, where each equivalence class contains the set of maps that 
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are continuously deformable into each other. The elements of z,(M; xo) are the 
equivalence classes of maps from $” to M with fixed base point. It is also possible 
to define (except for n = 0 as explained below) a suitable “product” of two maps: 
essentially the product of maps f and g (denoted by g - f) and is defined to be 
“f composed with g” or “ f followed by g.” Then it is easily verified that the prod- 
uct is closed, associative, an identity map exists, and every map has an inverse. In 
mathematical language, Y f, g, h € G, 


f-gEeG 

f(g -A=(F-g)-h 

Je e G suchthat f-e=e-f=f 

3 f! e G such that f.f '=f'-f=e (1.63) 


Thus all the group properties are satisfied and 7,(M; xo) is a group. 

Two homotopy groups with different base points, say 2,(M; xo) and z,(M; x), 
can be shown to be isomorphic and hence the reference to the base point is often 
dropped and the homotopy group simply written as z,,(7). Mathematicians have 
calculated the homotopy groups for a wide variety of manifolds and this makes it 
very convenient to determine if a given symmetry breaking leads to a topologically 
non-trivial vacuum manifold [145, 3, 171]. 

In the case of kinks or domain walls, the field defines a mapping from the 
points x = +00 to the vacuum manifold. Hence the relevant homotopy “group” is 
xo(M; xo), which contains maps from s? (a point) to M. Since the base pointis fixed, 
the image of either of the two possible $s (x = +00) has to be xo, and mo(M; xo) 
is trivial. Even if we do not impose the restriction that the maps should have a fixed 
base point, it is not possible to define a suitable composition of maps. Therefore zo 
does not have the right group structure and should merely be considered as a set of 
maps from S° to the vacuum manifold. The exception occurs if M = G/H is itself 
a group, which occurs when H is a normal subgroup of G, because then zo(M) can 
inherit the group structure of M. In this case, the product of two maps from S° to 
M can be defined to be the map from S° to the product of the two image points in 
M. Generally, however, 79(M) should simply be thought of as a set of maps from 
S’ to the various disconnected pieces of M. 

To connect the elements of the homotopy groups to topological field config- 
urations assume that the field, ®, is in the vacuum manifold on S%.. Therefore, 
P~ = P(x € Sh) defines a map from S” to the vacuum manifold and this map can 
be topologically non-trivial if 7„(M) is non-trivial. We want to show that if the map 
~ is topologically non-trivial, ® cannot be in the vacuum manifold at all points 
in the interior of S}. Consider what happens as the radius of S% is continuously 
decreased. If the field remains on the vacuum manifold, continuity implies that the 
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S” M 


Figure 1.4 The nth homotopy group consists of maps from the n-dimensional 
sphere, S”, to the vacuum manifold, M, such that the image of any map contains 
one common base point xọ € M. If two maps can be continuously deformed into 
each other, they are identified, and correspond to the same element of z,. If two 
maps cannot be continuously deformed into each other, then they correspond to 
distinct elements of zr,,. For example, this can happen if one of the maps encloses 
a “hole” in M, while the other encloses the hole a different number of times. 


map ®p from a sphere of radius R to M must also be non-trivial. Then as R — 0, 
the map would still be non-trivial, implying that the field is multivalued at the origin 
since the field must continue to map Sẹ} non-trivially as R — 0. However, a field 
(by definition) cannot be multivalued. The only way out is if the field does not lie on 
the vacuum manifold everywhere. Therefore non-trivial topology at infinity implies 
that the energy density does not vanish at some points in space. The distribution of 
energy density is the topological defect which, depending on dimensionality, can 
manifest itself as a domain wall (n = 0), string (n = 1) or monopole (n = 2) or 
texture (n = 3). 

The above argument establishes that topologically non-trivial boundary condi- 
tions imply non-vanishing energy in the field. However, it does not establish that 
a Static solution exists with those boundary conditions. These must be found on a 
case-by-case basis. Indeed there are examples of topologically non-trivial boundary 
conditions where no static solution exists.* Also distinct elements of z,(M) (n > 1) 
need not lead to distinct field solutions. Only those solutions that correspond to el- 
ements of z,,(M/) that cannot be continuously deformed into each other, if the maps 
are released from the base point, are distinct. The italicized remark is in recognition 
of the fact that there can be two maps that are mathematically distinct (i.e. cannot be 
deformed into each other) only because they are fixed at the base point. However, 
the analog of a “common base point” in field theory would be to restrict attention 
to field configurations for which the fields attain a certain fixed value at some point 
on S% . Such a restriction is generally unphysical and hence, we are interested in 


3 For example, in three dimensions, the boundary conditions corresponding to a charge two ’t Hooft-Polyakov 
magnetic monopole [124, 79] do not lead to any solution (for all but one value of model parameters). This is 
because any field configuration with those boundary conditions breaks up into two magnetic monopoles, each 
of unit charge, that repel each other and are never static. 
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maps that cannot be deformed into each other even if we release the restriction that 
all maps have a common base point (for a more detailed discussion, see [171]). 

In the case when the vacuum manifold has disconnected components, mo(G/H) 
is non-trivial since there are points (zero-dimensional spheres) that lie in different 
components that cannot be continuously deformed into one another. Therefore kinks 
occur whenever zo(G/#) is non-trivial. In the A®* model, G = Z2, H = 1 and 
xo(G/H) = Z2. In the sine-Gordon model G = Z, H = 1 and xo(G/H) = Z. If 
To = Zy, we name the resulting kinks “Zy kinks.” In these simple examples, zo 
forms a group because G is Abelian and so G/H itself is a group. An example 
in which zo is not a group can be constructed by choosing G = S3 (S, is the 
permutation group of n elements) broken down to H = Sp. 

The kinks in a model with disconnected elements in M can now be classified. 
Every element of z9(M) corresponds to a mapping from a point at spatial infinity 
to M and hence specifies a domain at infinity. Kinks occur if the domains at +00 
are distinct. Therefore pairs of elements of xo(M) classify domain walls. 


1.11 Bogomolnyi method revisited 


The Bogomolnyi method can be extended to include a large class of systems. Let us 
start with the general energy functional for a matrix-valued complex scalar field ® 


E= J dx[Tr|d,®|? + Trla D? + V(®, &*)] 
= J dx[ð, 0,3 Poa + ID dk Ory + V(O, DH] (164) 


where a sum over matrix components labeled by a, b is understood. As in Sec- 
tion 1.5, we would like to write the energy density in “whole square” form 


paf dx[Tr{13,®|? + |d, F U(®)|? £ (a, @'U) + (UTA, &)}] (1.65) 


where we are restricting ourselves to static solutions and U is some matrix-valued 
function of ® such that 


Tr(U'U) = V(®, &*) (1.66) 
The energy is minimized if 
0,2 = 0 (1.67) 
and 


Tr|d,® F U(®, &*)|* = 0 (1.68) 
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which in turn gives 


0,® F U(®, &*) = 0 (1.69) 
The energy of the kink is 
+00 
E= «| dx Tr(a, @' U + U'I D) (1.70) 
—oo 


There is a further special case — the “supersymmetric” case — in which the energy 
integral can be performed explicitly. This is if U is a total derivative 


ow 
U* = — 1.71 
PF (1.71) 
where W(®, ®*) is the “superpotential,” assumed to be real. Then 
ToS ow ow 
E= +f dx Tr | 4, @! —— + 0,67 — 
25 ad* oP 
+00 
= «| dx 0,W 
—oo 
= £[W(®(+00)) — W(®(—0o))] (1.72) 


Therefore we see that the Bogomolnyi method allows for first-order equations 
of motion provided that V can be written as Tr(U'U). The method also provides 
an explicit expression for the kink energy if V is given in terms of a superpotential 
W as 
oT ee eee Le i 

(®) = Tr(U'U) = Tr 


(1.73) 


1.12 On more techniques 


The kink solutions we have been discussing fall under the more general category 
of “solitary waves,” often discussed under the soliton heading. Strictly speaking, 
for a solution to classify as a “soliton,” it also has to satisfy certain conditions 
on its scattering with other solitons. The subject is incredibly rich, and has led to 
the development of very sophisticated mathematical techniques such as Backlund 
transformations, inverse scattering methods, Lax heirarchy, etc. In addition, solitons 
have found tremendous importance in physical applications, especially non-linear 
optics and communication. Readers interested in the mathematics and physics of 
solitons might wish to consult [1, 48, 56]. 

Strict solitons are usually discussed in one spatial dimension and have limited 
application in the context of particle physics. Nonetheless, there are equally so- 
phisticated techniques to study solitary wave solutions in higher dimensions. In 
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particular, the ADHM construction [12] is used to find instanton solutions in four 
spatial dimensions and the Nahm equations lead to magnetic monopole solutions 
in three dimensions [114]. 

The soliton analyses mentioned above consider equations with complicated non- 
linear terms and higher derivatives. In the context of particle physics, such terms 
and derivatives are rarely encountered. However, one complication that arises is 
due to larger (non-Abelian) symmetry groups. In the next chapter we will take 
the analysis of this section to such particle-physics motivated models. There we 
will find a spectrum of kink solutions with unusual interactions. As we proceed to 
further chapters, we will learn that the physics of such non-Abelian kinks can be 
quite different from that of the simple kinks discussed in this chapter. 


1.13 Open questions 


1. Discuss the conditions needed for a breather solution to exist. If an exact breather does not 
exist, can there be an approximate breather (see Section 3.1)? What is the approximation? 


2 


Kinks in more complicated models 


The Z, and sine-Gordon kinks discussed in the last chapter are not representative of 
kinks in models where non-Abelian symmetries are present. Kinks in such models 
have more degrees of freedom and this introduces degeneracies when imposing 
boundary conditions, leading to many kink solutions with different internal struc- 
tures (but the same topology). Indeed, kink-like solutions may exist even when the 
topological charge is zero. The interactions of kinks in these more complicated 
models, their formation and evolution, plus their interactions with other particles 
are very distinct from the kinks of the last chapter. 

We choose to focus on kinks in a model that is an example relevant to particle 
physics and cosmology. The model is the first of many Grand Unified Theories of 
particle physics that have been proposed [63]. The idea behind grand unification 
is that Nature really has only one gauge-coupling constant at high energies, and 
that the disparate values of the strong, weak, and electromagnetic coupling con- 
stants observed today are due to symmetry breaking and the renormalization-group 
running of coupling constants down to low energies. Since there is only one gauge- 
coupling constant in these models, there is a simple grand unified symmetry group 
G that is valid at high energies, for example, at the high temperatures present in 
the very early universe. At lower energies, G is spontaneously broken in stages, 
eventually leaving only the presently known quantum chromo dynamics (QCD) and 
electromagnetic symmetries SU (3)e x U(1)em of particle physics, with its two dif- 
ferent coupling constants. It can be shown [63] that the minimal possibility for G is 
SU(5). However, since Grand Unified Theories predict proton decay, experimental 
observation of the longevity of the proton (~ 5 x 10°? years) leads to constraints 
on grand unified models. The (non-supersymmetric) SU(5) Grand Unified Theory 
is ruled out by the current lower limits on the proton’s lifetime. Therefore particle- 
physics model builders consider yet larger groups G, or with an extended scalar 
field sector, or supersymmetric extensions of SU(5), and other models based on 
larger groups. Even if the symmetry group is larger than SU (5), it often happens 
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that after a series of symmetry breaking, the residual symmetry is SU(5), which 
then proceeds to break to the current symmetry group. Hence the study of SU (5) 
symmetry breaking is extremely relevant to particle physics, even if it is not the 
ultimate grand unified symmetry group. 

In this chapter we shall study kinks in a model with SU(5) x Z2 symmetry 
though almost all the discussion can be generalized to an SU(N) x Z2 model for 
odd values of N [163, 120]. The extra Z2 symmetry is explained in the next section. 
Since we only desire to study kinks in a particle-physics motivated model, it would 
seem simpler to choose a model based on the smaller SU(3) group. However, it 
can be shown that there is no way to construct a model with just SU (3) symmetry 
and with the simplest choice of field content, which is one adjoint field. Instead, 
the model must have the larger O(8) symmetry. Other fields need to be included 
so as to reduce the O(8) to SU(3), but that introduces additional parameters which 
make the $U(3) model more messy than the SU(5) model. 

Dealing with continuous groups such as SU(5) requires certain background 
material. The fundamental representation of SU(N) generators is described in Ap- 
pendix B. A summary of some aspects of the SU(5) model of grand unification is 
given in Section 5.5. 


2.1 SU(5) model 


The SU(5) model can be written as! 
1 
L = Tr(D, 0) — zX wX) — V(®) (2.1) 


where, in terms of components, ® is a scalar field (also called a Higgs field) 
transforming in the adjoint representation of SU (5), that is, > ®' = g@g" for 
g € SU(5). The gauge field strengths are Xu» = X/,T“ and the SU(5) generators 
T° are normalized such that Tr(7¢T”) = 5% /2. The definition of the covariant 
derivative is 


D, = 0, —ieX, (2.2) 
and its action on the adjoint scalar is given by 
D,® = 0,8 — ie[X,,, ®] (2.3) 
The gauge field strength is given in terms of the covariant derivative via 
—ieX,,, = [D,,, Dy] (2.4) 


1 We are using the Einstein summation convention in which repeated group and space-time indices are summed 
over. So, explicitly, ® = yi PIT”. See Appendix B for more details on the SU (5) generators T“. 
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and the potential is the most general quartic in ® 
V(®) = —m?Tr(@*) + hT] + ATr(@*) + yTr(?)— Vo (2.5) 


where Vo is a constant that is chosen so as to set the minimum value of the potential 
to zero. 

The model in Eq. (2.1) does not have any topological kinks because there are 
no broken discrete symmetries. In particular, the Z2 symmetry under P + —® is 
absent owing to the cubic term in Eq. (2.5). Note that P — —® is not achievable 
by an SU (5) transformation. To show this, consider Tr(?). This is invariant under 
any SU(5) transformation, but not under P — —®. However, if y = 0, there are 
topological kinks connecting the two vacua related by ® — —®. For non-zero but 
small y, these kinks are almost topological. In our analysis in this chapter we set 
y = 0, in which case the symmetry of the model is SU(5) x Z2. The philosophy 
underlying grand unification does not forbid discrete symmetry factors since such 
factors do not entail additional gauge-coupling constants. Indeed, model builders 
often set y = 0 for simplicity. Now a non-zero vacuum expectation value of ® 
breaks the discrete Z factor leading to topological kinks. 


2.2 SU(5) x Z2 symmetry breaking and topological kinks 
The potential in Eq. (2.5) has a (degenerate) global minimum at 
n 


2/15 
where n = m/4/2’ provided 


o = diag(2, 2, 2, —3, —3) (2.6) 


7 
A> 0, N=h+—Aa>0 D7. 
> ta (2.7) 
For the global minimum to have V(®g) = 0, in Eq. (2.5) we set 
Ae 
n=-7" (2.8) 


As discussed in Section 1.10, if we transform ®o by any element of SU (5) x Zo, 
the transformed ®ọ is still at a minimum of the potential. However, ®o is left 
unmoved by transformations belonging to 

[SU(3) x SUQ) x U(1)] 
G31 = (2.9) 
Z3 X Z2 
where SU (3) acts on the upper-left 3 x 3 block of Po, SU (2) on the lower-right 
2 x 2 block, and U(1) is generated by ®o itself. Hence, G32; is the unbroken 
symmetry group. 
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Figure 2.1 The vacuum manifold of the SU(5) x Z) model consists of two dis- 
connected 12-dimensional copies. Kink solutions correspond to paths that originate 
in one piece at x = —oo, denoted by ®(—), leave the vacuum manifold, and end 
in the other disconnected piece at x = +00. Topological considerations specify 
that (+) has to lie in the disconnected piece on the right, but not where it should 
be located within this piece. 


SU(5) has 24 generators while the unbroken group, G32;, has a total of 12 
generators, namely, 8 of SU(3), 3 of SU(2), and 1 of U(1). Therefore the vacuum 
manifold is 24 — 12 = 12 dimensional but in two disconnected pieces as depicted in 
Fig. 2.1 because of the Z3 factor. Kink solutions occur if the boundary conditions lie 
in different disconnected pieces. However, if we start at some point on the vacuum 
manifold at x = —oo, say ®(—0o0) = _, we have a choice of boundary conditions 
for ®;, the vacuum expectation value of ® at x = +00 (compare with the Z case 
where the path had to go from definite initial to definite final values of ®). 

We will narrow down the possible choices for ® very shortly. First we point out 
that the gauge fields can be set to zero in finding kink solutions [163]. To see this 
explicitly, the only linear term in the gauge field is ieTr(X‘[®, 0; ©]). However, our 
solution for ® satisfies [®, d;®] = 0 [120] and so the variation vanishes to linear 
order in gauge field fluctuations. A closer look also reveals that the quadratic terms 
of perturbations in the gauge fields contribute positively to the energy of the kink 
solutions and so the gauge fields do not cause an instability of the solutions [163]. 
Hence we set 


X, =0 (2.10) 


As we now show, the boundary conditions that lead to static solutions of the equa- 
tions of motion are rather special [120]. 


Theorem: A static solution can exist only if [®,, ®_] = 0. 


We only give a sketch of the proof here since it is of a technical nature. The 
essential idea is that if ©, (x) is a static solution, then the energy should be extrem- 
ized by it. By considering perturbations of the kind U(x)®,U'(x) where U(x) is an 
infinitesimal rotation of SU (5), one finds that the energy can be extremized only if 
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[®,, 39x Pk] = 0 for all x. Now at large x, we have Pk —> ® 4. In this region 0, ®, 
has terms that are proportional to ®_ as well, even if these are exponentially small, 
since ®(x) is an analytic function. Hence, a static solution requires [P,, ®_] = 0. 

The theorem immediately narrows down the possibilities that we need to consider 
when trying to construct kink solutions. If we fix 


n 
2/15 


®, can take on the following three values 


_ = Po = 


diag(2, 2, 2, —3, —3) (2.11) 


G0 = —-— fas. Oo. —3, —3 
reen A i 


pP = —— diag, 2, —3, 2, —3) 
24 we S 


Pa _dias(2, —3, —3, 2,2) (2.12) 
ETRE 
One can also rotate these three choices by elements of the unbroken group G321— 
that leaves ®_ invariant and obtain three disjoint classes of possible values of ®,. 
The three choices given above are representatives of their classes. 
The kink solution for any of the three boundary conditions is of the form 


pP = FOME + FO CM? + gO Q)M (2.13) 
where q = 0, 1, 2 labels the solution class, 
po ah po” eo” po” 
(q) + = (4) + — 
M{ = 7 MY’ = = (2.14) 


and M“ will be specified below. 
The boundary conditions for Fg are 


F900) =F1, FO Goo)=41, goo) =0 (215) 


The formulae for M@ and M® can now be explicitly written using Eq. (2.12) 
in (2.14) 


5 
(4) : 
My =n 7 J804 Las =l, 02-4) (2.16) 
1 
M” = n——diag(—413_,, 1,, 14, 615_,) (2.17) 
"1/55 8 3—q> 4q, 4q, 942-4 


M® = pdiag(q(2 — q)13-4, —(2 — 4)G — q)124, 9B — q)12-4) (2.18) 


with the normalization u given by 


u = nl2q(2 — p8 -g2 — 5q)? (2.19) 
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Figure 2.2 The profile functions FO (x) (nearly 1 throughout), Fo (x) (shaped 
like a tanh function), and g(x) (nearly zero) for the q = 1 topological kink with 
parameters h = —3/70, à = 1, and ņ = 1. 


If g = 0 org = 2 we set u = 0. We have used 0; and 1; to denote the k x k zero 
and unit matrices respectively. Note that the matrices M” are relatively orthogonal 


TrM?M) = 0 (2.20) 


but are not normalized to 7/2. 

Now we discuss the three kink solutions in the SU(5) x Z2 model. For q = 0, 
the solution is that of a Z2 kink that has been embedded in the SU(5) x Z2 model. 
The explicit solution is 


FO(x) =0, FOX) = -tann(*), g(x) =0 (2.21) 


where w = /2/m.Forq = 1, the profile functions have been evaluated numerically 
and are shown in Fig. 2.2. Approximate analytic solutions can also be found in 
[120]. For g = 2 the solution has also been found numerically. Here we describe 
an approximate solution which is exact if 
h 3 

= 2.22 
À 20 pen 
i.e. A’ = A/12. With this particular choice 


FO (x) = 1, Fa) = tanh), g(x) =0 (2.23) 
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where w = /2/m. This is also an approximate solution for h/A ~ —3/20. The 
energy of the approximate solution can be used to estimate the mass of the g = 2 


kink 

M (1 52 1/2 JP 

M? x 1+ = M® (2.24) 
6 l6 12! 6 


where M® denotes the mass of the q = 2 kink, and M = 2./2m3/32'. The 
expression for the energy is exact for h/A = —3/20. 

It can be shown for a range of parameters that the g = 2 kink solution is per- 
turbatively stable. Numerical evaluations of the energy find that the g = 2 kink is 
lighter than the q = 0, 1 kinks for all values of p. Equation (2.24) shows the q = 2 
kink is lighter than the q = 0 kink for a large range of parameters. This can be 
understood qualitatively by noting that only one component of ® changes sign in 
the q = 2 kink, while 3 and 5 components change sign in the q = 1 and q = 0 
kinks respectively. 


2.3 Non-topological SU(5) x Z2 kinks 


An interesting point to note is that the ansatz in Eq. (2.13) is valid even if pP are 
not in distinct topological sectors. These imply the existence of non-topological 
kink solutions in the model [120]. If we include a subscript NT to denote “non- 
topological” and T to denote “topological,” we have 


Ont, = FO EMG, + FOM + gO 0M (2.25) 


where the Mnrs matrices are still defined by Eq. (2.14) with the non-topological 
values of P+. Myris still given by Eq. (2.18). To consider a non-topological domain 
wall, we simply want to consider ®, to be in the same discrete sector as ®_. If r+ 
denotes a boundary condition for a topological kink, a possible boundary condition 


for a non-topological kink is: Pyr+ = —P7,. Then we find 
M% =M, MY =My), MQ = My? (2.26) 
Hence 
ON = FO OME + RP OME + ¢ OMY? (2.27) 


To get F® for the non-topological kink we have to solve the topological FY 
equation of motion but with the boundary conditions for Fa ) (see Eq. (2.15)). To 
obtain g for the non-topological kink, we need to interchange F se and F in the 
topological equation of motion. The boundary conditions for g are unchanged. 
Generally the non-topological solutions, when they exist, are unstable. However, 
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Table 2.1 The space of three topological kinks in the SU (5) model. 


G32; is the group SU(3) x SU(2) x U(1). The dimensionality of the space 
of each type of kink is also given. 


Kink Space Dimensionality 
qg=0 G321/G321 0 
q=1 G321/[SU(2) x U(1)?] 6 
q=2 G31 /[SU(2} x U(1)°] 4 


the possibility that some of them may be locally stable for certain potentials cannot 
be excluded. 


2.4 Space of SU(5) x Z2 kinks 


The kink solutions discussed in Section 2.1 can be transformed into other degenerate 
solutions using the SU(5) transformations. Hence, each solution is representative 
of a space of solutions. We now discuss the space associated with each of these 
solutions. 

If we denote a kink solution in the SU(5) x Z2 model by oe. another solution is 


go" = nent, hE Gz- (2.28) 


where G32;_ is the unbroken group whose elements leave ®_ unchanged.” The 
reason oom also describes a solution is that the rotation A does not change the 
energy of the field configuration, po, Therefore po” has the same energy and 
the same topology as © and hence it describes another kink solution. 

Of the elements of G321—, there are some that act trivially on po and for these 
h, po” is not distinct from pP. These elements form a subgroup of G321— that we 
call K4. Therefore the space of kinks can be labeled by elements of the coset space 
G321-/ K4. Since we are given the forms of the kink solutions in Eq. (2.13), it is not 
hard to work out K,. For example, for the q = 2 kink, K; is given by the SU(5) 
elements that commute with both G32;_ and G321+ and so K, = SU (2)? x UAP. 
Once we have determined K, the dimensionality of the coset space G32;_/Kg is 
determined as the dimensionality of G32;_, which is 12, minus the dimensionality 
of K,, which is 12, 6, and 8 for q = 0, 1, and 2 respectively. 

The three classes of kink solutions labeled by the index q in the SU(5) x Z2 
model have different spaces as shown in Table 2.1. 


2 We could also have included elements that change po as well as ®_. These would simply be global rotations 
of the entire solution and would be the same for every type of defect. 
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The dimensionality of the space of a given type of kink solution also corresponds 
to the dimensionality of the space of boundary conditions ®, for which that type 
of kink solution is obtained. As an example, there is only one value of P, namely 
P, = —@_, that gives rise to the q = 0 kink. While for the q = 1 kink, one can 
choose ®, to be any value from a 6-dimensional space. This means that, in any 
process where boundary conditions are chosen at random, the probabilities of get- 
ting the correct boundary conditions for a q = 0 or a q = 2 kink are of measure 
zero, since the space of boundary conditions for the q = 1 kink is two dimensions 
greater than that for the q = 2 kink. In any random process, the q = 1 kink is always 
obtained. Since this kink is unstable, it then decays into the g = 2 kink. Therefore 
the production of g = 2 kinks is a two-step process in this system. We will see 
further evidence of this two-step process in Chapter 6. 


2.5 S, kinks 


The SU(5) x Z2 model discussed above shows novel features because of the large 
non-Abelian symmetry. It is possible to see some of the richness of the model 
by going to a simpler model where the continuous non-Abelian symmetries are 
replaced by discrete non-Abelian symmetries (also see [92] for a similar model). 
If we truncate the SU(5) x Z2 model to just the diagonal degrees of freedom 
of ®, we get a model that is symmetric only under permutations of the diag- 
onal entries and the overall Z2. Hence the symmetry group is Ss; x Z2, where 
Ss is the permutation group of five objects. The model now has four real scalar 
fields, one for each diagonal generator of SU(5). With this truncation we can 
write 


® — fida + forg + fat3 + faY (2.29) 


where the f; are functions of space and time, and the generators A3, Ag, T3, and Y are 
defined in Appendix B. Inserting this form of ® into the SU(5) x Z2 Lagrangian 
in Eq. (2.1) we get 


1 4 
= 5 D OS +V Gi fo fr fa) (2.30) 
i=l 
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fe m2 
[AE + S3) +985 | fe + R (#- 2) + q" (2.31) 
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Figure 2.3 The vacuum manifold for the $5; x Z2 model contains two sets of ten 
points related by the Z} symmetry. Kink solutions exist that interpolate between 
vacua related by Z2 transformations and also between vacua within one set of ten 
points. The former correspond to the topological kinks in SU(5) x Z2 and the 
latter to the non-topological kinks in that model. 


This model has the desired $5 x Z2 symmetry because it is invariant under permu- 
tations of the diagonal elements of ®, that is, under permutations of various linear 
combinations of f;. The Z2 symmetry is under f; — — f; for every i. 

Symmetry breaking proceeds as in the SU(5) x Zz case. The Ss x Z2 symmetry 
is broken by a vacuum expectation value along the Y direction i.e. f4 4 0. The 
residual symmetry group consists of permutations in the SU (3) and SU (2) blocks. 
Therefore the unbroken symmetry group is H = S3 x S2. There are 5! x 2 = 240 
elements of S5 x Z2 and 3! x 2! = 12 elements of H. Therefore the vacuum mani- 
fold consists of 240/12 = 20 distinct points. Ten of these points are related to 
the other ten by the non-trivial element of Zz as shown in Fig. 2.3. If we fix the 
boundary condition at x = —oo, then a Zz kink can be obtained with ten different 
boundary conditions at x = +00. These ten solutions must somehow correspond 
to the kink solutions that we have already found in the SU(5) x Zz case. Counting 
all the possible different diagonal possibilities for P, in the SU(5) x Z2 model we 
see that there are three q = 2 kinks, six q = 1 kinks, and one q = 0 kink, making 
a total of ten kinks. In the Ss; x Z2 model there are ten more (one of these is the 
trivial solution) kinks that do not involve the Zz transformation (change of sign) 
in going from ®_ to ®,. These are the ten remnants of the non-topological kinks 
described in Section 2.3. 


2.6 Symmetries within kinks 


The symmetry groups outside the kink, G321+, are isomorphic (see Fig. 2.4). How- 
ever, the fields transform differently under the elements of these groups. As a result, 
there is a “clash of symmetries” [43] inside the kink, and the unbroken symmetry 
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Kink 


Figure 2.4 A kink and the symmetries outside denoted by H+. The groups H, 
and H_ are isomorphic but their action on fields may not necessarily be identical. 


group within the kink is generally smaller than that outside. This does not happen 
in the case of the Z2 kink in which the symmetry outside is trivial while inside it is 
Z> (since the field vanishes). We now examine the clash of symmetries in the case 
of the SU(5) x Z2 q = 2 kink. 

The general form of p? is given in Eq. (2.13) with the profile functions in 
Eq. (2.23). Then 


pP = 0) = MÊ « diag(0, 1,1, —1, —1) (2.32) 


The symmetries within the kink are given by the elements of SU(5) x Z2 that 
leave M? invariant. Hence the internal symmetry group consists of two SU (2) 
factors, one for each block proportional to the 2 x 2 identity, and two U (1) factors 
since all diagonal elements of SU (5) commute with MÊ. Therefore the symmetry 
group inside the SU (5) x Zz kink is [SU(2)]? x [U(1)]?. This is smaller than the 
SU(3) x SU(2) x U(1) symmetry group outside the kink. 

The conclusion that the symmetry inside a kink is smaller than that outside holds 
quite generally [164]. Classically this would imply that there are more massless par- 
ticles outside the kink than inside it. However, when quantum effects are taken into 
account this classical picture can change because the fundamental states in the out- 
side region may consist of confined groups of particles (““mesons” and “hadrons”) 
that are very massive [51]. If a particle carries non-Abelian charge of a symme- 
try that is unbroken outside the wall but broken inside to an Abelian subgroup, it 
may cost less energy for the particle to live on the wall. This is because it may be 


3 As in Section 2.4 we could have found the symmetry group inside the kink by finding those transformations in 
G321— that are also contained in G321+. 
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unconfined inside the wall where it only carries Abelian charge, while it can only 
exist as a heavy meson or a hadron outside the wall.* 


2.7 Interactions of static kinks in non-Abelian models 


The interaction potential between kinks found in Section 1.8 is easily generalized 
to kinks in non-Abelian field theories. Following the procedure discussed in that 
section, the force in the SU(5) x Z3 case is 


Fa =[-m)- We) +V@]" (2.33) 


where —a — R and —a + R are defined in Fig. 1.3. Evaluation of F yields an 
exponentially small interaction force whose sign depends on Tr(Q, Q2) [121] where 
Q; and Q, are the topological charges of the kinks. If the Higgs field at x = —oo 
is B_, between the two kinks is ®ọ, and is By at x = +00, then Q; x By — B_ 
and Q2 x 4 — Bg (See Eq. (1.8)). 

What is most interesting about the interaction is that a kink and an antikink 
can repel. Here one needs to be careful about the meaning of an “antikink.” An 
antikink should have a topological charge that is opposite to that of a kink. That 
is, a kink and its antikink together should be in the trivial topological sector. But 
this condition still leaves open several different kinds of antikinks for a given 
kink. To be specific consider a kink-antikink pair, where the Higgs field across the 
kink changes from ®(—oo) x +(2, 2,2, —3, —3) to ®(O) x —(2, —3, —3, 2, 2). 
(Here we suppress the normalization factor and the “diag” for convenience of 
writing.) There can be two types of antikinks to the right of this kink. In the first 
type (called Type I) the Higgs field can go from ®(0) x —(2, —3, —3, 2, 2) to 
(+20) x +(2, 2, 2, —3, —3), which is the same as the value of the Higgs field at 
x = —oo and thus reverts the change in the Higgs across the kink. In the second 
type (Type II), the Higgs field can go from ®(0) x —(2, —3, —3, 2, 2) to B(+00) « 
+(—3, 2, 2, —3, 2). Now the Higgs at x = +00 is not the same as the Higgs at 
x = —ox, but the two asymptotic field values are in the same topological sector. 

By evaluating Tr(Q; Q2), where Qı and Q% are the charge matrices of the two 
kinks, it is easy to check that the force between a kink and its Type I antikink is 
attractive, but the force between a kink and its Type II antikink is repulsive. The 
q = 2 kinks can have charge matrices Q“ that we list up to a proportionality factor 


Q® = (—4, 1,1,1,1), Q® = (1, —4, 1,1, 1), Q® = (1,1, —4, 1, 1), 
Q = (1,1, 1, —4, 1), Q® = (1,1, 1,1, -4) (2.34) 


4 Localization of particles to the interior of defects has led to the construction of cosmological scenarios where 
our observed universe is a three-dimensional defect or “brane” embedded in a higher dimensional space-time. 
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Stable antikinks have the same charges but with a minus sign. Then, one can take 
a kink with one of the five charges listed above and it repels an antikink that has 
the —4 occurring in a different entry because Tr(Q;Q2) > 0. Hence, there are 
combinations of kinks and antikinks for which the interaction is repulsive. Further, 
in a Statistical system a kink is most likely to have a Type II antikink as a neighbor 
and such a kink-antikink pair cannot annihilate since the force is repulsive. 

The result that the force between two kinks is proportional to the trace of the 
product of the charges extends to other solitons (e.g. magnetic monopoles) as well. 
In this way, the forces between certain monopoles with equivalent magnetic charge 
can be attractive whereas normally we would think that like magnetic charges repel, 
and between certain monopoles and antimonopoles can be repulsive. 


2.8 Kink lattices 


In this section we describe the possibility of forming stable lattices of domain 
walls in one spatial dimension and the consequences in higher dimensions. Our 
discussion is in the context of the Ss; x Z} model though similar structures have 
been seen in other field theory models as well [92, 43]. 

We know that Z2 topology forces a kink to be followed by an antikink. Then 
we can set up a sequence of kinks and antikinks whose charges are arranged in the 
following way 


-..0%99Q%Q%Q9IG® _.. (2.35) 


where Q“ and Q“ refer to a kink and an antikink of type i respectively (see 
Eq. (2.34)). Alternately, this sequence of kinks would be achieved with the following 
sequence of Higgs field vacuum expectation values (illustrated in Fig. 2.5) 


2.2 —(2, 2,2, —3, —3) — +(2, —3, —3, 2, 2) 

—> —(-3, 2,2, —3, 2) 

> +(2, —3, 2, 2, —3) 

—> —(2, 2, —3, —3, 2) 

> 4+(-3, —3, 2,2, 2) 
== (2,2, 2; —3, —3) > ... (2.36) 
The forces between kinks fall off exponentially fast and hence the dominant forces 
are between nearest neighbors. As discussed in the previous section, the sign of the 
force between the ith soliton (kink or antikink) and the (i + 1)th soliton (antikink 
or kink) is proportional to Tr(Q;Q;i+1) where Q; is the charge of the ith object. 


For the sequence above, Tr(Q; Qi+1) > 0 for every i and neighboring solitons repel 
each other. In particular, they cannot overlap and annihilate. 
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Figure 2.5 Inthe lattice of kinks of Eq. (2.36), the vacua are arranged sequentially 
in a pattern so as to return to the starting vacuum only after several transitions 
between the two discrete (Z2) sectors. 


The sequence of kinks in Eq. (2.35) has a period of six kinks. These six kinks 
have a net topological charge that vanishes since the last vacuum expectation value 
in Eq. (2.36) is the same as the first value. Hence we can put the sequence in a 
periodic box, i.e. compact space. This gives us a finite lattice of kinks. 

The sequence described above has the minimum possible period (namely, six). 
It is easy to construct other sequences with greater periodicity. For example 


00O 0O GOVO DD QI HO... (2.37) 


is a repeating sequence of ten kinks. 

The lattice of kinks is a solution in both the S5 x Z3 and the SU (5) x Z models. 
However, it is stable in the former and unstable in the latter. The instability in the 
SU(5) x Zz model occurs because a kink of a given charge, say Q®?, can change 
with no energy cost into a kink of some other charge, for example Q". Then, in 
the sequence of Eq. (2.35), the third kink changes into Q"”, then annihilates with 
the antikink with charge Q™ on its right. In this way the lattice can relax into the 
vacuum. In the S$; x Z2 case, however, the degree of freedom that can change the 
charge of a kink is absent and the lattice is stable. 

So far we have been discussing a kink lattice in one periodic dimension. This is 
equivalent to having a kink lattice in a circular space. Next consider what happens 
in a plane in two spatial dimensions. A circle in this plane can once again have a 
kink lattice since neighboring kinks and antikinks repel. However, when extended 
to the whole plane, the kink lattice must have a nodal point as shown in Fig. 2.6. In 
three spatial dimensions, the nodal points must extend into nodal curves.> 

We shall discuss kink lattices further in Chapter 6. 


5 This is very similar to the case where several domain walls terminate on topological strings, except that there 
are no topological strings in the model. 
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Figure 2.6 A domain wall lattice consisting of six domain walls can be formed 
in a one-dimensional sub-space (dashed circle) of a two-dimensional plane. This 
domain wall lattice is stable. Extending it to the two-dimensional plane, the differ- 
ent domain walls converge to a nodal point. This implies that the S5 x Z2 model 
contains domain wall nodes (or junctions) in two dimensions and nodal curves in 
three spatial dimensions. 


2.9 Open questions 


1. Discuss all topological and non-topological kink solutions in an SU(N) x Z2 model 
where N is even. In [163] the case with odd N is discussed.® 


© However, it is incorrectly stated that the Z2 symmetry is included in SU(N) when N is even, as can be seen 
from the Tr(®+) argument of Section 2.1. 
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Interactions 


In the previous two chapters, we have described kink solutions in several models 
but these solutions have mostly been discussed in isolation. In any real system, there 
is a variety of kinks and antikinks, in addition to small excitations (particles) of the 
fields. The interactions of kinks with other kinks and with particles play an important 
role in the evolution of the system. The motion of kinks is also accompanied by 
the radiation of particles. Ambient particles in the system scatter off kinks, and 
kinks collide with each other, and perhaps annihilate into particles. As discussed 
in Section 1.9, in some models a kink-antikink pair can bind together to form a 
non-dissipative solution which is called a “breather.” In other models, approximate 
breather solutions have been found, which play an important role in the scattering 
of a kink and an antikink. These topics are discussed in the following sections. 


3.1 Breathers and oscillons 


So far we have been considering kinks, which are static solutions to the equations 
of motion. In the sine-Gordon model of Eq. (1.51), a one-parameter family of non- 
static, dissipationless solutions is also known. These are bound states of a kink and 
an antikink and are called breathers. The breather solution was described briefly in 
Section 1.9 and can be re written as 


Polt, x; v) = 4 an”! anu Gti) l (3.1) 


B vcosh(./ax//1 + v2) 


where v is a free parameter (see Fig. 3.1). We will have more to say about breathers 
when we quantize kinks in Chapter 4 as they play a very fundamental role in the 
novel duality between the sine-Gordon model and the massive Thirring model (see 
Section 4.7). 

Breather solutions are not known to exist in the 4¢* model [135]. However, 
numerical studies of the scattering of a Z kink and antikink revealed the existence of 
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Figure 3.1 The sine-Gordon breather shown at various times during one oscil- 
lation period. At certain times, the field profile is that of a separated kink and an 
antikink. At other times, the kink and the antikink overlap and cannot be distin- 
guished. 


extremely long-lived, oscillating bound states of kinks and antikinks [4, 19, 26, 64]. 
The existence of kink-antikink bound states has been interpreted as a resonance 
phenomenon between the natural excitation frequency of the kink profile (shape 
mode) and the frequency of oscillation of the bound kink-antikink system. Radiation 
from a time-dependent scalar field configuration will be suppressed if the oscillation 
frequency of the configuration is small compared to the mass of the radiation quanta 
and this can be used to understand the longevity of oscillons (Farhi, 2005, private 
communication). 

The simplest hypothesis is that oscillons are approximate breather solutions 
since a region of the sine-Gordon potential and the A¢* potential have very similar 
shapes. We can compare the two potentials when the sine-Gordon potential has been 
shifted so that it has a maximum at @ = 0. The parameter 6 in the sine-Gordon 
model is chosen so that the first positive minimum is at ¢ = +7. @ is fixed by 
requiring that the masses of small excitations in the true vacua, given by the second 
derivative of the potential, are equal in the two models. Then the two potentials are 
given by 


A 
Vz(b) = 4-9 (3.2) 


Vic(®) = zz (1 — cos(B(d — n))) (3.3) 


a 
p? 
with 


a = 2dr’, p= (3.4) 


3/a 
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Figure 3.2 The A¢* potential (broken curve) and the shifted sine-Gordon poten- 
tial (solid curve) when the parameters are chosen so that the vacua occur at the 
same values of ġ and the curvatures of the potentials at the vacua are also equal. 


Figure 3.3 The profiles of the kinks in the 4¢* model (broken curve) and the 
shifted sine-Gordon model (solid curve) with potentials as shown in Fig. 3.2. 


The two potentials can be compared in the vicinity of their true vacuum at @ = 7. 
Then 


À 
Vz,(b) = An*(p — m? + anlo = n + 7G — ny" (3.5) 
and 
2n? 
4! 
In Fig. 3.2 we show these two potentials and in Fig. 3.3 we compare the kink 
profiles. 


We will return to the breather and its role in the quantum sine-Gordon model at 
the end of Section 4.7. 


À 
Valo) = An?(b — ny? — ($ — nf + Oo — n)°) (3.6) 
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3.2 Kinks and radiation 


By “radiation” we mean propagating excitations of small amplitude of a field, 
which in this chapter will be taken to be the same field that makes up the kink. 
Asymptotically, these excitations have the usual plane wave form: exp(i(wt + kx)). 
In the kink background, these “scattering states” are found as solutions to the 
equation of motion for fluctuations about the kink. If we denote the kink solution 
by x(x), the fluctuation field w(t, x) is 


W(t, x) = b(t, x) — px) (3.7) 


We will assume || < (p), where (@) is the vacuum expectation value of ø. To 
find the scattering states, we take y = f(x)e~'@’ where it is understood that the 
real or imaginary part should be taken — in other words, the physical modes are 
[fixe + f*(x)et@"]. Perturbing the Lagrangian for ¢ (first line of Eq. (1.2)), 
we find that f(x) satisfies the (linearized) equation of motion 


Hf =-f"+U(a)f = f (3.8) 


where 
2 


U(x) = Veo) = 
IP" | 6d 
The scattering states around a static kink are obtained by solving the Schrédinger- 
type equation, Eq. (3.8), which for some potentials, falls in the general class of 
equations discussed in Appendix C. 
We now consider the Z2 kink for which the potential U is obtained from Eqs. (3.9) 
and (1.2) to be 


(3.9) 


U(x) = (36; — n°) (3.10) 


We now list the eigenmodes of Eq. (3.8). (We will encounter them again in Chap- 
ter 4.) First, there are two bound states, also known as “discrete” modes: 


wo = 0, fo= sech?z (3.11) 
V3 


o1 = My, fi = sinhz sech?z (3.12) 


where z = x/w = myx /2. The w = 0 mode is called the “translation mode” and the 
second is the “shape mode.” Then there is a continuum of states for my < œw < oo 
which are the scattering states: 


fk = [3 tanh?z — 1 — w?k? — i 3wk tanhz] (3.13) 
The dispersion relation is 


or =k +m, (3.14) 
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We are now interested in processes that involve both a kink and the scattering 
states (radiation). For example, if a kink accelerates, it will emit radiation. What 
is the radiated power? The answer will depend on the forces that make the kink 
accelerate and whether or not these forces deform the structure of the kink.! We 
shall examine the radiation from kink shape deformations and other interactions of 
kinks and radiation after a brief diversion in the next section. 


3.3 Structure of the fluctuation Hamiltonian 


In this section we will show two interesting properties of the fluctuation 
Hamiltonian, H, defined in Eq. (3.8). The first is that the potential U(x) has a 
very special form that implies that the Hamiltonian can be factored. The second is 
that there exists a “partner Hamiltonian” with (almost) the same spectrum as the 
original Hamiltonian. 

The special form of U follows from the fact that the kink has a translation zero 
mode (see Section 1.1). Hence there exists an eigenstate with @ = 0. Denote this 
“translation mode” by yr. Hence 


Hyr = (38 + U(x) = 0 (3.15) 

Therefore 
U(x)= ve 3.16 
(x) a (3.16) 


which can also be rewritten as 


U= f+, f= y (3.17) 


For the particular cases of the Z2 and sine-Gordon kinks, not only is U(x) 
of the form in Eq. (3.17) but it is also reflectionless. Then, an incident wave is 
fully transmitted and the reflection coefficient vanishes. In this case, the only non- 
trivial characteristic of scattering states is that the waves get a phase shift owing to 
the presence of U(x). This property will be useful when we quantize the kink in 
Section 4.1. 

The Hamiltonian H with a potential of the form in Eq. (3.17) has the important 
property that it can be factored 


H = AtA = (+8 + fX—3 + f) (3.18) 
Therefore the equation for the eigenstates is simply 


Hf =AtAf =% f (3.19) 


1 In the case of domain walls in three spatial dimensions, the curvature of the wall is itself responsible for 
acceleration. This motion leads to the emission of scalar and gravitational radiation and will be discussed in 
Chapter 8. 
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The factorization has the consequence that one can readily construct a “partner” 
Hamiltonian, H_, that has almost an identical eigenspectrum as H. This partner 
Hamiltonian is H- = AAT. If f; is an eigenstate of H with eigenvalue w, then 
Af; is an eigenstate of H_ with the same eigenvalue. This argument works for all 
eigenstates except the one for which Af; = 0. Hence H has a single extra eigenstate 
with œ = 07 

The potential U(x) determines the spectrum of excitations around a soliton. 
The factorizability of the Hamiltonian is useful in the problem of reconstructing 
V(@) from the spectrum of fluctuations (i.e. the set of o’) using inverse scattering 
methods [165]. 


3.4 Interaction of kinks and radiation 


As remarked below Eq. (3.17) the potentials U(x) for both the Z2 and the sine- 
Gordon kinks are rather special since they are reflectionless. All that happens is 
that the transmitted wave gets phase shifted. This is equivalent to a time delay in 
the propagation of the wave through the kink. 

From the solution for the scattering states given in Eq. (3.13) for the Z2 kink we 
find a momentum dependent phase shift 


3wk 
= -1 
ôklz, = 2 tan (ss) (3.22) 
This corresponds to a time delay 
Ox 2 —1 3wk 
= = t Ss 3.23 
ala =” mn (5) 83.23) 


Z [k2 + mi, 


Similarly the phase shift and time delay in the case of the sine-Gordon kink are 


k 
ôklsg = x — 2 tan! (5) (3.24) 


= x — tan — (3.25) 
sG [k2 + m?, l My 


2 The two partner Hamiltonians can also be combined to form a supersymmetric Hamiltonian, Hss 
AtA 0 
Hs=(“9° gar) =(2,0%} =00*+0t0 (3.20) 


where 


o=(4 an pre a (3.21) 
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While there is no reflection of radiation of the same field that makes up the kink 
in the Z% and sine-Gordon cases, there can be reflection of fluctuations of other 
fields [53]. As an example [171], consider a second scalar field x included in the 
Z model so that the full Lagrangian becomes 


1 m? oO 
L = La + (0,7) — x? — -°x 3.26 
$ + 5u > X z? X (3.26) 


where Lọ is the Lagrangian for the Z2 model (Eq. (1.2)). Then the scattering modes 
of x in the presence of a Z2 kink are found by solving 


ax — aX +m? x +ogKx =0 (3.27) 
Substituting k = ntanh(x/w) and x = exp(—iwt) f(x), we get 
ay f + (v? — sech (X) f =0 (3.28) 


2 2 


where X =x/w, v = wo? — m4 — on’), 6 =0n?w?. (Recall that w = 


f2/An?.) 

Equation (3.28) is a special case of the differential equation described in Ap- 
pendix C. The scattering state is found for real values of v and has the asymptotics: 
f — e for x — œ, and for x > —oo: 


z r — ikw)P(—ik)ei* rd — ikw) (ikje 
Td/2+y -ikr (1/2-y-—ik) Td/2+y)Pd/2-y) 
where k = v/w and y = /o + 1/4. 


The reflection coefficient can be read off from the asymptotic behavior of f(x) 
as x — —oo and has been evaluated in Section 12.3 of [113] 


f 


(3.29) 


1 +cos(2xy) 


= 3.30 
cosh(2x k) + cos(27y) ( ) 
The transmission coefficient is 
2 sinh? (rk 
= eee) =1-—R (3.31) 
cosh(2x k) + cos(27y) 


From the asymptotic expression in Eq. (3.29), it is also possible to calculate the 
time delay of the reflected and transmitted waves owing to the kink. For example, 
if we write 

Td — ikw) (—ik) 
r(1/2+ y — ik) (1/2 — y — ik) 
where T is the transmission coefficient above, then the time delay of the transmitted 
wave is given by d,/a. 


= JT |e} (3.32) 
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3.5 Radiation from kink deformations 


A Static kink does not emit any radiation. Nor does it emit radiation if it is moving 
at constant velocity (see Eq. (1.10)). However, if the kink is accelerating (owing to 
some external force), or its shape is deformed, it can emit radiation in the form of 
scalar particles [106, 107]. In3 + 1 dimensions, acceleration and deformations arise 
since the kinks (domain walls) are moving under their own tension, except in the 
very special cases of static solutions. The radiation emitted from curved domain 
walls has not been calculated analytically, though the problem has been studied 
numerically [182]. In the case of 1 + 1 dimensional Z, kinks that are undergoing 
periodic deformations, the radiation has been found analytically in [110, 140], and 
we shall describe this calculation below. 

Following [110], we simplify notation by setting A = 2and7 = 1 inthe Zz model 
so that w = 1 in these units (see Section 1.1). Then the field @(x, t) is written in 
terms of the complete set of small excitations. This gives 


P(X, t) = P(X) + RO) fo) + AOAi@) + FO, t) (3.33) 


where x = tanh(x), fo and fı are the translation and shape modes respectively as 
given in Eqs. (3.11) and (3.12), R(t) and A(t) are their time-dependent amplitudes, 
and the function f(x, t) contains all the continuum states around the kink. The 
frequency of oscillation of R(t) is wọ = 0 and of A(t) is œw = /3. These values 
were derived for linearized fluctuations about the kink. Non-linearities will modify 
w = J3 but we assume that such modifications are small. 

We will work in the rest frame of the kink and so 


R(t) =0 (3.34) 


The idea now is to insert Eq. (3.33) in the equation of motion for @ with some 
choice of the amplitude A(t) which is assumed to be small, and then to find the 
solution for the scattering states, f(x, t), which form the radiation. 

Insertion of Eq. (3.33) in Eq. (1.4) gives 


(A+3A)fi + f — f” +2(392 — I) f = -6f + po 7A? 
—6(f + 2b) ffid — 2f? A? 
— 6, f? — 2f? (3.35) 


where the equations satisfied by ¢, and fı have been used. Assuming that A is small, 
and that f is O(A”) or smaller, the leading order equation is A+ 3A = 0+ O(A’). 
Then to order A’, the equation for f is 


(A+3A)fi+ f — f" +230 - I) f = —6p ff A’ (3.36) 
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The f-independent terms on the right-hand side of Eq. (3.35) are source terms 
which cause radiation. Hence f will not be zero at order A”. The terms will also 
cause the amplitude, A, of the shape mode to depart from the purely oscillatory 
behavior. To determine how much of the source affects radiation and how much 
affects the shape mode, note that fı and f are orthogonal 


J dx fi(x) f(x, t) = 0 (3.37) 


So we can decompose the equation into a direction parallel to fı in mode space and 
orthogonal to it. One assumption we have to make is that the back-reaction of the 
radiative modes on the shape mode is higher order in A. For example, Eq. (3.37) 
does not by itself imply that f” and fı are orthogonal. Then, multiplying Eq. (3.36) 
by fı and integrating over all space gives 


A+3A = -6A? J dxok f? = —6a A? (3.38) 
provided we have normalized fı so that 
[extreor =1 3.39) 
Explicit evaluation gives 
3a 
— 3.40 
oY ay Ceo 


The equation orthogonal to fı is 
f-—f" +2(3¢¢ — I) f = 66x f? A? + 6af A? (3.41) 


and this will determine the radiation from the deformed kink once we have specified 
A. 
The leading order solution for A is 


A = Ao cos( 73t) (3.42) 


Hence 
Az 
A= 5 leos(2v31) $i] (3.43) 


This form implies that the source for f in Eq. (3.41) has a time-dependent piece 
and another time-independent piece. Since the equation is linear in f, only the 
time-dependent piece proportional to cos(2./3r) is important. Setting 


f(x, t) = Ref F(x)) (3.44) 
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the equation that needs to be solved is 
3 
=F" + (6G, — 2-0) F = S (afi — bef) Ager (3.45) 


where wy = 24/3. Since the left-hand side is time-independent, this only has solu- 
tions for 


w = @ = 2V3 (3.46) 


and then all the solutions of the homogeneous equation are known (see Appendix C; 
[113, 126]). The solutions of the homogeneous equation with plane wave asymp- 
totics are 


F,(x) = Bp — 1 —q* — 3iqdx)e'* (3.47) 


where q = Vw” — 4. Knowing all the solutions of the homogeneous equation, it is 
possible to explicitly construct the (retarded) Green’s function suitable for outgoing 
radiation. 


Basie | —F_q(y)Fy(x)/W. (x <y) sae 
—F,(y)F_q)/W, (x > y) 
where W is the Wronskian 
Ws Fœ) F_q (x) — Fy (x) Fax) (3.49) 


The Wronskian is a constant and its value can be found by using the explicit solutions 
W = —2ig(q* + Iq? +4) (3.50) 


The solution of the inhomogeneous equation (3.45) is found by convoluting the 
source with the Green’s function 


+00 3 
F(x) = dy Gtx, y)5[afily) = PONFTO)]AG (3.51) 


oo 
With a little more manipulation, we obtain the radiation field in the x — +00 limit 
—3 A2ei(or—4x) +00 
2ig(2 — q? — 3ig) J- 
with œ = wy = 2V3 and q = Vw? — 4 = 2V2. The integral can be done explicitly 


leading to 
tepa TA E T T E, (3.53) 
7 32 sinhGrg/2)V q2 +1 ° 2 i 


f(x, t) = Re | dy KOOFO)] (3.52) 
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The phase 6 can be read off from Eq. (3.52) because the integral is purely imaginary 
and does not contribute 


3q 

5 = tan '( ——— 3.54 

on (z F z) E 

Now that we have the solution for the radiation field, we can find the energy 

flux by using the To; components of the energy-momentum tensor in Eq. (1.39). 

Including a factor of 2 to account for the radiation toward x — —o0o, we obtain the 
radiated power [110] 

dE 


x” —0.0204ź4 (3.55) 


The back-reaction of the radiation on the deformation amplitude can be estimated 
on the grounds of energy conservation. In [110] the results above are compared to 
the results of a numerical evolution of the deformations using the full non-linear 
equations with good agreement. 


3.6 Kinks from radiation 


By time reversing kink and antikink annihilation, it should be possible to obtain 
kink-antikink creation from incoming radiation. However, the stream of incoming 
radiation would have to be sent with just the correct phase relationship and energy. 
Such initial conditions occupy zero volume in the space of all initial conditions. A 
more physical problem is to identify the set, or a large subset, of initial conditions 
for the incoming radiation that will lead to kink-antikink creation. This problem is 
unsolved. Yet certain interesting results have been obtained in [110] in the “gradient 
flow” approximation in which the second time derivative terms in the equation of 
motion are neglected. 

Consider the collision of two kinks in the presence of a pre-existing kink [110], 
as depicted in Fig. 3.4. As the kink-antikink-kink (kkk) system evolves, a kink- 
antikink annihilate, and we are left with a kink whose shape is excited. Reversing 
this process, if we start with a kink whose shape is excited, it can produce a kink- 
antikink pair. In [110], this relation between the shape mode and the creation of a 
kink-antikink pair was explored. 


3.7 Scattering of kinks 


The sine-Gordon model is a famous example of a completely integrable system [48]. 
Sine-Gordon kinks are examples of “solitons” in the strict mathematical sense in 
which when two or more solitons (or anti-solitons) scatter, they simply pass through 
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Figure 3.4 A kink collides from the left with another kink coming in from the 
right in the presence of an antikink in the middle. The time evolution of the field is 
shown in succession by the solid, dotted, dashed, and dashed-dotted curves. The 
evolution shows that a kink and antikink annihilate leaving behind a kink whose 
shape mode is excited (dotted and dashed curves). With further evolution, the shape 
mode will dissipate and an unexcited kink will remain as seen in the dashed-dotted 
curve. [Figure reprinted from [110].] 


each other. The only consequence of the scattering is that there is a phase shift, or 
equivalently, a time delay. The time delay may be understood by realizing that the 
force between two kinks in the sine-Gordon model is attractive. Hence the kinks 
collide and form a bound state for some time. The time delay may be viewed as the 
time spent by the kinks in the form of a bound state. A crucial aspect of the scattering 
is that there is no dissipation. More details can be found, for example, in [48]. 

Kink scattering in the Z2 model has a more complex character. In this case, we 
cannot have kink-kink scattering because two Z2 kinks cannot be adjacent to each 
other. Instead, we need only consider kink-antikink scattering. This has been the 
subject of significant investigation [26, 4]. When a kink-antikink collide, the only 
possibilities are that they reflect back or they annihilate (see Fig. 3.5). 

As we might expect, at very low incoming velocities, a bound state is formed and 
annihilation inevitably occurs, while at very high velocities, reflection takes place. 
The remarkable discovery of numerical studies of kink-antikink scattering is that 
the change from annihilation to reflection does not happen at just one critical value 
of the incoming velocity. Instead there are bands of incoming velocity at which 
annihilation occurs, while at other values of the incoming velocity the kink and 
antikink are reflected. The plot in Fig. 3.6 shows these results. 

The unexpected dependence of kink-antikink scattering on the incoming velocity 
has been examined closely in [26, 4]. The behavior is understood as a resonance 
effect between oscillations of the mode that describes the shape distortions of the 
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Figure 3.5 A kink and an antikink with incoming velocity v; are shown in panel 
(a). The two possible outcomes of the scattering are shown in panels (b) and (c). 
In panel (b), the kinks scatter and reflect. Their outgoing velocity vy need not be 
equal to v;. In panel (c), the kink and antikink have annihilated and radiated away 
their energy, leaving behind the trivial vacuum. In both outcomes, the scattering 
is likely to be accompanied with radiation that has not been depicted. 
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Figure 3.6 The ratio of outgoing to incoming kink velocities after scattering 
versus the incoming velocity [26, 4]. When the outgoing velocity is plotted to be 
zero, a kink-antikink bound state is formed that decays to the vacuum by radiation. 
Notice that the kink-antikinks annihilate in certain bands in the initial velocity. 
[Figure reprinted from [26].] 
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Figure 3.7 Two curved domain walls collide and intercommute. At the collision 
point, there is lots of radiation owing to annihilation or owing to the formation of 
a closed domain wall that then collapses and decays into radiation. To imagine the 
walls in three dimensions, rotate the figures along the horizontal axis. In the initial 
state the two curved walls are disconnected from each other while in the final state, 
the wall is in the shape of a “wormhole,” with a sphere in the middle. 


kinks and the oscillations of the kinks as a whole owing to kink-antikink interactions. 
We shall not describe the details of the analysis here. 

The scattering of SU(5) x Zz kinks has been studied numerically in [121]. In 
this case, there is an additional degree of freedom, namely the non-Abelian charge 
of the colliding kinks (or “color”) and there are a variety of initial conditions that can 
be considered. For the stable variety of kinks — the q = 2 kinks (see Section 3.2) — 
the scattering of kinks and antikinks of the same color is qualitatively similar to that 
of Z2 kinks. If the colors are different, however, there is a repulsive force between 
the kinks and they are observed to bounce back elastically. 


3.8 Intercommuting of domain walls 


We finally consider the collision of two domain walls. The outcome is found by 
numerical evolution of the equations of motion. As the walls come together, they 
reconnect along the curve of intersection [136] as shown in Fig. 3.7. This process 
is called “intercommuting” or, simply, “reconnection.” 


3.9 Open questions 


ja 


. Suppose we want to create a well-separated Z2 kink-antikink pair by colliding small 
amplitude plane waves (particles) in the ¢ = +n vacuum. What conditions must be 
imposed on the incoming waves? What is the space of initial conditions that leads to 
domain wall formation? Can the initial conditions be implemented in a practical setting 
(e.g. accelerator experiments)? 

2. Can a domain wall lattice be generalized to other defects, e.g. a lattice of strings and 

monopoles? 
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. Study the interaction of domain walls and strings/magnetic monopoles in a model that 
contains both types of defects e.g. the SU(5) x Z2 model has walls and magnetic 
monopoles. 

. Construct a sine-Gordon-like breather field configuration in the 4+ model. This will not 
be an exact solution of the field equation of motion. Hence it will radiate. Calculate the 
radiated power. Are there circumstances in which the radiated power is very small? 

. Can the analysis of radiation from kink deformations be extended to the case of oscillating 
domain walls? The simplest procedure would be to decompose the field as in Eq. (3.33) 
and to include a suitable external (harmonic) potential that drives the translation mode 
only. This will cause the kink to oscillate as a whole without deformations. However, the 
oscillations will source the shape mode and the radiation, and an analysis of the kind in 
Section 3.5 seems feasible. 

. Can the analysis of radiation from kink deformations be extended to the case of spherical 
domain walls? 

. How can the radiation analysis be extended to vortex solutions in two or more spatial 
dimensions? 


4 
Kinks in quantum field theory 


A particle in a classical harmonic oscillator potential, mw*x*/2, has minimum en- 
ergy when it sits at rest at the bottom of the potential. Then the particle’s energy 
vanishes. The Heisenberg uncertainty principle however modifies this picture for 
the quantum harmonic oscillator. The particle cannot sit at rest (with definite mo- 
mentum) at the bottom of the potential (a definite location). Indeed, the quantum 
zero point motion lifts the ground state energy to w/2. Further, the excited states 
of the simple harmonic oscillator are discrete and occur at energies (n + 1/2)a, 
n=0,1,2,... 

Just as the classical harmonic oscillator is modified by quantum effects, any 
classical solution to a field theory is also modified by quantum effects. Quantum 
effects give corrections to the classical kink energy owing to zero point quantum 
field fluctuations. These quantum corrections are small provided the coupling con- 
stant in the model is weak. To “quantize the kink” means to evaluate all the energy 
levels of the kink (first quantization) and to develop a framework for doing quantum 
field theory in a kink background. This involves identifying all excitations in the 
presence of the kink and their interactions. The field theory of the excitations in 
the non-trivial background of the kink is akin to second quantization. Finally, one 
would also like to describe the creation and annihilation of kinks themselves by 
suitable kink creation and annihilation operators. This would be the elusive third 
quantization. 

Initially we calculate the leading order quantum corrections to the energy of the 
Z and sine-Gordon kinks. As these two examples illustrate, the precise value of 
the quantum correction depends on the exact model and kink under consideration. 
Yet there is one common feature — quantum corrections tend to reduce the energy 
of the kink. This result is quite general and we prove it using a variational argument 
in Section 4.5. 

The quantum corrections to the kink mass are obtained by using a perturbative 
analysis where the coupling constant is the expansion parameter, as first done 
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in [38, 42]. For fixed values of the masses of particles in the field theory, the 
energy of the classical solution is proportional to the coupling constant raised to 
a negative power (for example, see Eq. (1.20)) and so the perturbative analysis 
holds only if the kink is much more massive than the particles in the model. As the 
coupling constant is increased, quantum effects become stronger and eventually 
the perturbative scheme breaks down. Remarkably, the sine-Gordon model is still 
amenable to analysis in this regime and, at strong coupling, the sine-Gordon kinks 
become lighter than the particles. Indeed, there exists a weakly coupled description 
of the model in which perturbative methods can be used: this is the massive Thirring 
modelin which the particles (low energy excitations of a fermionic field) correspond 
to the sine-Gordon kinks at strong coupling (see Section 4.7). 

The phenomenon observed in the sine-Gordon and massive Thirring models, in 
which solitons of one model (Model 1) are identified with the particles of a second 
model (Model 2) and vice versa in certain regimes of the coupling constants, is 
known as “duality.” Model 1 is said to be dual to Model 2 if the particle-plus-soliton 
spectrum of Model 1 maps onto the soliton-plus-particle spectrum of Model 2 and 
vice versa. Both models describe the same physics, except that the light and heavy 
degrees of freedom are interchanged. 

The Z2 model does not share the remarkable symmetries of the sine-Gordon 
model and less is known about the Z3 kink at strong coupling. However, the mass 
of the Z, kink can be evaluated at strong coupling using lattice field theory. We 
describe these results in Section 4.8 and conclude, once again, that the kink becomes 
less massive as the coupling is increased and eventually becomes massless. 

In this book, we only describe quantization of the mass of the kink using canonical 
techniques. A more extensive discussion of various other techniques and issues can 
be found in [35, 126] and in the series of papers in [38, 42]. 


4.1 Quantization of kinks: broad outline 


In this section, we evaluate the contribution of the zero point fluctuations to the 
energy of the kink. Then we briefly discuss excited states. 
The quantization procedure can be outlined as follows: 


e Consider a field theory in two dimensions with compact spatial dimension of size 
L, assumed large compared to any other length scale in the problem. Periodic 
boundary conditions are imposed on the fields. Eventually take L — ov. 

e Consider small quantum fluctuations, y, about the classical kink background, 


Gx, 
P(t, x) = d(x) + Wt, x) (4.1) 
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n=-1l 
n=+41 
n=0 


Figure 4.1 A trivial potential on a periodic space with period L is shown on the 
left. The field modes are labeled by an integer n = 0, +1, +2, ... When there is a 
kink, the potential felt by the modes becomes non-trivial as depicted by the curved 
bottom of the figure on the right. What used to be the n = 0 mode in the trivial 
potential (on the left) becomes the lowest bound state, the zero mode, in the non- 
trivial potential. Similarly a linear combination of the n = +1 modes in the trivial 
box may become a second bound state (n = +1 in this illustration) and the other 
states remain unbound but shift in form and energy. 


Linearize the equation for yw and then quantize, that is, regard the field y as a 
quantum operator 


W(t, x) = lax felt. x) + af ARC, x) (4.2) 


where al and a, are creation and annihilation operators. The fk are mode functions 
i.e. orthonormal solutions of the linearized equations of motion for y in the kink 
background. The equation satisfied by f is 


0? fk — 2 fk + VAD fc = 0 (4.3) 


Find all the eigenmodes, fg, of the fluctuations and their eigenfrequencies wp. As 
shown in Fig. 4.1, in the presence of the kink the modes are displaced. Some of 
the low-lying modes without the kink become bound states in the presence of the 
kink, and the others become scattering states as L —> oo. 

Each eigenmode corresponds to a quantum harmonic oscillator with zero point 
fluctuations. Sum up the zero point energies of all the modes to get the quantum 
correction to the classical kink energy, Eq, 


% 1 
E = Ea + D i (4.4) 


In the L — œ limit, the sum over the modes becomes a sum over bound states 
and an integral over scattering states. Also note that Eq. (4.4) is only valid to 
leading order in the quantum corrections since we have ignored interactions of the 
fluctuation field, y. 
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In following this procedure, quantum field theoretic subtleties arise. 


e The zero point energy of the trivial vacuum (without the kink) must be subtracted 
from the zero point energy of the kink since we want to define the energy of the 
trivial vacuum to be zero. Therefore 


1 1 
E= Eat )) 5% - [Eso +>) Ta 


where Eeo is the classical energy of the trivial vacuum and is chosen to vanish 
(Eao = 0), and o are the eigenfrequencies of the modes in the trivial vacuum. 
¢ The energy must be expressed in terms of renormalized parameters. 


In the trivial vacuum, the energy eigenvalue for the mode with n nodes is 


oO = Jk? + mi, (4.5) 


where k, = 22n/L and n € Z, the set of all integers. Now suppose that the kink 
potential V” is turned on slowly, i.e. that the potential term in Eq. (4.3) is multiplied 
by a parameter that vanishes for the free field theory and is continuously increased 
to one to get to the kink case. As the parameter increases, modes in the trivial 
box evolve into modes in the kink background. Some of the low-lying modes in 
the trivial box become the bound states of the kink. Let us label these modes by 
the index b (for “bound”) and the remaining modes by c (for “continuum’’). (In the 
example of Fig. 4.1, b = 0, 1, and c is any integer except for 0, 1.) Then 


B= Bat 5 9 (œo) +5 | yP +m- Fm] (4.6) 


where we = ,/ p? + mi, and my denotes the mass of the w particles. In the limit 
L — œ, the sum over continuum states becomes an integral. 

The terms in Eq. (4.6) can be understood quite simply. The first term on the right 
is the classical kink energy, the second contains the excess quantum corrections 
owing to the zero point motion of the modes bound to the kink, and the third term 
is the excess energy in the zero point motion of the modes that are not bound to the 
kink. The wave numbers of the scattering modes in the background of the kink are 
denoted by p, while that of the modes in the trivial vacuum by k,,. 

In the trivial vacuum and when L — ox, the scattering states are plane waves, 
which are both energy and momentum eigenstates with k, = 27n/L. In the pres- 
ence of the kink, the scattering states are energy eigenstates but not momentum 
eigenstates and, in general, an incoming wave gets both reflected and transmit- 
ted. Without specifying the field theory, further progress is possible when the 
scattering potential, V’(¢,(x)), is reflectionless. This may seem very restrictive, 
but it holds for both the Z2 and sine-Gordon models and we assume it to be 
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true for the remainder of this analysis. Then, asymptotically, the scattering states 
behave as 


el(Px—a(P)/2) ag x = —L/2 > =œ (4.7) 
ellprta(P)/2) as x = +L/2 > +00 (4.8) 


where a(p) is a phase shift. Note that on multiplying by exp(ia(0)/2), the p = 0 
state can be chosen to be purely real at x = —oo. Since the scattering potential, 
V” (ġk(x)), is also real, this implies that the imaginary part of the wavefunction can 
be taken to be zero everywhere. Therefore w(0) = 0. 

The phase of the scattering states has a winding number given by the total phase 
change across the box. Since we have imposed periodic boundary conditions, the 
total phase winding, (pL + a(p))/2z, must be an integer. This quantizes p so that 


pL + a(py) = 2vr (4.9) 


where v € Z and we have denoted the vth wave-vector by py. 

Now that the scattering states in the soliton potential have been labeled by the 
integer v, and those when the potential vanishes by the integer n, correspondence 
must be drawn between v and n. To illustrate the problem, consider the p, = 0 
mode. As discussed above, a(0) = 0 and hence, from Eq. (4.9), v = 0 labels this 
mode. Further, this mode has the lowest energy of the continuum states. In the 
specific example of Fig. 4.1, this mode corresponds to the n = —1 mode in the 
trivial box since the n = 0, 1 modes have become bound and have dropped out of 
the set of scattering states. Therefore n = —1 corresponds to v = 0, in this example. 

With k, = 27n/L, we can write 


PL + A(py,) = 2nt =k, L (4.10) 
where 
A(py) = a@(py) + 27(n — v) (4.11) 


The shift in going from n to v is the change in the total winding of the phase as 
the potential evolves from the trivial box to the soliton potential (see Fig. 4.1). As 
long as there is no change in the relative ordering of the energy levels, the heirarchy 
of the energy levels is maintained, and the mapping between n and v is a constant 
shift. Since some of the low-lying states in the trivial potential have dropped out 
from the set of continuum states and have been converted into bound states, the set 
of integers n is partitioned into two subsets — one for the integers that lie above 
the would-be bound states and another for the integers that lie below the would-be 
bound states. The map from n in each subset to v is a constant shift but the shift is 
different in the two subsets. 
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Next we think ofn — v asa function of k. Inthe L —> oo limit, (n — v) is constant 
everywhere except at k = 0, meaning that the derivative of n — v with respect to k 
is a Dirac delta function at k = 0, 


a ff.ô(k oO) 4.12 

aE (k) + ae (4.12) 
To determine the coefficient of the delta function, let us denote by N, the number 
of states in the trivial potential that have dropped out as bound states in the kink 
potential. In a large interval n} — n_ (n+ is positive and n_ is negative), the cor- 
responding interval in v is smaller by N,, and hence the coefficient of the delta 
function is given by —2z7 N, 


dA da(k) 
a —N,27 (k) + dk 

For large momenta (and energy) the modes are unaffected by the deformation of 
the potential at the bottom of the well. Hence p, — kn in this region and A(p) > 0 
as |p| > œ. 

The phase shift A(p,) depends on the potential in the equation of motion for 
w(t, x), as in Eq. (4.3). As we explain below, the scattering potential created by 
the soliton background is non-perturbative. Therefore the phase shifts need not 
be small owing to factors of the coupling constant. However, note that A(p,,)/L 
is small as L — oo and we need only keep terms up to linear order in 1/L. 


Therefore 
ae ’) í 
JP +m, = (x - ” ) +m 


Kn ey +0(Z :) (4.14) 
L k2 + 


(4.13) 


Note that in the last line A(p,) has been replaced by A(k,,) since py = kn + O(1/L). 

We now want to express the energy of the kink in terms of renormalized param- 
eters. If we denote the renormalized mass of y by my R and the bare mass by my», 
then 


iy R = m? p = sm, (4.15) 


where dm,, denotes the quantum contribution of vacuum fluctuations to the mass 
of y, and ôm?, is due to the self-coupling of the field and hence is proportional to 
the coupling constant. 

The expression for the energy in Eq. (4.6) is valid to leading order in quantum 
corrections. The classical energy is inversely proportional to the coupling constant 
(e.g. Eq. (1.20)) and so the leading corrections are independent of the coupling 
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constant. Note that mẹ in the last two terms in Eq. (4.6) can be freely replaced 
by my r since we are only evaluating the lowest order (coupling constant inde- 
pendent) quantum correction to the energy and ôm?, is proportional to the coupling 
constant. Retaining only the terms that are of leading order in the coupling constant, 
expanding Ea in ôm}, and using Eq. (4.14), we get 


l l kn A(kn) 
E = EalmyRi Ar] + AEa +5) (oo) r (4.16) 
á ki + my 


where A Ea denotes the leading order change in Ea when replacing bare parameters 
by renormalized parameters. 
In the limit L — oo, the sum over n becomes an integral 


A ea 
> — J dk (4.17) 
2d 


Hence, 


1 N, 1 kA(k) 
E = Ealmy R; àr] + A Ea + 5 die = myr EF fa 


where we have made use of the fact that o = |k? + m? p > My R as kp X 
1/L > 0. 
On integration by parts 


kA(k +00 dA 
J pe RRO — J dkk? +m p—— (4.19) 
Je Fma Ze dk 


Since A(k) vanishes as k — +oo, the boundary term gives a finite contribution. 
The last term contains the derivative of A(k) and is given in Eq. (4.13). Therefore 
the final result is 


1 1 
E = Ealmy,xsdx] + AEa + 5 die aa [IKI ACK) It 


+7 fa etma t (4.20) 
— m? ; 
YR dk 
Our general calculation can be pushed a little further since, in one spatial dimen- 
sion, all divergences can be removed by normal ordering a “renormalized potential,” 
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Ve, which can be written in terms of the bare potential, V(@) [35] 


1 4A?\ æ 
Vr = exp an In me) age 


where A is a momentum cut-off, and m is the bare mass. The constant €o renor- 
malizes the vacuum energy, and is chosen so that the expectation value of the 
Hamiltonian in the ground state vanishes. For example, in 1¢* theory (Eq. (1.2)), 


| Vio) + €0 (4.21) 


2 1 n2, À 4 
Vr = [y (3yà — m^) + €o] + 5 6y — m^) + A (4.22) 
where 
1 i 4A?\ aoo 1 (7 dk (4.23) 
= — ln | —— = : 
Y Sir m2 8m —A af k2 + m2 


Then, the quantum correction to the mass is 5m? = 6Ay, while the quantum cor- 
rection to the mass of the excitations in the Z2 model is: 


5m), = 26m? = 12Ay (4.24) 


In the sine-Gordon model (Eq. (1.51)) 


Qa 


Va = 33 


and the quantum corrections to the parameters can be read off. 

Returning to the general expression in Eq. (4.21), the bare parameters occurring 
in V(@) can be chosen to absorb the cut-off dependent factors. Then the potential 
Vp is given entirely in terms of finite physical parameters. If the classical solution 
is found for the physical value of the coupling constant, denoted by Ar, then A Ea 
only depends on the correction to the mass term, 5m;,, 


[1 — eo VB cos(Bd)] + €o (4.25) 


EU [my R; Ar] S 


AEg = 
2My R 


? (4.26) 
where El denotes derivative of Ea with respect to the mass, my, r. At this stage, 
we are still left with the last two terms in Eq. (4.20) involving the phase shifts. 
However, there is no general prescription for finding the phase shifts, and each 
problem has to be dealt with individually. 

Equation (4.20) is our final general expression for the ground state energy of the 
quantized kink provided that the classical kink solution gives rise to a reflectionless 
potential. To make further progress one needs to find Ealmy, R; Ar], A Ec, œp, and 
the derivative of æ(k). These quantities are model specific and we shall find them 
in the à¢4 and sine-Gordon models in the next two sections. 
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Before going on to some examples, it is helpful to track the coupling constant 
dependence of the various terms in Eq. (4.16). We write the potential as 


m? 2 
V@)= =N + eES(ġ) (4.27) 


where m is the mass parameter, € is the small coupling constant, and S is some 
unspecified function of @, perhaps containing other parameters. The classical energy 
term in Eq. (4.16) is inversely proportional to the coupling constant. So the leading 
order correction is independent of the coupling constant. In the second term, ôm? 
is proportional to the coupling constant but E’, is inversely proportional to the 
coupling constant. Hence the product is independent of the coupling constant. Next 
we come to the coupling constant dependence of the energy eigenvalues and the 
phase shifts. The spectrum of excitations is found by solving for eigenmodes in 
the kink background. The kink background provides a potential with which the 
excitations interact. The important point here is that this potential is non-trivial 
even to zeroth order in the coupling constant. The vacuum expectation value of ¢, 
denoted by ġo, is found by minimizing V. Therefore 


4 2 
S(p e (4.28) 
Po € 
Then 
V” (po) = —m? + €S" (go) (4.29) 
and approximating $” (o) as S’(¢o)/¢o, 
E pel ~ ~m? (4.30) 
0 


Hence the scattering potential in Eq. (4.3) for the mode functions is independent 
of the coupling constant, and the phase shifts, w(k), are non-trivial even to zeroth 
order in the coupling constant. 

As we see in the specific examples given below, both smi, and the last sum in 
Eq. (4.16) are divergent. However, the divergences cancel, leading to a finite result 
for the energy. 


4.2 Example: Z kink 


We now find the energy of the quantized Z3 kink by evaluating explicitly the terms 
in Eq. (4.16). 
The classical energy piece is already known from Eq. (1.20) 
2/2 m? m? R 


EalMmy R; Ar] = a (4.31) 
R R 
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Then AE, is given by Eq. (4.26) and E/ is found by differentiating Eq. (4.31) 


2 
m? p 
ÀR 


E/ [my R; AR] = (4.32) 

The mass correction ôm?, arises owing to the interaction term 44/4 in this 
model. The calculation of 5m, is quite involved since it requires renormalization 
in a model with spontaneous symmetry breaking, which means that we should 
find the mass correction from the action in Eq. (1.5). Then there are both cubic 
and quartic interactions. This calculation can be found in quantum field theory 
textbooks, for example [119]. The end result is 


3A dk 
sm, = E J (4.33) 


The integral in Eq. (4.33) is divergent. However it is only one term in the expression 
for the quantum kink energy in Eq. (4.18). In particular, the last term with the 
phase shifts is also divergent, but the quantum kink energy is finite since the two 
divergences cancel. Note that we can replace my by myr in the final integral since 
we are only evaluating the leading order correction. 

Next consider the terms in Eq. (4.16) that involve the spectrum of fluctuations 
about the classical kink. To find the spectrum, substitute Eq. (4.1) in the field 
equation, Eq. (1.4), and expand to lowest non-trivial order in y. This was already 
done in Section 3.2 and we now summarize the results 


wo = 0, Xo = sech?z 
3 
w = he X1 = sinhz sechz 
My <w < œ, Xk = e!13 tanh?z — 1 — w7k* — i 3wk tanhz] 


where z = x/w = myx /2, and the dispersion relation is 
o =k +m, (4.34) 


Note that the eigenvalues w% are independent of the coupling constant because A 
does not occur in Eq. (3.10) if my = J2An is held fixed. However, this statement 
is only true to leading order in A because the mass parameter, the kink width, and 
indeed the form of the kink solution get modified owing to quantum corrections, 
and induce à dependence in the spectrum. Since we are only working to leading 
order in quantum corrections, the mass parameter my entering Eq. (3.10) and the 
definition of the kink width, w, are the same as my r. 
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The next step is to impose periodic boundary conditions with period L — oo on 
the scattering state. For this we find the asymptotic behavior of xx 


xn > e(2 — wk? Fi 3wk) x expfi(kz + a(k)/2)}, z— +00 (4.35) 


from which the phase shifts follow 


3wk 
= 2 tan™! | ——_—__ 4, 
a(k) an È =| (4.36) 
Hence 
da wk? +2 
— =-6 4.37 
dk Uke + Diwk2 + 4) oH) 
and 
6 
A(k) > —, |k|—> œ (4.38) 
wk 


Now we combine all the terms in Eq. (4.20) 


E m} R 7 v3 a 3 ee 3m, g D dk 1 
= ary ee RTZ 
Jie Oe VA Ma d eg m? pK +m 2/4 
(4.39) 


The last integral is done easily yielding the final result for the kink mass with 
leading order quantum correction 


p=" (3 1 ) Z 


De x 2/3) 2 
3 
My R 
= — 0.33 4.40 
Re My R ( ) 


Note the minus sign in front of the quantum correction to the energy. In Section 4.5 
we show that this is a general feature. 


4.3 Example: sine-Gordon kink 


To quantize the sine-Gordon kink of Section 1.9, we follow the same procedure as 
for the Z2 kink. The mode functions now satisfy 


2 2 
2 _ @ 
aya + tanh? X — Dy = ae (4.41) 
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where X = m,,x. The kink solution, from Eq. (1.52), is 


Px = tan (e) = sill E (e”"»*) 


(4.42) 
where A = af”. The classical energy (Eq. (1.55)) is 


va _ iy 
Eyc.ct = age = ra (4.43) 
The spectrum has only one bound state, the translational zero mode given by 
d 2m? 
oa, =0, Wo sa ele 


= = ——sech(myx 4.44 
ae Z pp honya) (4.44) 
The scattering state with wave-vector k can be written quite generally in terms of 


hypergeometric functions (see [113], Vol. II, Section 12.3, or Appendix C) 


1 3 e%* 
» = KX . . . 

We = N(cosh X)' F (ie = 1, —1kK + 7 + 5! = ikl) (4.45) 
where N is a normalization factor and x = k/my corresponds to the wave-vector. 
The phase shifts are found by taking the asymptotic forms of Eq. (4.45) 
Ye > Ne"¥, X > oœ 
> Nelit+2eikX X o (4.46) 

where tan 0 = x. Hence the phase shift is 


k 
a, = m — 2tan! (5) (4.47) 
My 
Therefore 


(4.48) 
v 
At large |k|, A(k) (as needed in Eq. (4.20)) is given by 


A(k) = ana 


, |kl —> œ (4.49) 
To find AZ, occurring in Eq. (4.20), we can use the renormalized potential in 
Eq. (4.25). The parameter 6, which occurs in the argument of the cosine function, 
is taken to be the physical (renormalized) value, while 


2 
dosak (1 + Ey) (4.50) 
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to leading order in 6. The subscripts refer to bare and renormalized quantities and 
y is defined in Eq. (4.23). Therefore 


+A dk 
AEa = R i / 
2m J-a [pay mi, 
Finally, with $` wp, = 0, we can put together all the various terms in Eq. (4.20) 
to get 


(4.51) 


3 
8m, R MyR tA dk My R a My R 


+A dk 
+0 J 
ÀR 2m J-A k? +m? T 2m J-A [k2 +m, 


AR T 
8m? 
= a — 0.32my.R (4.52) 


Once again the quantum correction is negative and, coincidentally, quite close to 
the Z value (see Eq. (4.40)). 


4.4 Quantized excitations of the kink 


So far we have only calculated the quantum correction to the mass of the kink in 
its ground state. Now consider the excited states of the kink. 

As in the second quantization of a free quantum field theory, particle creation 
and annhilation operators are introduced for each of the excitation modes of the 
kink. As we shall see, this is straightforward except for the zero mode. The end 
result is a procedure for doing quantum field theory with both particles and kinks 
included in the spectrum of states. Here we only give some introductory remarks. 
For a more extended discussion see [67, 126]. 

Let us denote the bound state mode functions by F;(t, x) and the scattering 
mode functions by f(t, x). The t dependence is of the form exp(—iw;t) where w; 
is the frequency of the bound or scattering mode. Then the second quantized field 
operator is 


p = pa) + Y [eo File. x) + FEE, 2] + YO [Ge felt. x) + a EO 2) 
b k 


(4.53) 
where œx is the classical kink solution, cl, cp are creation/annihilation operators 
for the bound states, and similarly al, dx are creation/annihilation operators for the 
scattering states. Now, for the zero mode, w = 0 and Fo(t, x) = Fy (t, x). Therefore 
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the zero mode contribution to the sum is 
[ĉo + ĉi Fo(x) (4.54) 


Since cp and ch are only present in the combination cg + ci let us define bọ = Cp + ĉi 
which is then the annihilation operator for the zero mode. However, note that 
bi = bo and so [bo, bi] = [bo, bo] = 0: the zero mode is classical as the operator 
bo commutes with all other operators. This is to be contrasted with [a,, a|] = 
27 5(k — p). 

Just as the translation mode is a bosonic zero mode, there can also be fermionic 
zero modes that we discuss in Chapter 5. In that case, the creation and annihilation 
operators satisfy anticommutation relations leading to {bo, bi} = 0. This relation 
has the remarkable consequence of leading to fractional quantum numbers as we 
discuss in Chapter 5. 


4.5 Sign of the leading order correction 


A striking feature of the leading order quantum corrections to the energies of the 
Z and sine-Gordon kink is that they are negative. In other words, quantum effects 
reduce the mass of the kink. A variational argument [104] (Coleman, S., 1992, 
private communication) shows that this observation holds true quite generally in 
one dimension. ! 

Let the Hamiltonian of the 1 + 1 dimensional system be 


H= fon = fa [Ho + V()] (4.55) 


where ¢ is a scalar field, 


H z 241% oy (4.56) 
TS 357 7 X . 


and z is the canonical field momenta. Written in this way, the parameters entering 
the Hamiltonian are bare parameters and subject to renormalization. In one spatial 
dimension, however, it can be shown that [35] 


H = Nm [Ho + Vr] (4.57) 


where N,, denotes normal ordering with respect to free particles of mass m, and 
the renormalized potential is (Eq. (4.21)) 
p 1 i 4A? \ d? 
= ex n 
R Pl 8x m? J do? 


1 The conclusion may not hold if the model also contains fermionic fields. 


| V$) +e (4.58) 
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where A is an ultraviolet momentum cut-off and € is a constant to be chosen such 
that (0||0) = O where |0) is the true ground state of the model. 

The energy of the kink, including the contribution of quantum fluctuations in the 
ground state, is 


E = (O|H[¢x + WI10)x (4.59) 


where |0), denotes the vacuum for the quantum fluctuations, y, around the classical 
one kink state Ox. 
Straightforward manipulation now gives the quantum correction to the kink mass 


E — Ear = (OH [k + W] — AL¢«]19)x 
= fo k(OINn (Hol gx + Y] — Holek] + Vrløk + Y] — Valdx)10)x 


where Ee,r is the energy of the classical solution obtained with the renormalized 
potential, Vp. Next we use the variational principle, which states that the ground 
state energy of a system is minimized in its true ground state, and the expectation 
of the Hamiltonian in any other trial state gives an upper bound to the ground state 
energy. If we denote the perturbative vacuum state — the state with zero particles of 
mass m — by |0, m), then 


E < Ear + fo (0, m|Nn(CHoldx + Y] — Hold] + Vealox + Y] — Valox])|0, m) 


= Ear 


The last line follows since there are no y independent terms in the expectation 
value under the integral,” and the annihilation operators of w occur to the right 
owing to normal ordering and annihilate the trial vacuum state. 

Note that Ear is the energy of the classical solution found by minimizing Hr[@], 
i.e. the Hamiltonian in Eq. (4.55) but with the potential given in Eq. (4.58). Since the 
true ground state of the system is not known, the constant € is not known either. The 
potential Vp can be minimized, but there is no guarantee that the minimal value of 
Vp will be zero. Therefore Ea, r might get an infinite contribution from integrating 
min(Vp) over all of space. Then the variational bound E < Ea, is not very useful. 
However, we do know the value of € to lowest order in the coupling constant and 
this is precisely so that (0, m| H|0, m} = 0. This coincides with choosing € so as to 
make min(Vp) = 0. Hence the bound 


E < Ear = Ea (4.60) 


2 To see this, note that the expectation value vanishes if y = 0. 
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where Ea denotes the classical energy without any quantum corrections, is mean- 
ingful to leading order in perturbation theory and it provides us with the completely 
general result that the lowest order correction to the soliton energy is negative. 


4.6 Boson-fermion connection 


A bosonic field, ¢, in quantum field theory satisfies the equal time commutation 
relation 


lox, 1), $, D] = ê — y) (4.61) 


Alternatively, a fermionic field, y, satisfies the anticommutation relations 


tala, D, Viiv, D) = 8 — Yab (4.62) 


where a, b = 1, 2 label the two components of the spinor in one spatial dimension. 
It is remarkable that one can construct explicitly a fermionic field w satisfying 
Eq. (4.62) in terms of a bosonic field ¢ that satisfies Eq. (4.61) [108]. 

The connection between Wa and ¢ is 


pa = Cres, W(x) = iC :e?™: (4.63) 


where the c-number C is defined in terms of a mass parameter jz and another cut-off 
parameter, €, 


1/2 
C= (=) elt/8e (4.64) 


The operators P contain a free parameter 6 and are defined by 


2a [* ; ip 

Px) = -i-> dé oE) F > O@) (4.65) 
P J-co 2 

The symbol :: in Eq. (4.63) denotes normal ordering with respect to the mass u. 
This means that the field @ is to be treated as a free field with mass parameter 
u and the quantum operator, ¢, is expanded in terms of creation and annihilation 
operators that create and destroy particles of this free field theory. A normal ordered 
operator contains various products of creation and annihilation operators with the 
annihilation operators always occurring on the right. It is understood that the in- 
tegral in Eq. (4.65) is cut off at large € by a factor exp(—eé). Note that normal 
ordering is a symbol and should be treated carefully — normal ordering of strings 
of operators should be done prior to commuting operators that occur within the 
string. 


66 Kinks in quantum field theory 
To check if Eq. (4.62) is satisfied for x 4 y, we use the identity (see Appendix D)? 


SCREEN Engh ee? i= ea) se? et (4.66) 


where A and B are any two operators that can be written as a linear sum of terms 
involving only creation or annihilation operators 


ASAT A, B = B* + B~ (4.67) 


The commutators [A*+, B7] and [B*, AT] are assumed to be c-numbers. Insertion 

of Eq. (4.66) in Eq. (4.62) gives the commutation relation in Eq. (4.61) for x Æ y. 
It is harder to check that the commutation relations in Eq. (4.61) hold when 

x = y. Since products of quantum operators at the same point are singular, the 

commutator must be evaluated at two different points in space, x and y, followed 

by the coincidence limit y —> x. We now outline the scheme employed in [108]. 
We want to check the anticommutation relation 


{Walx), VO = Zê — y) (4.68) 


where the constant Z, possibly infinite, has been introduced in recognition of the 
fact that the fields get renormalized. Rather than check Eq. (4.68), we can check 
the equivalent commutation relation 


U E), VO) =- G + Pi go z5) W(x)d(x — y) (4.69) 


where the current j” has been regularized using point-splitting and is defined by 


, | p 
j*(x) = lim TE Izer uE — y Py WO) + F" y) 
(4.70) 
where o is a regularizing parameter and F“(x — y) an unspecified c-valued func- 
tion. The y-matrices are defined by the algebra 


iy", y” aoe  yř=iy?y! (4.71) 


where g,,, = diag(1, —1) is the two-dimensional Minkowski metric. An explicit 
representation of the y-matrices is given in Eq. (5.15). In Eq. (4.69), €”” is the 
totally antisymmetric tensor. 

First the current j“ is evaluated with w, as given in Eq. (4.63). The evaluation 
requires 


[Ot (x, t), 6 (y, t)] = Ai — y)? — (dt + ie)’] (4.72) 


3 In the literature it is sometimes incorrectly stated that the identity e4+8 = el?-41/2e4e8 


is being used. 


e” (no normal ordering) 
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where dt = th — tı and A, is the propagator. For small x — y 


1 
Ay =- Inju? — y? — (dt + ie}] + Ox — yy’) (4.73) 


By differentiating Eq. (4.72) we can also obtain the commutators of time derivatives 
of ¢t and ¢~. These appear in the evaluation of j” since y is defined in terms of 
¢ in Eq. (4.63). 

The result for j” is singular in the limit y —> x except for a single choice of the 
regularizing parameter, o , occurring in the definition of j“. This single choice is 


2 2 
rae (1 2 =) (4.74) 


With this value of o, the commutator in Eq. (4.69) can be verified. Thus the 
w operator indeed satisfies the anticommutation relations of a fermionic field. 
Furthermore, the current can be explicitly calculated, leading to 


j= FP way (4.75) 
27 


To summarize, given a quantum scalar field in 1 + 1 dimensions, it is possi- 
ble to construct a fermionic field from it via the relation (4.63). Starting with a 
fermionic field, a bosonic field may be constructed from it via Eq. (4.75). Note that 
the transformations from bosons to fermions and vice versa hold at the quantum 
operator level and not just at the level of expectation values. Further, they hold 
for any choice of interactions in the bosonic or the fermionic model. However, in 
the case when the bosonic model is the sine-Gordon model, the fermionic model 
obtained by transforming to the fermionic variables is another well-known model, 
namely the massive Thirring model as we now describe. 


4.7 Equivalence of sine-Gordon and massive Thirring models 


The sine-Gordon model is given by the Lagrangian (Eq. (1.51)) 


1 2 Q 
Lsa = 50L — pz (1 — cos(Bo)) (4.76) 
while the massive Thirring model is 
Lar = i Jy — mb — Say Prv (4.77) 


where w is a two-component fermionic field. 
In [34] (also see [35]) it is shown that the sine-Gordon model does not have a 
well-defined ground state for B* > 87x. To clarify what this means, consider the 


68 Kinks in quantum field theory 


simple example of a free field theory 


1 ô 
Lire = 5 (ud) — 5? (4.78) 


This model has a well-defined ground state only in the range 6 > 0. The model does 
not have a well-defined ground state for ô < 0. Similarly the sine-Gordon model 
only has a ground state for a definite range of parameters, though the reasons are 
much more subtle.* The sine-Gordon only has a well-defined ground state if the 
parameter B? is restricted to lie in the interval (0, 877). 

In the range, 0 < 6? < 87, there is a one-to-one mapping between vacuum 
expectation values of a string of operators in the sine-Gordon model to those in 
the massive Thirring model. This means that any vacuum expectation value in the 
sine-Gordon model has a “corresponding” vacuum expectation value in the massive 
Thirring model. This strongly suggests that the two models are equivalent, even at 
the operator level [35]. 

As we have seen in the last section, there is indeed a two-component fermionic 
field, y, that can be constructed from a bosonic field ¢ (Eq. (4.63)). In [108] it was 
shown that y also obeys the equations of motion of the massive Thirring model 
if the bosonic field ø obeys the equations for the sine-Gordon equation with the 
coupling constant g written in terms of the coupling constant 8 as 


ae eee (4.79) 


Note that when the sine-Gordon model is weakly coupled (small 6), the massive 
Thirring model is strongly coupled and vice versa. Hence the sine-Gordon model 
and the massive Thirring model are completely equivalent as quantum field theories 
but one is a better description at small 6 (large g) and the other at large 6 (small g). 

What has the equivalence of the sine-Gordon and massive Thirring models got 
to do with kinks? Consider the commutation relations between @ and y. Using 
Eq. (4.63) and the identity (see Appendix D) [A,:e? : ] =: [A,e?]: with A = 
(y) and : e? := y we find 


2 
[o(y). Wx) = Fe (y <x) (4.80) 


ILO) Wx) =0, (> x) (4.81) 


Now consider the action of w(x) on an eigenstate, |s) of the field operator ¢. Let 
us choose this eigenstate to be such that 


ols) =0 (4.82) 


4 For example, in contrast to the model in Eq. (4.78), the classical sine-Gordon model has well-defined global 
minima for all values of the coupling constant £. 


4.7 Equivalence of sine-Gordon and massive Thirring models 69 


If we write |s’) = y(x)|s), the relation in Eq. (4.80) gives 


/ 20 1 
dys) = g" b <x) (4.83) 


and Eq. (4.81) gives 


pls) =0, >x) (4.84) 


Therefore the state obtained after action by Y (x) is one where the value of ¢ is 27x /B 
for y < x andOfor y > x. In other words, the field w(x) creates a step-change of @. 
The step-function profile is viewed as a “bare kink” which gets dressed by quantum 
effects that produce a smooth kink profile with some finite width. So the field w(x) 
is the creation operator for a (bare) soliton at location x. In the Thirring model, the 
field y(x) is interpreted as the creation operator for a fermion located at x. Hence 
the sine-Gordon kink is identified with the fermion in the massive Thirring model. 
The topological charge on a sine-Gordon kink is 


Qk = [ow = [ost (4.85) 


where the fermionic current is defined in terms of the bosonic current in Eq. (4.75). 
Therefore the fermion in the massive Thirring model carries the topological charge 
of the sine-Gordon kink. In other words, the kink of the strongly coupled sine- 
Gordon model is better described as a weakly coupled fermion of the massive 
Thirring model. Here we see the duality between particles and solitons. 

Can we also interpret the bosonic particles of the sine-Gordon model in terms of 
“solitons” of the massive Thirring model? The massive Thirring model only contains 
fermions, and classical solutions of the Dirac equation do not have the interpretation 
of solitons. This is because the fermionic fields anticommute and fermions obey the 
Pauli exclusion principle. Instead a classical solution of the Dirac equation is a state 
that one (and only one) fermion can occupy. However, there can be bound states 
of two or more fermions since the force between a fermion and an antifermion is 
attractive for g > 0. A bound state of two fermions can be shown to correspond to 
a particle of the sine-Gordon field ¢. If the fermions in the weakly coupled massive 
Thirring model have mass m, then the bound state energy is approximately 2m 
since it involves two fermions. However, the binding energy decreases (becomes 
more negative) with increasing interaction strength, g, and eventually the bound 
state becomes lighter than a single fermion. At this stage, a suitable description of 
the system is in terms of the bound state being the fundamental degree of freedom 
as in the sine-Gordon model. 

The bound state of two massive Thirring fermions is also a bound state of 
two sine-Gordon kinks i.e. a breather. Hence it should be possible to interpret 
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the breather as a particle in the sine-Gordon model. This is seen to be true when the 
breather is quantized [38—41, 35]. Then, to lowest order, the energy levels of 
the quantized breather are equal to the mass of one, two, three, etc. particles of 
the sine-Gordon particle. 


4.8 Z, kinks on the lattice 


Lattice field theory provides another tool to probe the quantum nature of solitons 
and, in particular, the variation of mass with coupling constant. 

The starting point is the action for the Z model defined in Eq. (1.2). The action 
is to be inserted in the Feynman path integral, which can then be used to find 
expectation values for any quantum operator. In the Feynman path integral, it is 
necessary to integrate over field configurations, and this is done numerically on a 
discretized Euclidean space-time. The reader is referred to the lattice literature for 
details [37, 112, 141]. Here we shall give the results relevant to the Zz kink. 

The mass of a Z; kink is defined as the expectation value of a suitable operator 
defined on the lattice in the limit that the lattice spacing, a, goes to zero. One 
important issue is that there are several different candidate operators on the lattice 
that all go to the correct limit as a — 0 and, in practice, it is not possible to take the 
limit all the way to a = 0. At best, the numerical analysis gives the expectation of 
the operator on the lattice for several different values of a and then some scheme 
must be found for extrapolating the results to a — 0. In [32], the authors evaluate 
the mass of the Z2 kink using two different lattice operators. The results are shown 
in Fig. 4.2. We note that the kink mass decreases monotonically as the coupling 
constant increases and remains bounded by the classical mass. At a certain coupling, 
the kink mass goes to zero, and the kink, not the @ quanta, is the lightest degree of 
freedom in the model. 

The mass of the sine-Gordon kink has been calculated analytically for a range 
of parameters in [156] (also see [35]). 


4.9 Comments 


Several researchers have taken alternate paths to studying quantized kinks. In 
supersymmetric theories there is greater control over quantum corrections and the 
mass can, in some cases, be evaluated exactly [51]. Alternate methods to quantize 
supersymmetric kinks have also been developed in [66]. Variational methods to 
study the A#* theory have been developed in [49]. The scattering of kinks in classi- 
cal and quantum theory has been studied in [153]. Kink masses and scattering have 
also been calculated in [132] using the Hartree approximation. 
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Figure 4.2 The figure shows how the mass of the Z2 kink depends on the inverse 
coupling constant, By = 1/A9, where 49 = 6Aa? is the coupling constant in the 
discrete theory on a 48 x 48 lattice and a is the lattice spacing. (The factor of 6 


is due to our choice of 1/4 in the à term in Eq. (1.2) as opposed to 1/4! in [32].) 


The lattice mass parameter, ro = —m?a?, is held fixed at rọ = —2.2. From the 


plot we see that the classical value of the kink mass is larger than the quantum 
value. The one-loop corrected mass (see Section 4.2) and the mass found by using 
two different choices of the lattice mass operator are also shown. The kink mass 
vanishes at By = 0.0804. [Figure reprinted from [32].] 


The construction of fermion operators from boson operators and vice versa has 
been discussed and used extensively in condensed matter applications under the 
name of bosonization. A review, in addition to an historical introduction, may be 
found in [44]. Finally, the technique of bosonization has also been applied to thermal 
systems in [69]. 


4.10 Open questions 


1. The quantum corrections to the Z2 and sine-Gordon kinks were calculated explicitly using 
the phase shifts. However, the phase shift approach only works if the potential U(x) is 
reflectionless. What are the conditions necessary for a potential to be reflectionless? Are 
reflectionless potentials always in factorizable form (see Section 3.3)? (The example of 
a step-function potential shows that the converse is not true.) 
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. We have shown that the leading order quantum correction to the kink mass is always 


negative. Can this statement be generalized to all orders? Can one show that the mass 
of a kink goes to zero in the strong coupling limit? Or perhaps that it is monotonically 
decreasing as a function of increasing coupling constant? 


. If the Z3 kink at strong coupling is to be viewed as a particle, then the particle must obey 


unusual statistics because two kinks cannot be next to each other. Discuss this statistics 
and its implications for the dual model. 


. From the SU(5) x Z2 example we learned that a classical kink may be embedded in 


many different ways in “large” models. On quantization, do the different embeddings 
correspond to distinct degrees of freedom? 


. Does the addition of fermionic particles change the conclusion that quantum corrections 


always reduce the energy of a kink? 


. For the sine-Gordon model we have explicitly seen that there is a relation between kinks 


and particles. It seems reasonable that the connection holds in other models too. In3 + 1 
dimensions, we could expect the connection to exist between magnetic monopoles and 
observed particles (e.g. [162, 103]). Construct a model that has families of solitons, 
similar to the electron, muon, and tau families observed in Nature (see [122]). 


. In Section 3.1 we have discussed the existence of quasi-breather solutions called 


“oscillons” in the Z model. Can quantum oscillons have an interpretation that is similar 
to quantized breathers as discussed at the end of Section 4.7? 


5 


Condensates and zero modes on kinks 


In this chapter we study the effect of a kink on other bosonic or fermionic fields that 
may be present in the system. Under certain circumstances, it might be energetically 
favorable for a bosonic field, denoted by x, to become non-trivial within the kink. 
Then we say that there is a “bosonic condensate” which is trapped on the kink. 
On a domain wall, the condensate has dynamics that are restricted to lie on the 
world-sheet of the wall. 

The situation is similar for a fermionic field though there are subtleties. For a 
fermionic field, denoted by y, the Dirac equation is solved in the presence of a 
kink background made up of bosonic fields. This determines the various quantum 
modes that the fermionic excitations can occupy. In several cases, there can be 
“zero modes” of fermions in the background of a kink and this leads to several 
new considerations. (Fermionic zero modes were first discovered in [27, 84] in the 
context of strings.) In addition to the zero mode, there may be fermionic bound 
states. The high energy states that are not bound to the wall are called “scattering 
(or continuum) states.” 

A difference between bosonic and fermionic condensates is that bosonic solu- 
tions can be treated classically but fermionic solutions can only be interpreted in 
quantum theory. For example, while there may be a bosonic solution with x = 0, 
the solution w = 0 of the Dirac equation has no meaning because this solution 
is not normalizable. Solutions of the Dirac equation are only meant to supply us 
with the modes that fermionic particles or antiparticles can occupy, and as such 
are required to be normalizable. It is a separate issue to decide if the modes are 
occupied or not. A mode contributes to the energy of the soliton only if it is occu- 
pied. This is quite different from the bosonic case in which there can be a classical 
condensate, on top of which there are modes that may or may not be occupied 
by bosonic particles. Fermions can form a classical condensate only after they 
have paired up to form bosons (“Cooper pairs”), and this leads to superfluidity or 
superconductivity. 
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Fermionic zero modes can lend solitons some novel properties such as fractional 
quantum numbers (see Section 5.3). 


5.1 Bosonic condensates 


Consider the model 


1 
L = Lilo] + 5.x) — UC, x) (5.1) 


where L,[@] is the Lagrangian that leads to a kink solution in ¢. For example, Ly 
can be the Lagrangian for the Z3 or sine-Gordon models discussed in Chapter 1. x 
is another scalar field that interacts with @ via some general interaction potential 
U(¢, x). Note that U(¢, x) does not contain any terms that are independent of x — 
those are included in the potential, V(@), that occurs in Ly. As an example, we 
could have 
2 
U(b. x) = — hx + Exh 4 92? (5.2) 
2 4 2 
We are assuming that the parameters in the model are chosen so that the minimum 
of the full potential, V + U, is at o Æ 0 but x = 0. This requirement also excludes 
terms that are linear in x (e.g. o”x). 
In the fixed background of the kink, x satisfies the classical equation of motion 


07x — 32x + Uy (dx(x), x) = 0 (5.3) 


where U, denotes the derivative of U with respect to x and øx is the kink solution. 
Far from the wall, the lowest energy solution is x(00) = 0. 

A solution to Eq. (5.3) is x(x) = 0 and the energy of this solution is equal to the 
kink energy in the model Ly. However, the trivial solution may not be the one of 
lowest energy. To show that a lower energy solution exists, we need only show that 
the trivial solution, x = 0, is unstable. Then we consider linearized perturbations 
of the form x = cos(@t) f(x) around the trivial solution. Inserting this form into 
Eq. (5.3) leads to the Schrödinger equation 


—82 f + Uy, (bx(x)) f =o f (5.4) 


where U,, denotes the second derivative of V with respect to x. If this equa- 
tion has solutions with œ? < 0, it implies that there are solutions for x on the 
kink background that grow with time as cosh(+|@|t), denoting an instability of 
the state with x = 0. This means that the solution with least energy must have a 
non-trivial x configuration. The lowest energy x configuration is non-zero inside 
the kink and vanishing outside and is called a “bosonic condensate” (or simply 
“condensate’’). 
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5.1.1 Bosonic condensate: an example 


A simple example in which there is a bosonic condensate on a Z> kink can be found 
in the model of Eq. (5.1), or explicitly, 


1 1 à m? A o 
L= 8 2 z(a a RD 2N2 X 2 UX 4 Dig 55 
EO) PaO a aK aK = OX (5.5) 
Ignoring the condensate field x, the kink solution is 
x 
ox = ntanh (—) (5.6) 
w 
and the Schrödinger equation corresponding to Eq. (5.4) is 
—3? f +|- m? + on tanh? X] f =? f (5.7) 


where X = x/w. 

With on? > m?, we see that the Schrödinger potential is asymptotically positive, 
and hence f (+00) = 0. This is consistent with the requirement that x not have a 
vacuum expectation value. At the origin, U,, < 0, and hence the Schrödinger 
potential is a well that is centered at the origin. Since a potential well in one 
dimension always has a bound state [139], it follows that there is at least one bound 
state for x. For a deep enough well i.e. large enough m? , the bound state has negative 
energy eigenvalue (œ? < 0), and the trivial solution x = 0 is unstable. Hence there 
is a range of parameters (m?) for which a x condensate exists. 

To determine the range of m? for which there is an instability, consider the 
critical case when there is a zero eigenvalue solution, fo, of Eq. (5.7). Then we can 
write 


o 2 o 2 
—3? fo + 2 tanh? X — 1] fy = [m z a] ifs (5.8) 


This is exactly the same form as Eq. (3.8), together with the potential in Eq. (3.10), 
provided we identify 3A with o, and œ? with the term within square brackets on 
the right-hand side. Since the lowest energy eigenvalue is zero for Eq. (3.8), there 
is a zero eigenvalue for Eq. (5.7) if 

2_ om 


kk = 


r= (5.9) 


For a larger value of m? , Eq. (5.7) has a negative eigenvalue, signaling an instability. 
Therefore a condensate solution exists in the range 
2 
on 2 


3 <m < on (5.10) 
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The exact profile for the condensate can be found by solving the full coupled 
equations of motion for @ and x. This includes the non-linear terms in x and the 
back-reaction of the condensate on the kink, and, in most cases, has to be done 
numerically. Let us denote the solution obtained in this way by (k(x), xo(x)). 
Then 


X(t, x,y,z) = xo(x) cos(wt — kyy — k:z + 0o), w= NLS +k? (5.11) 
where 6p is a constant, is also a solution. The reason is simply that 
(3? — a5 — 87) cos(wt — kyy — kzz + o) = 0 (5.12) 


The trigonometric form of the solution in Eq. (5.11) was chosen so as to obtain 
a real solution. An identical analysis in the case when x is a complex field leads to 


x(t, x, y, Z) = Xox ete’ Hy-kz+ Hp) (5.13) 


The solution represents waves propagating in the (k,, k,) direction in the plane of 
the domain wall. 


5.2 Fermionic zero modes 


Fermionic fields may be coupled to the kink via terms that respect the discrete 
symmetries in the bosonic sector that are responsible for the existence of the kink. 
In the case of the Z} model, the coupling can be a Yukawa term and the Lagrangian 
may be written as 


L=Lo+iv py —govy (5.14) 


where Lẹ denotes the scalar part of the Lagrangian, J = y“3,, g is the coupling 
constant, y is a four-component fermionic field, and y” are the Dirac matrices 
that satisfy {y”, y”} = 2g”” with g,,, = diag(1, —1, —1, —1) being the space-time 
metric. The explicit representation of the Dirac matrices that we adopt is 


0 _ 0 1 iE 0 o! 
4 a pge- eg 26 i) 


wherei = 1, 2, 3 (also sometimes written asi = x, y, z) and the Pauli spin matrices 
are defined as 


ı_/0 1 2_(0 =i afd. 0 
oa a = <a E 1) (5.16) 


The Yukawa interaction term in the model in Eq. (5.14) respects the 6 —> —@ 
symmetry of the Z2 model provided we also transform the fermion field by 
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y > wv’ = y>wv where 
. -1 0 
ysiyy'yy = ( p 1) (5.17) 


This can be seen by using the properties (y°)' = y5, {y>, y“} = 0 and (y5? = 1.) 
If d(x) denotes the kink solution, the Dirac equation in the kink background is 


ify — sox(x)y =0 (5.18) 


Let us first try and solve Eq. (5.18) explicitly. Recognizing that øg does not 
depend on f¢, y, and z, we write the ansatz 


Y = ft, y,2&() (5.19) 


where f(t, y, z) is a function while &(x) is a four-component spinor. With this 
ansatz, the Dirac equation separates 


iy 0q f = —Y“kaf (5.20) 
iy" ðr — gkg = +y “kag (5.21) 


where y“%k, is the constant matrix of separation and the index a runs over t, y, Z. 
Requiring that the fermion be localized on the wall, we get the boundary conditions 


E(+œ)=0 (5.22) 


These boundary conditions are valid only for bound states. If we wish to consider the 
scattering of fermions off a domain wall, we would choose incoming and reflected 
plane waves at x = — 00. 

The Dirac equations have an infinite number of solutions, corresponding to all 
the fermion eigenmodes in the domain wall background. These include fermionic 
bound states and scattering states. There is one state, however, which is novel 
because it leads to some very interesting properties of the soliton, described in the 
sections below. This state is the one with zero energy eigenvalue, also called the 
“zero mode.” 

Equation (5.20) can be solved 


f = explikax®) = exp(i(wt — kyy — kzz)) (5.23) 


Zero energy is obtained by setting w = 0 = k, =k, and then f = 1. Let us first 
look at this case (ka = 0). 

Multiplying Eq. (5.21) by iy* we see that iy*é satisfies the same equation of 
motion as €. Therefore if € is a solution, then so is iy“. Hence solutions to the 


1 The Yukawa term does not respect the @¢ > $ + 27/6 symmetry of the sine-Gordon model and hence our 
discussion of fermion zero modes cannot be used for that case. Nor do we consider the case of fermions with 
Majorana mass terms [147]. 
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Dirac equation come in distinct pairs unless £ is an eigenstate of iy*, in which case 
the two solutions € and iy*é are not distinct. The zero mode solution is found by 
choosing é to be an eigenstate of iy* 


iy*E =c& (5.24) 
and, since (iy*)* = 1, we must have c = +1. The £ equation now becomes 


dx = cog (5.25) 


and the solution is 
E(x) = &(0) exp ef aoar | (5.26) 


Assuming $,(+00) > 0 and g > 0, the boundary conditions in Eq. (5.22) are only 
satisfied if c = —1. (The boundary condition at x = —œ is also satisfied provided 
dx(—oo) < 0.) Therefore the zero mode solution is 


&(x) = §(0) exp |-« | aoda | » (g>9) (5.27) 


If (+00) < 0 and g > 0, the solution is obtained by choosing c = +1. 
To determine (0), we solve the eigenvalue equation iy*&(0) = —&(O) and find 


a 


B 
if 


1a 


(0) = (5.28) 


where œ, 6 are any complex constants. Therefore there are two basis zero modes 
(with coefficients œ and £) and the general zero mode is a linear superposition 
of these two modes. The constants, œ and £, can be fixed by normalizing the 
wavefunction. 

Next consider the case with ka Æ 0. Then Eq. (5.27) is still a solution to Eq. (5.21) 
provided k,y“&(0) = 0. By explicitly substituting the y“ matrices and &(0), this 
leads to the two equations 


kya + ilo + k) =0 (5.29) 
il@ — ka — kf =0 (5.30) 


A solution for a and £ exists only if 


w= k +k (5.31) 
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which is the dispersion relation for a massless particle (see Fig. 5.1). With this 
relation, the solutions fix the ratio of œ and £ to obtain 


væ + k 
re VESS T hada! sgn(ky )iv@ — kz 
2/@ —sgn(ky)./@ — k; 
ivo + k; 
where N is a normalization constant where sgn(ky) = ky/|ky|. 
So far we have not specified the exact form of the kink profile x and Eq. (5.32) 
holds for any model in which the Yukawa interaction term respects the symmetries. 


Next, as an example, we use the solution for the Z kink (see Eq. (1.9)). Then the 
integral over ¢, can be done explicitly to yield 


(5.32) 


vo + k: 


Neiot—kyy-kz) x \78Vv27 | sgen(ky)i/@ — k 
y = ——__ [sech (-)| ? : , (g>0) 
2J/o@ —sgn(ky )v w— kz 
ivo + k: 
(5.33) 


where w is the width of the kink as defined in Eq. (1.21). This is our final expression 
for the zero mode on the Z> kink. 

In the asymptotic vacuum, where ¢ is constant, the Dirac equation derived from 
Eq. (5.14) yields four solutions all with the same momentum. These four states are 
referred to as spin up and down states for the particle and hole (or antiparticle). On 
the domain wall, however, there are only two zero mode solutions for fixed value 
of the momentum (k,, k,). One of these has positive energy (w) and the other has 
negative energy. Therefore the two states may be called particle and hole states but 
the spin degree of freedom is not present. Consider the special case when k, = 0 
and k; Æ 0. Then we have œ = +k, and the spinor is proportional to (1, 0, 0, i)” if 
w = +k,, and to (0, i, —1, O7 if% = —k,. If we also take k, = 0, both these two 
states have w = 0 and become degenerate in energy. 

The two-fold degeneracy of the zero mode (w = 0) occurs since we are work- 
ing in three spatial dimensions where the Dirac spinors have four components. 
If we find the zero modes in one spatial dimension, the fermions are described 
by two-component spinors, and then there is only a single zero mode. If we use 
four-component spinors in one spatial dimension, it amounts to having two two- 
component spinors that do not interact with each other. Hence the degrees of free- 
dom are doubled. 

Note that the boundary conditions in Eq. (5.22) can only be satisfied if x changes 
sign in going from x = —oo to +00. So the topological nature of the kink is essential 
to the existence of the zero mode. 
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@ (0) 


Figure 5.1 The dispersion curve for fermions in the vacuum is shown on the left 
and for fermion zero modes on the domain wall on the right. 


In constructing the zero mode, we have postulated that € be an eigenstate of 
iy*. Therefore there is a possibility that there might be other zero mode solutions. 
However, it is possible to prove that this is not the case and the zero mode(s) that 
we have found are the only such solutions. The proof proceeds by choosing a set 
of orthogonal basis spinors 


1 0 i 0 
0 1 0 i 

x= ]o]> =], [> =] ]> x= 1, (5.34) 
i 0 1 0 


The first two spinors are eigenstates of iy* with eigenvalue — 1. These are the spinors 
that occur in the general solution we have already found subject to the condition 
that iy“ = —&. Since the Dirac equation is linear, any new solution must be a 
linear combination of x3 and x4. However, both these basis spinors are eigenstates 
of iy* with eigenvalue +1 and we have seen that such states cannot be part of the 
solution since the boundary conditions cannot be met. Therefore there are no other 
zero mode solutions beside the ones that we have already constructed. 

As mentioned in the introduction to this chapter, the interpretation of fermionic 
zero modes is quite distinct from that of bosonic condensates. Fermionic modes 
should be viewed as states in which the fermions can reside. A mode by itself 
does not carry energy density or charge or some other physical quantity. Only if 
the mode is occupied, can it contribute to the energy. However, the zero mode is 
special in some ways since, even if it is occupied, the fermion occupying the zero 
mode contributes zero energy. Likewise, if the zero mode is unoccupied, it also 
contributes zero energy and so the system has a degenerate ground state. Indeed 
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the occurrence of a zero mode leads to some novel and important quantum field 
theoretic consequences that we shall outline in Section 5.3. 

In the discussion of fermion zero modes above we have considered only a Yukawa 
interaction between the fermion field and the field that makes up the domain wall. 
More generally, there can also be Majorana interactions. Zero modes of Majorana 
fermions on domain walls have been discussed in [147]. 

Just like scalar field condensates and fermionic zero modes on domain walls, 
there can also be gauge field (or spin-1) condensates. These arise when the model 
has broken gauge symmetries in addition to broken discrete symmetries. This is 
precisely the situation in the SU(5) x Z2 model discussed in Chapter 2 and the 
kinks in the model have condensates of spin-1 fields as we describe in Section 5.5. 

Finally we close this section by remarking that there are several mathematical 
“index theorems” that can be used to obtain information on the number of zero 
modes on a soliton [176]. In the case of domain walls that we have been discussing, 
however, the index theorems do not lead to a useful result. 


5.3 Fractional quantum numbers 


To quantize a fermionic field we find all the modes (solutions of the Dirac equation) 
and then associate creation and annhilation operators with each of the modes.” The 
same procedure may be followed in the presence of zero modes [83]. As discussed 
in the previous section, there is a single zero mode on the Zz kink (in one spatial 
dimension), which is denoted by wo. Then the expansion of the field operator in 
modes is 


y =aoWo+ Y [bY + ahve] (5.35) 
P 


The second term is the usual sum over the positive energy modes, y,+, and fermion- 
number conjugates of the negative energy modes, ye? There is no sum over spin 
because there is no spin degree of freedom in one spatial dimension. The first term 
in Eq. (5.35) contains the zero mode, Wo, and ag is the operator associated with 
this mode. The term may seem strange because the zero mode does not have a 
corresponding conjugated term. This is because Y = Wo and so the mode func- 
tions associated with ap and its conjugate operator are identical. However, one still 
has the usual equal time anticommutation relations for the field and its canonical 
momentum 


{Walx), VIOI = 8O — Y)bab (5.36) 


2 We work in one spatial dimension in this section and hence spinors have two components. 


3 That is, Yp- is the wavefunction of a hole in the Dirac sea that has momentum p. 
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while other anticommutators vanish. Using the expansion in terms of creation and 
annihilation operators this gives 


{ao, ao} = {aj al} =0, — {ao, aj} = 1 (5.37) 


Since the Dirac Lagrangian in Eq. (5.14) is invariant if the fermion fields are 
rotated by a phase, the model has a conserved fermion number current. The Noether 
current is given by Yy” y (u = 0, 1). In the quantum theory the physical current 
operator needs to be normal ordered. This is equivalent to defining the fermion 
number operator as 


1 
Qr = fa gee 5 f exile — Wa) (5.38) 


We can act by this operator on any state to determine the fermion number of that 
state. Let us denote the state with no positive energy particles and empty zero mode 
by |0;—) and the state with no positive energy particles and filled zero mode by 
10; +). Then the fermion numbers of these two states are 


1 
Q;|0;+) = zlaga — aoa} ]|0; +) 


1 
= 5[2ayao — 1110; +) 


1 
= £510; =) (5.39) 


Therefore the kink carries a half-integer fermion number of either sign. If the 
fermion carries electric charge, the electric charge on the kink is also half-integral. 

It is critical to not think of the kink as being “kink plus fermion.” Instead the kink 
is made of both the bosonic and fermionic fields. Then there are only two states for 
the kink: one with filled zero mode and the second with the zero mode empty. 

Surprising as the half-integer fermion number is, further work in [150, 68] ob- 
tained different fractional charges in other systems (see Section 9.1). Indeed, [68] 
shows that the charges can even be irrational. 


5.4 Other consequences 


If the bosonic condensate is electrically charged, the domain wall becomes super- 
conducting. To see this in some more detail, consider the case of a complex, elec- 
trically charged, scalar field, x, interacting with the field ¢ that forms a domain 
wall 


2 


1 m À o 
L = Li] + LIA] + 51D;x? — A T= (5.40) 
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The first term is the Lagrangian for the Z) model and the second is the usual 
Maxwell Lagrangian for the gauge field A,,. The covariant derivative is defined by 


D, = 0, — igAy (5.41) 
The propagating modes of the condensate are 
X = xolg)» (5.42) 


where x(x) is the condensate profile and k is the wave-vector restricted to lie in 
the plane of the wall, the yz-plane. Since x carries electric charge, q, the electric 
current is 


om, iq t t 
j= 5 (x Vx -xVx') (5.43) 
Inserting Eq. (5.42) into (5.43), we find that the current is along the k direction 
jx = 4lxol"k (5.44) 


The simplest way to see that the wall with the condensate is superconducting is 
to write 


X = xol(xe” (5.45) 


where xo is the condensate solution and 0 is the phase variable. Then the low 
energy Lagrangian for the 0 degree of freedom can be derived by integrating the 
full Lagrangian density, Eq. (5.40), over x to get 


L(0) = (00 —eA,) (5.46) 


where we have omitted an overall constant factor. This effective Lagrangian is 
the relativistic generalization of the Lagrangian in the Ginzburg-Landau theory of 
superconductivity. Assuming that the relativistic generalization does not make any 
qualitative difference, results from the Ginzburg-Landau theory can then be applied 
directly to the present case. In particular, the domain wall with charged condensate 
can be expected to carry persistent electric currents, have magnetic vortices, and 
exhibit the Meissner effect (expulsion of magnetic fields) [61]. 

We now discuss fermionic superconductivity on domain walls. The relevant 
modes are given in Eq. (5.33) and the (normal ordered) current is 


jy =4: Yyy: (5.47) 
Using the expansion of wW in terms of creation and annihilation operators 
(Eq. (5.35)), the current in any Fock state of fermions can be evaluated. Similarly, 
the electric charge on a domain wall can also be evaluated. 


Fermions on domain walls can only make the wall superconducting if they form 
Cooper pairs and condense. It is believed that the slightest attractive interaction 
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between the fermions on the wall will lead to condensation below some critical 
temperature. On a domain wall, there are possible channels for attractive inter- 
actions. For example, the fermions interact with each other via exchange of @ 
quanta and this can lead to an attractive force. The problem of rigorously showing 
fermionic superconductivity of domain walls has not been investigated. In partic- 
ular, the Meissner effect, which is the hallmark of superconductivity, has not been 
shown. Indeed, the response of non-interacting fermion zero modes to an external 
magnetic field has been discussed with the conclusion that the walls are diamagnetic 
[173] (also see [82, 172]). 

In the particle-physics/cosmology literature, the existence of electrically charged 
zero modes is simply assumed to imply superconductivity (though see [15]). A 
reason for this assumption is that a current on a wall persists even without the 
application of an external electric field. Once the current carrying fermionic zero 
mode states have been populated there are very few processes by which these states 
can be emptied [184]. Two such dissipative processes are the scattering of counter- 
propagating fermion zero modes, and the scattering of fermion zero modes with 
fluctuations of the domain wall field itself. Generally these processes occur at a very 
slow rate, at least in astrophysical situations of interest. Hence, strictly speaking, 
domain walls in particle physics/cosmology have only been shown to be excellent 
conductors and not superconductors. 

The equilibrium current on a domain wall in any setting depends on the balance 
of the rates of current increase owing to an external electric field and decrease 
owing to dissipation. Note that an external magnetic field in which a domain wall is 
moving is, effectively, an electric field in the rest frame of the wall. Since magnetic 
fields are ubiquitous in astrophysics, any cosmological domain walls with fermion 
zero modes can be expected to be current carrying. Superconducting domain walls 
in realistic grand unification models have been discussed in [98]. 

The fermion zero mode states that we have discussed above are single particle 
eigenstates. The true states of the domain wall are also affected by fermion-fermion 
interactions. The many-body problem falls in the class of two-dimensional systems 
of interacting fermions. In the presence of a strong external magnetic field, so that 
the Landau level spacing is large compared to other energy scales, the fermions on 
the wall are similar to electrons in a quantum Hall system. 


5.5 Condensates on SU(5) x Z2 kinks 


In Chapter 2 we have discussed kinks in an SU (5) x Z2 model, which is the simplest 
example of a Grand Unified Theory. Even though SU (5) grand unification is known 
not to be phenomenologically viable, the model is still pedagogically useful. 
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The Lagrangian for the model is 
L=Lyl®, x, Xu] + Lelx, Y, Xu] (5.48) 


where the SU(5) adjoint field, ©, does not couple directly to the fermionic fields 
(denoted generically by y). Only an additional SU (5) fundamental field, x, couples 
to the fermions. The vacuum expectation value of x is responsible for electroweak 
symmetry breaking and the masses of all the observed quarks and leptons arise 
from this symmetry breaking. The SU(5) symmetry breaking has no consequences 
for the fermionic sector, except via the x field. This indirect effect can be important 
in the presence of kinks, since x can form a condensate on the kink, which can then 
interact with some of the fermions. We will discuss this further below. 
The bosonic part of the Lagrangian is 


1 
Ly = Tr(D,, ©)” + |Duxl? — V(®, x) — STX wX") (5.49) 


The covariant derivative is defined by D, = ð, — ig X „ and its action on the scalar 
fields is 


D, ,® = ð ® — ig[X,, ®], Dux = (0, — igX,)x (5.50) 


The potential is given by 


V(®, x) = V(b) + V(X) +A (Tr) x1 x + Asx) (5.51) 

with 
V(®) = =m (Tr?) + A, (Trb?)? + A2(Tr*) (5.52) 
V(x) = —m3x'x + A(x! x)? (5.53) 


Successful grand unification requires that the global minimum of the potential 
leaves an SU(3) x U(1) symmetry unbroken.‘ As already described in Section 2.2, 
the minimum of the potential with x set equal to zero, occurs at 


n 
2715 


(up to SU(5) x Z2 transformations) in the parameter range 


o = 


diag(2, 2,2, —3, —3) (5.54) 


7 
A >= 0, Aaah a0 (5.55) 


The vacuum expectation value of ® breaks SU(5) x Z2 to [SU(3) x SU(2) x 
U(1)]/Ze. If we assume that the back-reaction of a vacuum expectation value of x 


4 Symmetry breaking patterns have been discussed quite generally in [99]. 
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on that of ® is small, we can write down a reduced potential for x alone 


As 
Vied(X3 Po) = (-m3 + AgTr@) + sy) Xi Xa 


3 
+(-m3 + AgTr®? + ser) xixe + rs(xix (5.56) 


where a = 1, 2,3 and b = 4,5. The symmetry is broken to [SU(3) x U(1)]/Z3 
only if the vacuum expectation value of x is along the x4 or xs directions. This 
further restricts the range of parameters to 

Sy > m5 — A4Tr@z > 2s p, A3 > 0 (5.57) 
which also implies As < 0 and m? < AgTr®5. We assume that these conditions on 
the parameters are satisfied. Then a minimum of the reduced potential occurs at 


x! = New(O, 0, 0, 1, 0) (5.58) 
where 
1 3A5 
2 2 2 2 
= ——( m — Astro} — 5.59 
New 213 (mi 4 r 0 20 n ) ( ) 


is the electroweak symmetry breaking scale. The final [SU(3) x U(1)]/Z3 sym- 
metry corresponds to the color and electromagnetic symmetries present today. 

Next we describe the fermionic sector.’ There are two fermion fields: y, which 
is in the fundamental (5-dimensional) representation of SU(5) and ¢, which is in 
the antisymmetric 10-dimensional representation. The known quarks and leptons 
fit within the components of these fields. With the choice of vacuum expectation 
values in Eq. (5.58) 


(WL = (d, d?, d°, e7, —ve)L (5.60) 

(Wr = (di, do, dz, e", —vS), (5.61) 
0 ues —y2 uy di 
i —u®? 0 uc! uz də 

Gi =| u? -u 0 m d (5.62) 
v2 —u; —uü? —u3 0 et 


—d 1 —d> —d3 —et 0 L 
(see Eq. (14.9) in [30]). The numerical index on the u and d fields refers to color 
charge, and the placement (subscript or superscript) depends on the representation 


5 Actually we describe only one of the three families of the standard model fermionic sector, and then too the 
neutrino is taken to be massless. 
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(fundamental or fundamental conjugate) in which the field transforms under the 
unbroken SU(3). The c superscript denotes charge conjugation: 


wo = iy y" (5.63) 
The L and R subscripts refer to left- and right-handed spinors 
pa 1+; 
ar ea ea N (5.64) 


The Dirac y matrices are defined in Eqs. (5.15) and (5.17). 
Now we are ready to describe the interactions of the various fields with the 
SU(5) x Zz kink, described as the q = 2 kink in Chapter 2. 


e In the presence of a (q = 2) kink, the vacuum expectation values are 
n 


2/15 
x" (—00) = new(0, 0, 0, 1, 0) 


Nn a 
(+œ) = — diag(2, —3, —3, 2, 2 
(+00) 2/15 g( ) 


X7 (+00) = new(0, 0, 1, 0, 0) 


P(—00) = + 


diag(2, 2,2, —3, —3) 


Note that the non-trivial entry of x has to coincide with one of the —3 entries of 
® since this is what minimizes the potential V(®, x). Therefore x must rotate 
through the kink. Inside the kink, the fields are not pure rotations of the asymptotic 
values. 

The component ıı goes from +2 to —2 as the wall is crossed. Hence it must 
vanish in the wall. This is very similar to the Z2 case, where the field vanishes 
at the center of the wall. The field x interacts with ® as given by the potential 
in Eq. (5.51). Note the interaction term 5Tr( x'@?x), which directly couples x; 
to ®,,. (The other term couples all components of x to Tr only.) By explicit 
construction it can be seen that x; can condense on the wall for a certain range 
of parameters [146]. Hence the SU(5) x Z2 model allows for scalar condensates 
on the wall (see Section 5.1). In addition, since x; is a complex scalar field, the 
condensate has an associated phase. The choice of phase on different parts of the 
wall may be different, leading to vortices in x; that can only exist on the wall. 
Since x; transforms non-trivially under the unbroken SU (3), the vortices carry 
color magnetic field. This is similar to our discussion below Eq. (5.46). 

Next we consider fermion interactions with the wall [146]. The fermions do not 
couple directly to ®. Hence the only coupling to the wall is due to the rotation of x 
in passage through the wall and to the condensate in the x; component. Consider 
the scattering of the fifth component, Ys, which corresponds to a neutrino on the 
left side of the wall but a d-quark on the right. This fifth component has non-zero 
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reflection and transmission probability. If it reflects, the particle is still a neutrino. 
If it transmits, it must change into a d-quark. If a neutrino becomes a d-quark in 
passing through the wall, it must pick up electric and color charge from the wall. 
Hence we are forced to conclude that there must be electric and color excitations 
that live entirely on the wall. If a x; condensate is not present, the only available 
excitations are the charged gauge field components. Hence charged gauge fields 
must condense on the wall. 

To see the presence of a charged gauge field condensate, it is most convenient 
to go to a gauge where the scalar field vacuum expectation values are oriented in 
the same directions on both sides of the wall, as we now discuss. 

Consider a very thin wall, so that 


n ; 
(x <0)=+ diag(2, 2,2, —3, —3) = ® 
( ) NG g( ) 0 
OCs 0S T 3 = 99,9) (5.65) 


2/15 


Now we perform a local gauge transformation that rotates ® into the direction of 
Po (up to a sign) everywhere. Such a gauge rotation is local since it is equal to the 
identity transformation for x < 0 but is non-trivial for x > 0 since it exchanges 
the 23- and 45-blocks of ®. In both regions, x < 0 and x > 0, the gauge rotation 
is constant. The rotation is non-constant only at x = 0 i.e. on the wall. Hence the 
gauge fields after the rotation vanish everywhere except on the wall itself and there 
are gauge degrees of freedom residing on the wall. A more explicit calculation 
shows that the gauge fields living on the wall carry electric and color charge. 


5.6 Possibility of fermion bound states 


In addition to fermionic zero modes on a kink, there may also be fermionic bound 
states. Such bound states would have a non-vanishing energy eigenvalue w with 
0 < œ < m. Since the energy eigenvalue is less than the asymptotic mass, the 
fermion would be bound to the wall. We examine whether the model in Eq. (5.14) 
leads to fermionic bound states on a Z> kink. 

For convenience we work in one spatial dimension. Then spinors have two 
components and there are only two gamma matrices, which can be taken to be 


03 {1 0 bow (f0 A 
y =o =(j a ysi SNo (5.66) 


Then the Dirac equation id — gd, Ww = 0 together with y = exp(—iwrt)é gives 


ðs = —(@ + Bdx)&2 
0x62 = +(@ — gdx)éi (5.67) 
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where 


_ ( &1) 
D= i) (5.68) 


and we are interested in solutions with 
0<w<m = gn (5.69) 


The boundary conditions at the origin for €; and & may be determined by 
noting that we are free to rescale both &; and & by a constant factor. So we can 
set €,(0) = +1. Further, using the symmetry (x) = —¢(—x), we find that the 
equations are invariant if we replace £1 (x) by c&.(—x) and &(x) by c&;(—x), where 
c is aconstant. Hence 


&1(x) = c&2(—x), §2(x) = c§1(—x) (5.70) 
This gives 
E) = c&(—x) = c°& (x) (5.71) 
Since (x) cannot vanish for all x, we get 
c= +1 (5.72) 
Therefore there are two possible boundary conditions at the origin 
§2(0) = +é, (0) = +1 (5.73) 


At infinity we require £; —> 0 and & — 0. 

A numerical search for a solution with non-zero œw did not reveal any bound 
states for the range of parameters 0.1 < mew < 20, where w is the width of the 
kink. However this does not exclude the existence of fermion bound states (beside 
the zero mode) on kinks in other systems, and it remains an open problem to find 
systems where such bound states exist. 


5.7 Open questions 


1. Explore the classical and quantum physics of a domain wall with electrically charged 
bosonic and fermionic zero modes placed in an external magnetic field. What happens 
if the domain wall is moving? 

2. Calculate the reflection of photons off a superconducting domain wall. Is the wall a good 
mirror? (See [184].) 

3. Construct a system in which the kink has both a zero mode and a fermionic bound state. 
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Formation of kinks 


In this chapter we study the formation of kinks and domain walls during a phase 
transition. We start by describing the effective potential for a field theory at finite 
temperature. This sets up a useful framework for discussing phase transitions and 
defect formation. 


6.1 Effective potential 


The effective potential is a tool that is often used to study phase transitions in 
field theory [89, 179, 90, 47, 100]. The idea is to consider the interaction of 
a scalar degree of freedom (“order parameter”) with a thermal background of 
particles. Such processes induce additional temperature dependent terms in the 
potential for the order parameter, leading to an “effective potential.’ The shape 
of the effective potential varies as a function of temperature and new minima 
might appear. The global minimum defines the vacuum of the model. If a new 
global minimum appears at some temperature, it indicates that the system makes 
a transition to a new expectation value of the order parameter and there is a 
phase change. We now describe the (one loop) effective potential in a little more 
detail. 
We consider a field theory of scalar, spinor and vector fields 


L = Lg + Lp (6.1) 
with the bosonic Lagrangian 
Lg = 5(Dy®)D"®, — V(®) — Fi, Fe (6.2) 
where ®; are the components of the scalar fields, 
Dy = ðu — ie A; T? (6.3) 


90 
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the T° are group generators, and 
Fa, = ð AS — 3p A! + ef” A’ AS (6.4) 


where A‘) are the gauge fields. 
The Lagrangian for a fermionic multiplet W is 


Lp = i¥y“D,Y a WT; UO; (6.5) 


where T; are the Yukawa coupling matrices. The quantity W denotes a collection of 
fermionic fields and the Yukawa coupling term may be written more explicitly as 
we id A Yi ®; where a, £ label the various fermionic fields, the superscripts o, p 
on the fermion fields are spinor indices, and i labels the interaction term with the 
scalar field ®;. T; has spinor indices because it could contain the unit matrix (vector 
coupling) and/or the y> matrix (axial coupling) defined in Eq. (5.17). 

If the expectation values of the scalar fields are denoted by ®o;, then the mass 
matrices of the various fields are written as 


j 3? V 
Hij = z , scalar fields (6.6) 
i dD,00; b=) 
m =T;®o;, spinor fields (6.7) 
M?, = e°(TaTo)ij Po; Doj, vector fields (6.8) 


where a, b are gauge field group indices. 
Then the finite temperature, one loop effective potential is! 


M ao Re os 
ett(Po, T) = ® — T^ — —NT ; 
Vere(Po, T) = V (Po) + z4 T (6.9) 
where 
7 
N =M; + z^ (6.10) 


is the number of bosonic and fermionic spin states, and 
1 
M? = Tr(u® + 3Tr(M?) + 5T"°my’m) (6.11) 


where y° is defined in Eq. (5.15). Note that M? depends on the expectation value 
Po through the defining equations for the mass matrices given above. For example, 
M? contains a term proportional to Tr(®%). 

An important feature of the effective potential is that it can show the presence 
of phase transitions. If there are scalar fields with negative mass squared terms in 


l Radiative corrections and spontaneous symmetry breaking are discussed in [33, 178]. 
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Figure 6.1 Sketch of effective potential for first-order phase transition (left) and 
second-order phase transition (right). In the first-order case, the global minimum of 
the potential at high temperature (o = 0 in illustration) becomes a local minimum 
at low temperature. In the second-order case, the global minimum of the potential 
at high temperature becomes a local maximum at low temperature. The effective 
potential at Po = 0 decreases with increasing temperature because of the last term 
proportional to -M T4 in Eq. (6.9). 


V(®), the contributions from the M?T? term in the effective potential, Eq. (6.9), 
can make the effective mass squared positive for these fields if the temperature 
is high enough (see Fig. 6.1). Therefore when the system is at high temperature, 
the effective mass squared can be positive and the minimum of the potential at 
o = 0. As the system is cooled, the effective mass squared becomes negative and 
the minimum of the effective potential occurs at non-zero values of ®ọ and the 
lowest energy state has shifted from o = 0 to By Æ 0. The order parameter, ®, 
acquires a non-zero “vacuum expectation value” at some critical temperature. This 
is the phenomenon of spontaneous symmetry breaking and manifests itself as a 
phase transition. The phase at high temperature had a certain symmetry dictated by 
the invariance of the field theory with o = 0 and at low temperature the symmetry 
is changed because now ®ọo Æ 0. 

As a simple example of an effective potential, consider the Ad* model of Eq. (1.2) 
with 
À 
4 


mõ 2 À 4 4 
ae aa (6.12) 
Then u? of Eq. (6.6) is given by 
uR? = —m2 + 3a (6.13) 


and since there is only one scalar field in this model 


M? = —m? + 31H? (6.14) 
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Therefore, up to a term that is independent of ġo, the effective potential becomes 


-2 À 
mM” 9 4 
Vireo, T) = z F 4% (6.15) 
with 
-2 2, Àm 
m = -m + A (6.16) 


Note that the masses of small excitations around the true vacuum are given by 
(Po) with ġo being the vacuum expectation value. By minimizing Vet (Eq. (6.15)) 
we get 


Pomin T) =0, m >0 (6.17) 
—m2 2 
=; ——, m<0 (6.18) 
À 
leading to the mass squared for small excitations (particles) in the true vacuum 
À 
men = Vel Pomin) = F(T? - T) T > Te (6.19) 
À 
= L -T’, T<T (6.20) 


where T, is the critical temperature 


T 2mo 


T, 
c Ji. 


(6.21) 


In cosmology, since the universe is expanding, it is also cooling. Therefore 
we can have one or many cosmological phase transitions and the particle-physics 
symmetries at high temperatures (early universe) and low temperatures (recent 
universe) are different. The symmetry after the phase transition can be smaller 
or larger than the symmetry before the phase transition. In other words, lowering 
the temperature can spontaneously break or restore a symmetry. We will mostly 
consider symmetry breaking during the phase transition but examples of symmetry 
restoration are also easy to construct. A system in which symmetry restoration is 
observed is Rochelle salt [85, 179]. 


6.2 Phase dynamics 


The effective potential Veg(®o, T) is calculated for a system that is in thermal equi- 
librium, assuming a homogeneous vacuum expectation value of the order parameter 
Po. Yet thermal equilibrium is not maintained during the phase transition and also 
the phase change occurs in an inhomogeneous manner. The dynamics are clear 
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Figure 6.2 A schematic diagram of bubbles nucleating in a first-order phase 
transition. The two types of bubbles correspond to the two different values of the 
order parameter. The bubbles grow and collide, and new bubbles nucleate as well. 
Eventually the whole system is in the new phase. 


for a first-order phase transition in which the high temperature phase becomes a 
metastable state (see Fig. 6.1) at some critical temperature. Now the system can be 
stuck in this metastable state even when the temperature drops significantly below 
the critical temperature. An external perturbation can cause the system to transition 
to the global vacuum. In the absence of an external perturbation, quantum tunneling 
can trigger the transition. In either case, bubbles of a critical size of the true vacua 
(®o Æ 0) nucleate in the false vacuum (®p = 0) background (see Fig. 6.2). These 
bubbles grow and eventually merge thus filling space and completing the phase 
transition. Clearly this process is not homogeneous and cannot be described by an 
effective potential. 

In a second-order transition, in contrast to a first-order transition, there is no 
metastable state in which the system can be trapped. Thus ®ọo evolves continuously 
(“spinodal decomposition”) from o = 0 to Po Æ 0. However, different spatial 
regions evolve at different rates owing to thermal and quantum fluctuations, and ®o 
is not spatially uniform. Once again, since the effective potential assumes constant 
®o, it can indicate the existence of a second-order phase transition but cannot 
be expected to accurately describe the dynamics of the transition. Since defect 
formation crucially depends on the inhomogeneities of the order parameter during 
the transition, new ideas have been needed to predict the statistical properties of 
defects formed in a second-order phase transition. 

In one spatial dimension, the distribution of kinks is described by the number 
density of kinks, and correlators of kink locations. In higher dimensions, the prob- 
lem becomes richer because domain walls are extended and can curve and have 
complicated topology. In addition to the mass density in domain walls, we are 
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interested in the statistical distribution of shapes and sizes of domain walls formed 
at the phase transition. 


6.3 Kibble mechanism: first-order phase transition 


At a first-order phase transition, the order parameter has to change from ®p = 0 to 
its non-zero vacuum expectation value. We are interested in the case when there is 
more than one possible non-zero value for ®y. Then the dynamics in a small spatial 
region select a vacuum. However the vacuum selected in different spatial regions 
can be different. For example, in the case of the Z2 model, the field in a certain 
region might relax into the @ = +7 vacuum, whereas in another region it might 
relax into ġ = —n (see Fig. 6.2). 

In a first-order phase transition, each bubble is filled with constant ®o i.e. a 
fixed vacuum is chosen within a bubble but it can be different for different bubbles. 
With time, the bubbles grow and collide and fill up the volume. Let us denote by 
€ the characteristic size of a region where the same vacuum is selected, after the 
phase transition is over. Then & is the typical size of bubbles when they percolate. 
If I denotes the bubble nucleation rate per unit volume and v is the velocity of 
the growing bubble walls, then we can define a length scale and a time scale on 
dimensional grounds (in D spatial dimensions) by 


- pvi) E j y AREN 
E= (=) t= (sr) (6.22) 


The domain size £ is a numerical factor times Ẹ and in practice we take £ ~ Ẹ. 
Similarly 7 is related to the time that it takes to complete the phase transition.” 

The process of bubble percolation has been studied both analytically and nu- 
merically [95]. Taking the centers of the bubbles as the vertices of a lattice and 
connecting only the centers of bubbles that collide, we obtain a random lattice (see 
Fig. 6.3). We would like to determine the characteristics of such a random lattice 
since this plays a role in determining the network of defects that form. For example, 
the typical number of bubbles with which any given bubble collides, also known 
as the “coordination number” of the random lattice, plays a role in the fraction of 
closed topological defects (closed domain walls, loops of string, or closely paired 
monopole-antimonopole pair) that are formed. 

In one spatial dimension, every bubble trivially collides with two other bubbles. 
In two spatial dimensions the average number of collisions is the same as the 
coordination number of a fully triangulated lattice of infinite extent. From purely 


2 Equation (6.22) relates £ to the nucleation rate I’, but it is very hard to measure I in any experiment. In fact, it 
may be easier to measure properties of the defect network, then £, and from it infer T. 
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Figure 6.3 If two bubbles collide, their centers are joined by straight lines. The 
figure then shows the “random bubble lattice” expected in a first-order phase 
transition in two spatial dimensions. 


geometrical constraints that we describe next, the coordination number is six (see, 
for example, [129]). 

The lattice is infinite in extent and by identifying the points at infinity we can 
view the lattice as lying on a two-dimensional sphere. Then Euler’s formula relates 
the number of vertices (V), edges (E£) and faces (F) of the lattice 


V—-E+F=2 (6.23) 


Let the coordination number be n. Therefore for every vertex there are n edges 
but every edge is bounded by two vertices. This relates the number of edges to the 
number of vertices: E = nV/2. Also, every face is a triangle, giving three edges 
to every face. But an edge is shared by two faces. So E = 3F'/2. Putting together 
these relations in Euler’s formula gives 


V —nV/24+nV/3=2 (6.24) 


In the limit of V —> ov, this yields n = 6. 

In three spatial dimensions, similar arguments have been given [95] to show 
that the average coordination number is 13.4. This result is not completely fixed 
by geometrical constraints as in two dimensions and the result can vary a little 
depending on the details of the bubble size distribution. 

Returning to the 4¢* model, each bubble either has the phase pọ = +7 or œo = 
—n within it. If bubbles of different phases collide, a domain wall forms between the 
centers of those bubbles. If bubbles of the same phase collide, a wall does not form, 
though it is possible that a closed domain wall or a wall-antiwall pair forms owing 
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to the energetics of the bubble collision. We expect small closed walls and closely 
separated walls and antiwalls to annihilate. Hence the distribution of domain walls 
after the phase transition is simply described by the locations of bubble collisions 
when the bubbles carry different phases. Since the phase in the bubbles is +n with 
equal probability, the phase transition is simulated by assigning +7 to each of the 
vertices of the random bubble lattice as in Fig. 6.3. We shall further discuss the 
properties of the wall network at formation during a first-order phase transition in 
Section 6.6. 

As we have seen, a first-order phase transition is relatively simple to conceptual- 
ize. A second-order phase transition is harder to understand. To discuss second-order 
phase transitions, it is useful to first define an equilibrium correlation length. 


6.4 Correlation length 


The “equilibrium correlation length,” Æ, is defined as the distance over which field 
correlations are significant. Generally the field correlations at two spatial points 
fall off exponentially with increasing separation between the points, exp(—r/é), 
and the exponent defines the equilibrium correlation length, £. Hence we need to 
evaluate the correlation function 


G(r) = (Tiot, YP, yT) (6.25) 


where G only depends onr = |x — y| because the system is translationally invariant. 
The thermal state is denoted by |T} and is defined as the state containing the 
equilibrium number density distribution (Fermi-Dirac or Bose-Einstein) of particles 


IT) = \{nx})r (6.26) 


1 
(6.27) 


Me efo% +1 


where 6 = 1/T and ax is the energy of particles in the k mode.* 
With the Z model in mind, we have only one scalar field and the quantum field 
operator can be expanded in modes about the true vacuum 


dk 1 
(27)? Vor 
where ġo(T) is the vacuum expectation value of the field at temperature T, and a, 
and al are annihilation and creation operators. The dispersion relation is that for 


a free particle with temperature dependent mass m(T) (see Eq. (6.20) for the A¢* 
model; we have dropped the subscript “eff” for convenience) 


[emia tilexgy, 4 etiontikex gt] (6.28) 


b(t, x) = P(T) + J 


o =k? +m? (6.29) 


3 The chemical potential vanishes in the present case. 
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The thermal state |T) contains a Bose-Einstein distribution of the particle exci- 
tations and the number of scalar particles at momentum k is given by 
1 


= obo — | (6.30) 


Nk 
where 8 = 1/T (Boltzmann constant has been set to 1). 
By inserting the expansion in Eq. (6.28) in the correlator, we find 
dk 1 dp 1 
CrP Vor) CY /2op 


where K is a constant which is independent of temperature and proportional to 
6°)(x — y). Then 


G(r) = elfkx-PY)(TIalay|T) + K (6.31) 


wk e ik (x—y) 


a 
(27)? ef% — 1 ty 


G(r) = 


T 
A A (6.32) 
Air 


where, in doing the integral, we have assumed m(T) «< T. From the final expression 
we get the equilibrium correlation length 


1 


E = —__ 6.33 
3 mT) (6.33) 
For the Z2 model (Eq. (6.20)) 
m? = T —T2), T>T (6.34) 
À 
= e -T?°), T<T, (6.35) 


Therefore the equilibrium correlation length is 


T >T, (6.36) 


7 4 1 
T) = ,/-——= 
age : T <T. (6.37) 
— i T2- T? < le . 
or 


E(T) x |T — T7"? (6.38) 


The essential feature in £ is the singularity at T = T,, that occurs at a time that 
we denote te. Assuming that the cooling (quench) occurs at a constant rate Te /Text 
(in a range of temperature around T,), we have 


T, 
T —T, = —— (t — t.) (6.39) 
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Figure 6.4 Sketch of equilibrium correlation length as a function of time as given 
in Eq. (6.40). 


hence we write 
E~n |T = T4” « |t t|” (6.40) 


for T close to T.. The exponent v is called a “critical exponent” and the mean 
field theory calculation described above gives v = 1/2. However, the mean field 
theory ignores particle interactions and renormalization group methods give v = 
2/3, which is closer to experiment. A sketch of the shape of & (Eq. (6.40)) is shown 
in Fig. 6.4. 

One shortcoming of the mean field calculation of E is that we have quantized 
the field ¢ in a fixed true vacuum so that ¢@0(T) in Eq. (6.28) is independent of x. 
This assumes that the same vacuum is chosen everywhere below T. On the other 
hand, we are precisely interested in the spatial extent of a region ina single vacuum. 
Hence a more suitable expansion of ¢@ would be 

aa : i] (6.41 
Ons Fag l it x)ak + fet, ay] (6.41) 
instead of Eq. (6.28). The vacuum expectation value, ġo, is now allowed to depend 
on both ¢ and x since the background domain walls may be non-static. The second 
term in the expansion describes small fluctuations (particles) with mode functions 
J, in the classical background @o(t, x, T) at x. 

The expansion in Eq. (6.41) is not as obvious as might appear at first sight. We 
have seen in Chapter 4 that a kink can itself be written in terms of particles via the 
Mandelstam operator. This was done in the sine-Gordon model but it is conceivable 
that such an operator also exists in other models. So the above expansion only makes 


ot, x) = polt, X, T) + 


sense for a state in which there is a clear separation between the particle degrees 
of freedom appearing in the sum and the soliton degrees of freedom included in 
Qo. For example, if the walls are very close to each other the separation of the two 
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terms may not be justified. Hence the phenomenon of defect formation is closely 
tied to the separation of classical (soliton) and quantum (particle) variables. 
We are interested in 


éo = (T|do(t, x, T)go(t, y, DIT) (6.42) 


However, we have no way of calculating this equal time “domain correlation func- 
tion” since (i) the thermal state refers to a thermal distribution of particles, not of 
domains, and (ii) the defects do not remain in thermal equilibrium with the par- 
ticles. This impasse is made less severe by realizing that the calculation of & for 
T > T, does not suffer from this problem since then ġo = 0 is the unique vacuum. 
We expect the correlations for T < T, to be determined by those for T > T, and so 
it might be sufficient to know the correlation length for T > Te. (We discuss this 
further in Section 6.5.) To emphasize this point, the two branches of the sketch in 
Fig. 6.4 correspond to two different quantities — the equilibrium correlation length 
for t < t is for excitations in a different vacuum from that for t > te. So, while 
Eq. (6.40) describes the correlations of particle excitations in a given true vacuum 
for T < T, (t > te), we cannot expect that this has anything to do with the size of 
domains of constant vacuum. 

The next subtlety is that the divergence of the equilibrium correlation length at 
T = T, should not be taken too literally. The reason is that there is an external 
agency (refrigerator) driving the phase transition on a time scale given by Text. As 
the system gets closer to the critical temperature, it takes longer for equilibrium 
to be established, while the external agency continues to cool the system at a rate 
determined by external factors. At some temperature above the critical temperature 
the time taken to maintain equilibrium becomes larger than the time scale at which 
the external conditions are changing. 

Assume that the external temperature is being lowered at a constant rate 


T=T, (1 = =) (6.43) 


Text 
where T, is the critical temperature and we have chosen T, = 0 for convenience. 
The equilibrium correlation length has the form 


ECT) = nle” (6.44) 
where v is a critical exponent, 7 is some unspecified length scale, and 
e=1- 2 (6.45) 
T, 


(6.46) 
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This equation shows that as € > 0—, Ẹ must change at an ever faster rate if equi- 
librium is to be maintained. 

The relaxation rate can be obtained by perturbing the system and finding 
how long it takes for the perturbation (“sound”) to equilibrate. The result is the 
“relaxation time” 


Trel = tole|-" (6.47) 


where u is another critical exponent. Then the “speed of relaxation” is the sound 
speed 


$. = I jeje 
Trel To 


c(T) = (6.48) 


Note that Tre, diverges as T approaches T}. This is called “critical slowing down.” 
When the system cannot keep up with the external changes, equilibrium is lost. 
Denoting the temperature at which Tre, becomes equal to Text by T; we find 


VT 1/(u+1) 
y en el ik ( ) (6.49) 
Text 


VTo 1/(u+1) 
tx = ~tx( 2) (6.50) 


Text 


which occurs at 


So we expect the correlation length £ for T > T, to be equal to the equilibrium 
correlation length £ until time tą, after which € departs from € and grows more 
slowly (see Fig. 6.5). The behavior of Ẹ between t, and tẹ is not known and it 
is generally assumed that £ does not change very much in this interval. After ts, 
there are two distinct vacua, and we need to consider both the correlation scale of 
chosen vacua (denoted by £o) and the correlations of excitations within a chosen 
vacuum, £. As time goes by, walls annihilate and the domain size with a given 
vacuum grows. We discuss & in the next section. 


6.5 Kibble-Zurek mechanism: second-order phase transition 


The domain correlation length, &), over which the same vacuum is chosen, is dif- 
ferent from the equilibrium correlation length denoted by & (Eq. (6.40)). It is also 
different from the correlation length € obtained for particle excitations, including 
the phenomenon of critical slow down, since &) has nothing to do with particle 
excitations. We now discuss different approaches to estimating £o (for a review see 


[10}). 
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Figure 6.5 The correlation length at high temperature (t < t) increases as the 
critical temperature is approached, departing from the equilibrium correlation 
length when critical slowing down becomes important at tą. Below the critical 
temperature, there are two correlation length scales of interest. The domain cor- 
relation length, £o, relates to the extent of the spatial domains that are in the same 
vacuum. This is precisely the spacing of domain walls. Wall-antiwall annihilations 
cause & to grow with time. The particle correlation length, £, however, decreases 
with time since the mass of the particles grows, and eventually approaches the zero 
temperature value. The dynamics of how & separates from & are not understood 
and are denoted by the shaded region. 


To estimate & in the case of a second-order phase transition, Kibble [87] used 
two different criteria. First, he obtained an upper bound to £o in the cosmological 
context based on causality considerations. If the phase transition takes place at a 
certain cosmic time T, then the vacua at points separated by more than ct, where c 
is the speed of light, must have been selected independently since ct is the size of 
the cosmic horizon. Hence £ < ct. This is the “causality bound.” 

The second estimate is based on finding the Ginzburg length. This is the length 
over which the choice of vacuum cannot change owing to thermal fluctuations. For 
concreteness, let us imagine that there is a domain of size / in which ¢ = +¢@9(T) 
(T < T.) in the Zz model. In one spatial dimension this corresponds to a wall- 
antiwall separated by a distance /, and in three dimensions it corresponds to a 
closed domain wall of characteristic size l. The idea is that, if l is small, thermal 
fluctuations can spontaneously change the phase within the domain from ¢ = +¢9 
to $ = —¢o. However, if / is large, the phase in the domain is frozen, and the 
distribution of defects does not change spontaneously owing to thermal fluctuations. 
The smallest length / for which a domain is frozen defines the distance between 
closest defects and hence predicts the number density of defects. 

The energy required to change the phase in a volume R? is given by R?AV(T) 
where AV is the free energy density difference between the minimum and maximum 
of the potential at a temperature 7. The thermal fluctuation energy available per 
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excitation mode is T according to equipartition. Equating the required and the 
available energies gives 


R'AV(T) ZT (6.51) 


Therefore, at temperature T < Te, a region that is smaller than 


T 1/3 


will have enough thermal energy to fluctuate from one vacuum to the other. For 
example, in the Z model (see Eq. (6.15)) 


—4 


m 


Therefore, at temperature T < T,, the length scale below which regions are still 
fluctuating are 


4 T 3 yB T7723 
R(T) = 7173 (Ges ry = ATT, f E z (6.54) 


where the last approximation holds for T ~ Ty. 

For a region to fluctuate from one vacuum to another, not only does it need 
the energy to jump over the barrier, but all different parts of the region need to 
jump together. This means that all the particles in the domain should be activated 
coherently. The particle coherence scale is described by the correlation length, 
which is approximated by the equilibrium correlation length, £. Therefore, at a 
temperature 7’, regions of size less than /; = min(R(T), &(T)) (subscript “f” stands 
for “fluctuating”) can actively change vacua. The Ginzburg temperature, Tg, is 
defined by the condition R(Tg) = &(Tc), and the Ginzburg length is defined by 


lg = &(Tg). For the Z} model, this gives 


A (T? = 72)? (6.53) 


c 


T, — To ~ ÀT, (6.55) 


lo = (To) © 


6.56 
IT, (6.56) 


Early estimates took the Ginzburg temperature to be the epoch when domain walls 
are formed. The number density of walls then follows by dimensional analysis as 
~ 1/12. 

The relevance of the Ginzburg temperature for defect formation is not clear. As 
discussed in the previous section, the correlation length & is calculated for particle 
excitations in a given vacuum, whereas we are interested in the correlation length 
of the vacuum domains denoted by & . In fact, experiments in He-3 find that defects 
are produced at a temperature below T, but above Tg, implying that the Ginzburg 
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criterion is not a necessary condition for defect formation. A discussion of the 
relevance of the Ginzburg criterion in the context of vortex formation in He-3 and 
He-4 may be found in [86]. 

Zurek estimated the domain size, £z, by considering the time scales involved 
during the phase transition [187, 188, 93]. As discussed at the end of Section 6.4, the 
system cannot keep up with external changes at t = tą (Fig. 6.5). Zurek postulated 
that the correlation length at the instant when the system can no longer keep up 
with the external changes determines the size of the domains that get frozen. This 
in turn determines the number of defects. 

To estimate & at te we know £ at the time critical slowing down becomes im- 
portant. To this we add the distance that a perturbation can propagate from the 
slow-down time, tą, to the phase transition time, f, (see Fig. 6.5). That gives us 


te 
fut) = &.) f dret (6.57) 

ty 

Lage fra e 
= n~ | — 6.58 
My +u-v (=) 28) 
The crucial part of this relation is 
Text v/(+H) 

Eo(tc) X ( ) (6.59) 

To 


This relation gives the dependence of the number density of domain walls in D 
spatial dimensions on the external time 


Text —vD/(+p) 
na ( ) (6.60) 


To 


We can control Text in experiments and hence this is a testable prediction. 

The above analysis can be improved yet further. For example, we have calculated 
éo att = te. Yet thermal fluctuations after t = te (i.e. T < Tę) may be important and 
the domain structure may freeze out at yet lower temperatures, as in the discussion 
of the Ginzburg length scale above. So the relevant time at which & is stable to 
thermal fluctuations is somewhat after to, in agreement with the analysis of [8]. 

There is yet another view of defect formation at a phase transition first proposed 
in [5]. In numerical simulations of a U(1) field theory, the authors found that 
there is a distribution of vortices even at temperatures above the phase transition. 
However, these vortices are small, closed structures. At the critical temperature, 
the vortices link up and form infinite, open structures. Thus the phase transition is 
coincident with a percolation transition of the vortices. If this feature is generally 
true, we expect a population of small, closed domain walls to exist above the critical 
temperature. As the temperature is lowered, the walls connect and grow larger and 
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Table 6.1 Size distribution of black clusters found by simulations 
on a cubic lattice. 


Cluster size 1 2 3 4 6 10 31082 


Number 462 84 14 13 1 1 1 


at the critical temperature, the walls percolate, giving walls of infinite extent. The 
percolation picture has not been checked by simulating the domain wall forming 
phase transition. However, we can still study the statistical properties of the network 
of walls formed after a phase transition using some simple arguments that we now 
describe. 

The topic of defect formation and, more generally, phase transition dynamics is 
still under active investigation. 


6.6 Domain wall network formation 


The previous sections focused on the density of domain walls that can be expected 
to form during a suitable phase transition. In this section we focus on a somewhat 
different aspect of the problem: what are the statistical properties of the domain 
walls formed at a phase transition? Are the domain walls formed as little closed 
spherical structures? Or are they infinite and planar? First we discuss the simple 
case of a network of Zz walls and then the more complicated case of SU(5) x Z2 
walls. 


6.6.1 Z network 


The properties of the network of Z2 domain walls at formation have been determined 
by numerical simulations implementing the “Kibble mechanism.” The vacuum in 
any correlated region of space is determined at random. Then, if there are only two 
degenerate vacua (call them black and white), there are spatial regions that are in the 
black phase with 50% probability and others in the white phase. The boundaries 
between these regions of different phases are the locations of the domain walls 
(see Fig. 6.6). 

Numerical simulations of the Kibble mechanism on a cubic lattice gave the 
statistics shown in Table 6.1 [74, 159]. The data show that there is essentially one 
giant connected black cluster. By symmetry there is one connected white cluster. 
In the infinite volume limit, these clusters are also infinite and their surface areas 
are infinite. Therefore the topological domain wall formed at the phase transition 
is infinite. 
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Figure 6.6 The distribution of two phases (black and white) on a square lattice 
in two spatial dimensions. Domain walls lie at the interface of the black and white 
regions. 


6.6.2 SU(5) network 


What does the Kibble mechanism predict for SU (5) x Z domain walls? Just as in 
the Z2 case, we have to throw down values of the Higgs field on a lattice, assuming 
that every point on the vacuum manifold is equally likely, and then examine the 
walls that would form at the interface. In Section 2.2 we have found that there are 
three kinds of wall solutions in this model and we have labeled the walls by the 
index q, which can take values 0, 1, or 2. Each kind of wall has the same topology 
but they have different masses. Each wall type is formed with some probability. 
Based on the Kibble mechanism, the probability that a certain wall forms is directly 
related to the number of boundary values that result in the formation of that kind of 
wall. So we need to evaluate all the boundary conditions that lead to domain walls 
with a certain value of q. 

The space of boundary conditions leading to a given type of domain wall is 
discussed in Section 2.4. However, similar considerations occur in simpler models 
and it is helpful to think of the problem in a discrete case, for example the Ss x Z2 
kinks described in Section 2.5. Take a fixed (discrete) vacuum in one domain. The 
neighboring domain can be in any other vacuum state with equal probability. 
There are ten possible states for the neighboring domain. Only one of these gives 
the q = 0 wall, six give the q = 1 wall, and three give the q = 2 wall. Then the 
Kibble mechanism implies that the network contains g = 0, 1, 2 walls and their 
number densities are in the ratio 1 : 6 : 3. This means that the network is domi- 
nantly composed of the q = 1 wall. However, since the q = 2 wall is the lightest 
wall, the g = 1 walls formed by the Kibble mechanism during the phase transi- 
tion subsequently decay into g = 2 walls. We will show some evidence for this 
two-stage process in Section 6.7. 
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Similarly we can identify the space of boundary conditions that lead to a particular 
kind of kink in the SU(5) x Z2 model. We considered this problem in Section 2.4 
and listed the spaces in Table 2.1. From the table we read off that the space of 
boundary conditions leading to the g = 0 kink is zero dimensional, for the g = 1 
kink it is six dimensional, and is four dimensional for the g = 2 kink. Since a six- 
dimensional space is infinitely bigger than a 0- or a four-dimensional space, the 
probability of a kink being of the q = 0 or q = 2 variety is zero, and the probability 
of the kink being of the g = 1 variety is 1. 

A subtlety that has not been discussed above is that there is also the possibility 
that if we lay down Higgs fields randomly, we may get [P_, 4] # 0 (see theorem 
in Section 2.2). In this case, as described in Section 2.2, there is no static solution 
to the equations. Then the field configuration evolves toward a static configuration. 
Our discussion above assumes that such a configuration has been reached, and 
neighboring domains always have values of ® that commute. This is not completely 
satisfactory since there are time scales that are associated with the relaxation and 
these must be compared to other time scales characterizing the phase transition. 
This is why a numerical study, such as that in Section 6.7, is needed. 

To summarize, the Kibble mechanism predicts that only g = 1 domain walls 
are formed at the SU(5) x Z2 phase transition. However, we know that the stable 
variety of walls have q = 2, and hence the q = 1 walls decay into them. The 
formation of walls and the conversion of q = 1 walls into q = 2 walls during 
a phase transformation in the SU (5) x Z2 model has not been studied. However, 
these questions have been addressed in the related $5 x Z2 model as we now discuss. 
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As discussed in the last section, the q = 1 domain wall of the $5 x Z2 model 
occupies the largest volume in the space of boundary conditions but the g = 2 wall 
has least energy. Hence there is a tension between “entropy” (number of states) and 
“energy” (mass of wall). In a phase transition, based on the Kibble argument, we 
might expect the entropy to be more important. However, the higher energy walls 
q = 1 cannot survive indefinitely and eventually decay into the q = 2 walls. One 
way to study these processes is by direct simulation of the fields as a function of 
temperature [123, 6, 7]. 

The simulations are based on a Langevin equation where thermal effects are 
treated as a noise term in the classical equations of motion together with a damping 
term. For the Ss x Z2 model (Eq. (2.30)) with its four scalar fields, the equations 
are 


(a? -V’) f#t+Vit Aas =T; (6.61) 
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Figure 6.7 Energy density distribution in space at an early time at high temper- 
ature (top panel) and then at various times after the phase transition in the lower 
panels. The last panel shows that the system has relaxed into a stable lattice of 
kinks. 


where i = 1,...,4 and V; denotes the derivative of V with respect to f;. If A = 
0 =T;, these equations are simply the classical equations of motion for the fj. 
In a thermal system, we imagine that the fields are in contact with a heat bath at 
temperature T with which energy can be exchanged. Then there can be dissipation 
which is represented by the A term and thermal noise which is represented by the 
T; term. The dissipation constant A is taken to be independent of the temperature 
but the T; are stochastic and taken to be Gaussian distributed with the following 
correlation functions 


(Tix, 1) = 0, 
(T(x, OV Cy, t) = 2A T 3;j6(x — yt — t’) (6.62) 


The procedure is to solve Eq. (6.61) with any initial condition. The noise and 
dissipation eventually drive the system to a thermal distribution at temperature T. 
To mimic the phase transition, the noise is then set to zero. All of a sudden the 
system has to find a new equilibrium state. This equilibrium state has domain walls 
and these are located and tracked in the subsequent evolution. 

In one spatial dimension, the results are shown in Fig. 6.7. At high temperature 
the energy distribution is very noisy. After the phase transition, the presence of 
kinks is clear. During the evolution, some of these kinks annihilate. In the end we 
are left with a kink lattice. 
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Figure 6.8 Length in walls (denoted by N) in two spatial dimensions against time 
for the q = 2 walls in the Ss; x Z2 model (upper solid curve) and the single field Z2 
case (lower solid curve). The dashed curve corresponds to the total Ss; x Z2 wall 
length measured by counting zeros of the diagonal elements of ®(x) and hence 
includes walls with any value of q. The difference between the solid and dashed 
curves shows that the initial network consists of a large fraction of q 4 2 walls but 
then later all the walls decay into the q = 2 walls. Comparison with the Z2 case 
shows that the S5 x Z2 network decays more slowly. (The upturn at the very end 
in the Z2 case is due to the finite simulation box.) 


Similar numerical simulations have also been done in two spatial dimensions. 
The total energy in all kinds of walls is plotted as a function of time in Fig. 6.8. The 
figure also shows the energy in only the g = 2 walls as a function of time. The 
difference of these curves shows that not all walls are of the q = 2 variety at 
formation. Other kinds of walls are present immediately after the phase transition 
but they must then decay into the least massive g = 2 wall. 

As discussed in Section 2.8, the Ss; x Z2 kinks can have nodes in two spatial 
dimensions (see Fig. 2.6). So we expect a network of domain walls to form after a 
phase transition in which six or more domain walls are joined at junctions. This is 
exactly what is seen in simulations (Fig. 6.9). Another feature that is apparent on a 
closer look at the network is that there are many pairs of walls that are very close 
to each other. These pairs occur because the unstable g = 1 walls eventually decay 
into two q = 2 walls. The forces separating the q = 2 walls are exponentially small 
and so they stay close-by during further evolution. 

We have seen that the final state of the Ss x Z2 phase transition in one spatial 
dimension is a lattice of domain walls (Fig. 6.7). In one dimension, it can be argued 
that a lattice forms with unit probability provided the size of the simulation box is 
much larger than the wall thickness. In two dimensions, if the spatial extent in one 
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Figure 6.9 Network of S5 x Z2 walls in two spatial dimensions soon after the 
phase transition. The picture looks very similar to the network of (one-dimensional) 
walls connected to a network of (point-like) strings studied in [130]. 


direction is smaller than that in the other direction, so that the simulation box is 
rectangular with periodic boundary conditions, the evolution is very much like in 
one dimension and a lattice forms once again (see Fig. 6.10). Even on a square two- 
dimensional simulation box, a domain wall lattice is seen to form with a probability 
~ 0.05 [7]. 


6.8 Biased phase transitions 


The existence of domain walls relies only on the existence of discrete vacua. Then 
it is possible to imagine situations where the degeneracy of the discrete vacua 
is slightly broken (see Fig. 6.11).4 Now the probability that the higher energy 
vacuum is selected during the phase transition in some region is less than 1/2 
and the probability that the lower energy vacuum is selected is larger than 1/2. 
This process can again be simulated on a square lattice by throwing down black 
squares with probability p < 1/2. If p is very small, there are only a few black 
squares and these are disconnected from each other. So the domain walls are small 


4 Or perhaps the vacua are exactly degenerate but the likelihood of being in one particular vacuum is slightly 
larger because of the way in which the system was prepared. 
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Figure 6.10 Three stages for the domain wall network evolution in a toroidal do- 
main, with dimensions L, = 500 and L, = 150. The different shades correspond 
to the five possible charges of the domain walls (see Section 2.7). Note that in the 
bottom figure there is a pair of neighboring wall and antiwall of the same type (the 
walls just before and after the 300 mark). These later annihilate, leading to a final 
stable lattice consisting of ten walls. 


Figure 6.11 An asymmetric well in which the degeneracy of the vacua is slightly 
broken. 
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and closed. At a critical value of p, call it p., the black squares connect and the 
distribution is dominated by one infinite cluster of black squares. Then the black 
squares are said to percolate. Therefore the domain wall formation problem reduces 
to the classic problem of “percolation theory” [143, 96, 36] where we are interested 
in the critical probability and also the critical exponents that appear in various 
correlation functions as the critical point is approached. On a triangular lattice in 
two dimensions, the critical probability is known to be 0.5 and on a cubic lattice in 
three dimensions it is 0.31. The problem may even be studied on a random lattice 
as discussed in [95]. 

The analysis for biased domain walls implies that even if the potential is slightly 
asymmetric, infinite domain walls can form. For the SU (5) x Zz potential described 
in Eq. (2.5), the asymmetry is due to the cubic term with coupling constant, y. For 
small but non-zero values of y, infinite domain walls form. 


6.9 Open questions 


1. What is the number density of domain walls formed in a second-order phase transition? 
The question may need to be sharpened since the density keeps changing with time. 
Also, while the number density is important in cosmology, condensed matter physicists 
are mainly concerned with scaling laws (critical exponents) since these are expected to 
be universal. So a sharper question would be in terms of a critical exponent related to 
the number density (e.g. Eq. (6.60)). 

2. Is there a condensed matter system which gives a domain wall lattice? 

. Is there any role for domain wall lattices in (higher dimensional) cosmology? 


A Ww 


. If a domain wall lattice can be generalized to strings and monopoles, do string and 
monopole lattices form during a phase transition? 


7 


Dynamics of domain walls 


In this chapter we discuss the dynamics of kinks and domain walls first in the zero 
thickness approximation, and then briefly in the full field theory. The zero thickness 
approximation can be expected to be valid in the case when all other length scales, 
such as the radii of curvature of a domain wall, are much larger than the wall 
thickness.' We start by deriving the action for a kink in 1 + 1 dimensions as this 
is the simplest case and contains the essential features of the higher dimensional 
cases. Then we derive the action for a domain wall in 3 + 1 dimensions and some 
consequences. In this chapter we ignore gravitational effects which can be quite 
important in certain situations (see Chapter 8). 


7.1 Kinks in 1 + 1 dimensions 


In 1 + 1 dimensions, if we ignore the structure of the kink, then we expect the kink 
to behave simply as a massive point particle. Its dynamics are then given by the 
usual action for a massive point particle 


Siyi = -m far (7.1) 


where M is the mass of the kink and dr is the line element which may also be 
written as 


dX” aX” \ 1? 
) (1.2) 


eee (eu dt 


where g,,, is the metric of the space-time background and X"(t) is the location of 
the kink at time f. 

While the action in Eq. (7.1) seems reasonable on physical grounds, there should 
be a systematic way to derive it starting from the original field theory action of which 


l This expectation is not completely correct since the wall velocity is also important [183, 73]. 
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Figure 7.1 The world-line of the kink is represented by the curve. The kink frame 
coordinates yf = (t, £) are defined in the instantaneous rest frame of the kink and 
are functions of the background coordinates x“ = (t, x). 


the kink is a solution. Such a derivation should lead to Eq. (7.1) plus corrections 
that depend on the internal structure of the kink. 

To derive the effective action (Eq. (7.1)), the key assumption is that the field 
profile of the kink is well-approximated by that of the known static kink solution in 
the instantaneous rest frame of the kink. To proceed with the derivation, we work in 
“kink frame coordinates” which are denoted by y* = (t, £), a = 0, 1, as illustrated 
in Fig. 7.1 (t is also called the kink world-line coordinate). These coordinates are 
functions of the background coordinates that are denoted by x” = (t, x), uw = 0, 1. 

The kink world-line is given by the position 2-vector X” = (t, X(t)). Therefore 
the vector tangent to the world-lineis 7“ = Nr(1, 0,X) where Nr is anormalization 
factor chosen to enforce 


&wT“T”=1 (7.3) 
The unit vector, N“(t), orthogonal to the world-line is found by solving 
Suv * NY =0 (7.4) 
together with the normalization condition 
8uvN? NY =—1 (7.5) 


In the special case of a Minkowski background, guv = Nav = diag(1, —1), we find 
T” = y(1, V) where V = 0,X, and N” = y(V, 1) where y = 1/V1 — V?. 

The coordinate t is along T” and £ is along N”. Therefore, in the neighborhood 
of some fixed point on the world-line, any space-time point can be written as 


x" = X} +17) + ENG = X" (1) + EN" (1) (7.6) 
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where the subscript 0 refers to the fixed point on the world-line. Since the energy 
density in the fields vanishes far from the kink, only the neighborhood of the world- 
line is relevant for deriving the effective action. Hence & is small and to lowest order 
we can replace T in the last term by T to get 


x” = X*(t) +EN*(t) (7.7) 


With the coordinate transformation in Eq. (7.7), the world-line metric can be 
written in the y“ coordinate system 


hab = Euvõax” Opx” (7.8) 
Therefore 


hoo = 8yv(0,X" + §0,N")(0,X" + EIN”) 
= 8yv0,X"0,X" + OF) 

hoi = 8yv(0,X" + §0,N")N” = OC) 

hii = 8u NN” =-1 


where we have used the orthogonality of 9 X” « T” and N“, and the normalization 
of N”. So the determinant of hap is 


h = —gy,(X")0,X"0,X" + OE) (7.9) 


where we have also expanded the background metric around the kink location. 
Next we write, 


P(x") = boy") + WO) (7.10) 


where ¢ is the static kink profile function in the kink frame coordinates. For 
example, in the case of the Z2 kink, pọ = ntanh(./A/2 nė) (see Eq. (1.9)). The 
field y is the departure of the true field configuration from the static kink profile 
po. The assumption is that the contribution of y to the action is small and hence y 
can be used as a parameter for a perturbative expansion. 

Now the field theory action is 


f= J EREE A (7.11) 


in terms of the Lagrangian density L and g = Det(g,,,). The metric is taken to be 
a fixed background and the gravitational effects of the wall are ignored. The full 
problem of gravitating domain walls is significantly more complicated at a technical 
level [21]. 

Now we write this action in the kink frame coordinates to get 


S = So + OE, y) (7.12) 
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where M is the mass of the kink. The last equality follows since the solution is 
static and hence the Lagrangian density equals the energy density up to a sign. 
The integration of the Lagrangian density then gives the —M factor. The effective 
action is therefore the action for a point particle, simply given by the length of the 
world-line. This result can easily be extended to walls (and strings) propagating in 
higher dimensions, and the leading term in the effective action is proportional to the 
world volume. Such geometric effective actions are often referred to as “Nambu- 
Goto actions.” Even if the self-gravity of the domain wall is taken into account, the 
dominant contribution to the effective action is still the Nambu-Goto action [21]. 

The next-to-leading order terms in the effective action, denoted by O(&, y) in 
Eq. (7.12), have been discussed for domain walls in [138, 21, 73, 28], building on 
the earlier analysis for strings [57, 105, 72]. The first-order corrections in both y 
and & vanish because the field ġo is a solution of the equation of motion and hence 
the action is an extremum at ġo. The lowest non-trivial corrections come at second 
order in € and y. An alternative approach to studying domain wall dynamics has 
been developed in [9]. 

Finally we remark that the parameter t can be chosen arbitrarily. Any other 
world-line coordinate, t’(t), leaves the effective action invariant. This fact is called 
“reparametrization invariance” of the action. 


7.2 Walls in 3 + 1 dimensions 


The location of a domain wall, X“(t, ¢, x), is described by three world-volume 
coordinates y* = (T, ¢, x). Any point, x”, can now be written in terms of the “wall 
frame coordinates” (see Eq. (7.6) and Fig. 7.2) 


x" = X" (t, 6, X) FENN, 6, x) (7.14) 


where N” is the normal to the wall. 
The derivation of the Nambu-Goto action proceeds exactly as in the kink case 
of the last section and we get 


So = -o f ep h] (7.15) 
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Figure 7.2 A curved section of a domain wall is shown. The world-sheet coordi- 
nates are (T, €, x, €) while those in the ambient space (“bulk”) are (t, x, y, z). 


where o is the energy per unit area (tension) of the wall, the integral is over the 
wall world volume, h = Det(h,,), and the world-volume metric is 


hab = 8yv(X")dgX" Op X” (7.16) 


where a, b = t, ¢, x. Note that the determinant of hap is positive for the kink in 
1 + 1 dimensions and also the domain wall in 3 + 1 dimensions. 

The major difference between the kink in 1 + 1 dimensions and the domain wall 
is that the wall can be curved, and so the profile øo, which only applies to planar 
walls, does not solve the equation of motion. For example, as the wall moves, it 
accelerates and emits radiation. The radiation part must be treated as a perturbation. 
However, the analysis is conceptually the same as for the kink and the derivation 
may be found in [138, 21, 73, 28]. 

From the Nambu-Goto action for the domain wall, we can derive the equations 
of motion. The variation of So involves the variation of h = Det(hyg). This follows 
from the identity (see Appendix E) 


dInDetM = Tr(M~'sM) (7.17) 
valid for any invertible matrix M. Applying this identity to the matrix h,, we get 
dh = hh” shan (7.18) 
where h?? is the inverse of hap so that 
h®hpe = ô (7.19) 


Therefore the variation of So is 


= g. 3 ab 
So = fa py |h|h® dhap (7.20) 
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We obtain the wall equation of motion by requiring 5S) = 0 together with the 
definition of hap in Eq. (7.16) 
1 
VIAl 


where the Christoffel symbol is defined by the background metric g „v 


daly AJh XT) = Te Aer agX" 3X” (7.21) 


op 


o & 
Foy = 5 (9r8 en + Ou8pv = Oo8uv) (7.22) 


In the special case of a Minkowski background metric, the Christoffel symbol 
vanishes and 
1 
Vint 
Using Eq. (7.17), the determinant h can be eliminated and the equation of motion 
can be written as 


da(/|h|h2?a,X7) = 0 (7.23) 


1 
3a (h 3p XT) + ah ðahca h*?d,X° = 0 (7.24) 


The equation of motion for a wall is highly non-linear because hap itself is 
defined as a quadratic in derivatives of X“. One way to simplify the equations is to 
choose convenient coordinates. This is possible because the equations of motion 
of the wall are reparametrization invariant, i.e. we are free to choose any world- 
volume coordinates (t, ¢, x). A similar situation occurs for strings that havea 1 + 1 
dimensional world sheet. There, by a choice of coordinates, the equation of motion 
can be converted to a simple wave equation in 1 + 1 dimensions together with some 
quadratic constraints that can be solved quite generally. In the case of the domain 
wall, however, no such convenient choice of coordinates is known and the equations 
have not been solved in general. Only a few special solutions are known. Of these, 
static solutions subject to suitable boundary conditions have minimal surface area, 
and these have been extensively studied in the mathematics literature e.g. [115]. 

In a realistic setting, the dynamics of the walls are affected by inter-kink forces, 
by the interaction of any surrounding particles, the gravitational field of the wall, and 
the evolution of the background space-time. In addition, there are collisions between 
different walls, leading to intercommuting (Section 3.8), and annihilation of walls 
and antiwalls. If there are zero modes on the walls as described in Chapter 5, they 
could also carry charges and currents and this would introduce other interactions. 


7.3 Some solutions 


In 1 + 1 dimensions the kink moves like a point particle of mass M. The dynamics 
are richer in 3 + 1 dimensions where a closed domain wall can oscillate and move 
in complicated ways. The Nambu-Goto action is valid when the radii of curvature of 
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the wall and the separation of different sections of wall are both large compared to 
the thickness of the wall. In addition, the velocity of the wall (in the center of mass 
frame) should be small. (See Section 7.3.3 for the criterion in the case of collapsing 
spherical domain walls.) When these conditions are not met, the only way to proceed 
is to consider the dynamics using the underlying field theory. In this section, we 
ignore field theory effects and describe some solutions to the Nambu-Goto action. 


7.3.1 Planar solutions: traveling waves 


A planar domain wall in the z = 0 plane is given by 


X" (t, 6, xX) = (T, č, xX, 0) (7.25) 
Next consider a planar domain wall with some ripples 
X" (T, &, x) = (8,6, X, Z(T, f, X) (7.26) 


The function z describes the ripples and we would like solutions for z. 
For the wall in Eq. (7.26), the world-volume metric is 


hab = Nab — 9aZ pZ (7.27) 
where 


Inverting hap is not simple, but inverting nab is. So consider the “trial” inverse 
metric 


he = n”: a ne daz n° ez (7.29) 


Then by evaluating hap’, it can be seen that h”° is the correct inverse metric 
provided 


n dqzdpz = 0 (7.30) 


Now we can use Eq. (7.24) and the constraint (7.30) to get the equation of motion 
for the function z(t, ¢, x) 


0°0qz =0 (7.31) 


Hence any function that satisfies Eqs. (7.31) and (7.30) extremizes the Nambu-Goto 
action for a domain wall. 

Solutions of Eqs. (7.31) and (7.30) have been discussed in [58]. The constraint 
condition implies that the solution must necessarily be time-dependent. A class of 
solutions is obtained by noting, for example, that z = f(t — ¢) solves the equation 
of motion and also the constraint for any choice of function f. This corresponds 
to a pulse of arbitrary shape on a planar domain wall that propagates in the +x 
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Figure 7.3 Sketch of a traveling wave on a planar domain wall. The pulse propa- 
gates at the speed of light along the wall. 


direction at the speed of light. Similarly 
z = f(t £(m16 + n2x)), n? + n=l (1.32) 


is a solution for any unit vector (n1, n2). These solutions are known as “traveling 
waves” (see Fig. 7.3). 

Other solutions of the wave equation (Eq. (7.31)) are also known — for exam- 
ple, circular waves — but these do not satisfy the constraint equation and/or have 
singularities. 


7.3.2 Axially symmetric walls 


Here we look for a static wall solution in a Minkowski background. The (Cartesian) 
coordinates of the wall take the form 


X"(t, 6, A) = (Tt, R(A) cos 0, R(A) sing, à) (7.33) 
with na» = diag(1, —1, —1, —1). The wall metric is seen to be 
hap = diag(1, —R?, —(1 + R”)) (7.34) 


where R’ is the derivative of R with respect to à. The equation of motion, Eq. (7.21), 
then leads to 


d R d RR’ 
we — 0, — | —— ES eg TR (7.35) 
dà ( /1 “| dà ( /1 =) 


with the solution 


R(à) = , cosh(@A) (7.36) 
a 
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Figure 7.4 Sketch of a catenoid solution. 


where œ is a parameter, T = t, A = z and @ is the angle in cylindrical coordinates. 
Equation (7.36) describes a one-parameter family of static, axially symmetric, do- 
main wall solutions (see Fig. 7.4). 

The solution in Eq. (7.36) is a catenoid that is seen in soap films which, like 
domain walls, also minimize their surface area [23]. Experiments with soap films 
are done with two parallel circular rings, each of diameter D, placed a certain 
distance, d, apart. Then the soap film forms a catenoid for d/D < 0.66 [117]. 
Actually there are two catenoid solutions for d/D < 0.66 since the relation «D = 
cosh(ad /2) has two solutions for «œ for fixed values of D and d in this regime. A 
third solution, which consists of two disconnected disks circumscribed by each of 
the circular rings also exists. For larger values of the separation-to-diameter ratio, 
d/D, the two-disk solution has less surface area than the catenoid solutions, and 
the catenoid can pinch off and minimize its area by transforming to the two disks. 
It seems reasonable to assume that the soap film analysis also applies to the domain 
wall. 

The catenoid is a static solution of the Nambu-Goto equations of motion. It 
could happen that the catenoid is not a solution of the field equations. A simple 
example of a solution to the Nambu-Goto equations that does not solve the field 
equations can be constructed quite easily. Two parallel planar walls (a wall and 
an antiwall) form a solution to the Nambu-Goto equations but, since these walls 
have an exponentially small attractive force, they do not form a solution to the 
field equations. However, by fixing the boundary conditions (as in the soap film 
case by the rings), the catenoid solution for domain walls has been constructed 
numerically by solving the equations of motion for the scalar field in the Z} model 
(Sutcliffe, P., 2005, private communication). The stability of the catenoid solution 
to the Nambu-Goto equations is an open question (Section 7.7). 
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Quite complicated static domain wall solutions have also been studied in the 
context of quasicrystals [137] and microemulsions [70]. 

In addition to static solutions, we could seek time-dependent solutions with axial 
symmetry. The simplest such case would be a cylindrical domain wall whose radius 
is a function of time. The radius would contract, pass through zero, and then grow 
again. A similar solution is obtained for spherical walls which we discuss more 
explicitly in the next section. 

To obtain the cylindrical solution, we note that energy is conserved during col- 
lapse. The energy per unit length of a cylindrical wall is 


o2nR 
A = ——— = constant (7.37) 


ta 


where o is the energy per unit area of the wall, R is the radius of the cylinder at 
time ¢, and an overdot denotes differentiation with respect to t. The square root 
factor in the denominator takes care of the Lorentz boost. 

The conservation of energy (i.e. constancy of A), immediately leads to the 
solution 


R(t) = Ro cos (x) (7.38) 


where Ro = A/o2z is the radius when the wall is at rest. 


7.3.3 Spherical walls 


Our final example of domain wall solutions is with a spherical ansatz 
X"(t, 0, 6) = (t, R(T) (7.39) 
where 
t =t, fF = (sin 0 cos ġ, sin 8 sing, cos 0) (7.40) 


and 6,¢@ are spherical angular coordinates. The space-time metric is yyy = 
diag(1, —1, —1, —1). 
We now find 


hap = diag(1 — R*, —R?, — R° sin’ 6) (7.41) 


where overdots denote derivatives with respect to t. After some algebra, from 
Eq. (7.21) we obtain the equation of motion 


R= ae k? 7.42 
= st ) (7.42) 
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Figure 7.5 Radius of a collapsing spherical domain wall against time in the thin 
wall approximation. The coordinates in the plot are in units of the maximum radius 
of the wall. 


For R Æ 0, R? Æ 0, 1 this can also be written as 
d R? 
— | —— | = 0 (7.43) 
dt 1— R2 


Vv1— R? 
where M is a constant of motion, to be identified with the mass of the spherical 
domain wall (o is the mass per unit area of the wall). 
The solution can be written in terms of the elliptic integral of the first kind 


which implies 


f o aa On (1.45) 
Xy 1— x4 Ro 
where 
R= M genie (7.46) 
Ato Ro 


Ro has the interpretation of being the radius when the wall is at rest and x, is the 
value of x at some initial time to. The sign in Eq. (7.45) is chosen according to 
whether one is interested in the expanding or contracting solution. The radius of a 
collapsing spherical domain wall is plotted in Fig. 7.5. 

The behavior of perturbations on the spherical domain wall has been studied in 
[182]. The result is that at late times the ratio of the perturbation amplitude divided 
by the radius of the spherical wall, grows as 1/R as the wall collapses. 
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In the Nambu-Goto description, the spherical domain wall oscillates about the 
center. However, the solution is only valid as long as the thin wall approximation 
holds. By comparing various terms in the field equations of motion, the thin wall 
approximation is seen to break down when [183, 73] 


R 1/3 
ee ee (7.47) 
Ro Ro 


where w is the wall thickness. This relation is also confirmed by numerically solving 
the equation of motion in the field theory [183]. In [73], the leading order corrections 
owing to the thickness and gravity of the spherical domain wall are included, with the 
conclusion that both these effects tend to slow down the dynamics. The Nambu-Goto 
action also becomes inadequate owing to radiative losses. As the wall collapses, we 
expect energy losses owing to radiation and eventually annihilation of the domain 
wall into radiation. We discuss these processes further in Section 7.5. 

The collapse of a zero thickness spherical domain wall is prevented if the back- 
ground space-time is expanding. Static solutions are obtained if the background is 
expanding at a constant rate, as in de Sitter space. In a particular coordinate system, 
the line element for de Sitter space becomes time independent 


ds? = f(r)dt* — f7'(r)dr? — r?(d0? + sin? 0de?) (7.48) 


where f(r) = 1 — H’r? and H is a constant corresponding to the expansion rate. 
Following the analysis of [16] for a circular string, the action for a spherical domain 
wall in the zero thickness limit is 


S=-4n0 J dtR? | f — S (7.49) 


where R(t) is the radius of the spherical wall and f = f (R). Extremization of this 
action leads to the first integral 


R? — f? +e? RIF? = R? +V(R)=0 (7.50) 


where e = E /4xro and E is a constant (the first integral). For a static solution we 
need both V (R) = 0 and V’(R) = 0 where prime denotes derivative with respect 
to R. These conditions give the static solution 


R= a? (7.51) 
= : ; 


_ 40 2 
= a NE 


with 


(7.52) 
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The potential V(R) is a maximum at the location of this solution and therefore 
the solution is unstable. The instability can be understood without calculation. If the 
radius of the wall is perturbed to be a little smaller than the value at the solution, 
the effects of Hubble expansion are weaker while the force owing to curvature is 
stronger, and so the wall collapses. On the other hand, if the radius is perturbed to 
be a little larger than the solution value, the expansion effect is stronger while the 
curvature force is weaker, and the wall expands to yet greater radii. 

Planar and spherical domain walls in de Sitter space have been considered in the 
full field theory in [17, 18]. It is found [17] that instanton solutions describing the 
nucleation of spherical domain walls exist only when the thickness of the wall is 
less than H~!/,/2. This result is also relevant to the problem of finding static thick 
spherical domain walls in de Sitter space, since an instanton solution can exist only 
if the static domain wall solution exists (though the converse may not hold). Hence 
spherical domain wall solutions of the field theory in de Sitter space exist if the 
domain wall thickness is less than H~!//2. 


7.4 Solutions in field theory: traveling waves 


The traveling wave solutions discussed in Section 7.3.1 in the zero thick- 
ness approximation are also exact solutions to the field equations of motion 
[160, 161]. 

Consider the field 


olt, x) = polz = Zo(t, x, y)) (7.53) 


where ¢(z) is the classical solution for a domain wall in the z = 0 plane. We now 
insert this ansatz in the field theory equation of motion. A little algebra shows that 
the ansatz is a solution provided 


3aðzo =0, (8aZ0)* = 0 (7.54) 


where a = t, x, y. These are the same equations obtained above for planar solu- 
tions to the Nambu-Goto equations (Eqs. (7.31) and (7.30)). As discussed there, 
the only non-singular solutions to these equations have the form of traveling 
waves e.g. 


zo, x, y) = f(t +x, y) (7.55) 
Hence traveling waves are solutions to the field equations and do not dissipate 


owing to radiation.? 


2 Tt can be shown that traveling waves do not dissipate even when they are considered in quantum field theory 
[46]. 
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Figure 7.6 Collapse of a spherical sine-Gordon domain wall. The curves in the 
left-hand plots show the field as a function of radial distance for several different 
times. The right-hand plot shows the corresponding energy density distributions. 
[Figure reprinted from [183].] 


7.5 Spherical domain walls: field theory 


We have seen in Section 3.7 that the collision of a kink and an antikink in 1 + 1 
dimensions leads to chaotic dynamics. The kinks bounce back for certain velocities 
while for other velocities, both smaller and larger, they annihilate. So we might 
expect the dynamics of a collapsing spherical domain wall to show similar features. 
Numerical simulations of the sine-Gordon model show that a collapsing spherical 
domain wall does not radiate very much energy until it becomes very small (of order 
the thickness of the wall), then emits a large amount of radiation, then bounces back 
to form an expanding spherical domain wall (though with less energy than the initial 
configuration), which then reverses and collapses again (see Fig. 7.6). Simulations 
of a A¢* spherical domain wall, however, do not show any bounce back [183]. 


7.6 Kink lattice dynamics (Toda lattice) 


In Section 6.6.2 we have seen that a phase transition can lead to the formation of a 
lattice of kinks (Fig. 6.7). What happens if one of the kinks in a lattice collides with 
a neighboring kink? The interaction potential between neighboring kinks decays 
exponentially with distance and energy conservation implies that the collision is 
perfectly elastic. The momentum of the incoming kink is transferred to the target 
kink [123]. These properties are exactly those assumed for a chain of masses in 
what is called a “Toda lattice” [155]. The many beautiful properties of a Toda 
lattice apply to the (one-dimensional) lattice of kinks as well. For example, there 
are soliton solutions that run along the Toda lattice. So there are also solitons in the 
dynamical modes of the kink lattice i.e. solitons in the dynamics of solitons! 
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7.7 Open questions 


. Are there closed domain walls in three dimensions that do not self-intersect as they 
oscillate? What happens in higher dimensions? 


2. Can one show analytically that walls must intercommute on intersection? 


. When traveling waves on domain walls collide, they dissipate some of their energy. Find 
the energy radiated. Find the energy that goes into excitations of the bound state in the 
case of the Z> wall. 

. Analyze the catenoid domain wall solution and its stability. 
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Gravity and cosmology of domain walls 


Domain walls resulting from a symmetry breaking in the early universe could have 
novel and dramatic gravitational and cosmological consequences. 

We first derive the gravitational effects of a planar domain wall, describing the 
different ways to view the system. Then we discuss spherical walls as an example 
of curved domain walls. To discuss the cosmological consequences, it is necessary 
to have a picture of domain wall formation in the cosmological context. With 
the background of Chapter 6 we discuss the formation of the wall network in 
cosmology, then the evolution and cosmological implications. We end by reviewing 
the cosmological constraints on domain walls and the few possible ways around 
the constraints. 


8.1 Energy-momentum of domain walls 


The energy-momentum tensor for a scalar field with potential V(@) is given in 
Eq. (1.39) 


1 
Tv = .69)6 — TE -— væl (8.1) 


In the thin-wall limit, varying the Nambu-Goto action (Eq. (7.15)) gives the energy- 
momentum tensor 


re |h| nea, X IX” SO(x# — X”) (8.2) 


dp 
To 
where X“ is the location of the wall. For a planar wall located at x = 0 in flat 
space-time, this gives 


TE = 0(1, 0, —1, —1)d(x) (8.3) 


NG, plane 
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For a planar wall, including self-gravity, the energy-momentum tensor can be ex- 
plicitly written once we have chosen a suitable ansatz for the metric (see Eq. (8.4) 
below). 


8.2 Gravity: thin planar domain walls 


The gravitational effects of a planar domain wall have been found in the thin-wall 
limit in [167, 169, 80] using the vacuum solutions found in [154]. The thin-wall 
limit simplifies the analysis because then there is no need to solve the field equations 
of motion. All the energy-momentum is localized on the thin domain wall and so 
only the vacuum Einstein equations need to be solved on either side of the wall. 
The presence of the wall shows up in matching the vacuum solutions on the two 
sides of the wall i.e. implementing the “junction conditions.” Such a matching is 
facilitated by using the Gauss-Codazzi formalism [81] and this has been done in 
[80]. Here we derive the metric of a domain wall without going through the general 
Gauss-Codazzi formalism, following the derivation in [169] instead. 

A planar domain wall located in the x = 0 plane has rotational symmetry in 
this plane. Further we expect space-time symmetry under x — —x. Under these 
conditions the form of the line element can be taken to be [154] 


ds? = e™ (+d? — dx”) — e” (dy? + dz’) (8.4) 


where u and v are functions of ¢ and |x|. Note that the possibility that the metric is 
time-dependent has been retained. 

In the thin-wall limit, there is no energy-momentum off the wall and so the 
energy-momentum tensor, T», vanishes everywhere except on the wall. Therefore 
only the vacuum Einstein equations, R,,, = 0, where R» is the Ricci tensor, need 
be solved. The solution for x > 0 is 


e™ = f(t +x)+ g(x -— t) (8.5) 
1 
n= AEE ERA ET k = (8.6) 
where the functions f, g, h, and k satisfy 


fe x g” = 2 f'h' + 2g'k' Z 0 (8.7) 
f” + g" =, 2 = 2g'k' = 0 (8.8) 


where primes denote derivatives with respect to x. The solution for x < 0 can be 
obtained by symmetry since u and v are functions of |x|. 
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Next we solve the Einstein equations, Tay = G,,/8"7G, where G,, is the 
Einstein tensor calculated for the metric in Eq. (8.4). This leads to 


1 
T? = —— yves x 
0 G” (x) 


AnG 

T}=0 

T? = T? = ae + ve s(x) (8.9) 
2 3 8G 0 0 t 


where up = u(t, x = 0+), vo = v(t, x = 0+). 

In general, uo, Up; and Up are time-dependent, and so these expressions for T” are 
also time-dependent. However, the energy-momentum tensor for the wall should 
be time-independent. This gives us the constraint that the functions f, g, h, and k 
must be chosen so that wo, up, and up are time-independent. Then the only possible 
choice for the functions (for x > 0) that also satisfy Eq. (8.8) is 


f = 0, g= e560 


h = ee ) py ) (8.10) 
7 x), z x : 
where 
K = 4r Go (8.11) 
The corresponding functions for x < 0 are 
f = e5 0t») g=0 
K K 
h = —(t + x), k = —— (t — x) (8.12) 
2 4 
Then the domain wall line element is 
ds? = e7¥ Ml [dt? — dx? — e¥' (dy? + dz°)] (8.13) 


which can also be put in the commonly encountered form 
ds? = = KX) d — dX? — (1 — K| X) e (dy? + dz (8.14) 
where k = 2x Go via the coordinate transformation 


1x] = 4d a, (8.15) 
K 


8.3 Gravitational properties of the thin planar wall 


On spatial slices of constant X (X = Xo) the metric of Eq. (8.14) takes the form 
ds? = dP? — e™' (d5? + dz’) (8.16) 
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where overbars denote that the coordinates have been rescaled by the factor 
(1—x«|Xo|) and k =«/(1—«|Xo|). The three-dimensional line element of 
Eq. (8.16) shows that space-like slices of constant X are expanding exponentially 
fast, just as in an inflationary space-time. 

The inflationary nature of the metric can be understood from the viewpoint of 
an observer living on the wall who is blind to the coordinate normal to the wall. 
From such an observer’s perspective, the space-time is filled with vacuum energy, 
as given by the energy-momentum tensor of Eq. (8.3), and hence is inflating. 

Next we examine the metric on spatial slices obtained by setting y = yo, z = Zo 


ds? = (1 — x|X|)"dt? — dx? (8.17) 


This is the metric of 1 + 1 dimensional Rindler space-time, which is Minkowski 
space-time written in the rest frame coordinates of a uniformly accelerated observer 
with acceleration a = 1/« away from the wall which is located at X = 0. To see 
this, use the coordinate transformation 


t= a — XD e = e “') 
2K 

p= SID a + e *‘) (8.18) 
2K 


In these coordinates the Rindler line element is of Minkowski form 
ds? = dt? — dé? (8.19) 


Now note that 


2 
PoP = (G - x1) (8.20) 
K 


Therefore the world line of a particle at fixed X is ahyperboloid in Minkowski space- 
time, which describes a particle moving at constant acceleration. In particular, the 
wall located at X = 0 has acceleration 1/«. Therefore an inertial observer sees the 
wall accelerating away with acceleration 1 /«. From the perspective of an observer 
on the wall, all particles are repelled from the wall. 

In the Rindler space metric there is a horizon at |X| = 1/k. It is clear from the 
coordinate transformation given above, this is a coordinate singularity since the 
space-time is equivalent to Minkowski space-time. 

As discussed in [80] the full domain wall metric in Eq. (8.14) can also be brought 
to Minkowski form. This shows explicitly that the domain wall space-time is flat 
everywhere except on the wall itself. As in the reduced metric of Eq. (8.17), in the 
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Minkowski coordinates (tw, xm, ym, Zm) the wall is located at (see Eq. (8.20)) 
1 
xu + yu t+ ZM = tw + 5 (8.21) 


Hence, in the coordinates where the metric is Minkowski, the wall is spherical with 
time-dependent radius that decreases (for ty < 0) until it gets to 1/« at ty = 0 and 
then bounces back. This behavior does not depend on which side of the wall the 
observer is located. Both see the wall accelerating away from them with constant 
acceleration 1/k. 

There is an intuitive way to see that the wall’s gravity must be repulsive. In the 
weak field approximation, the gravitational potential of the wall is proportional to 
p + pi + p2 + p3 where p is the energy density of the wall and p; are the pressure 
components of the energy-momentum tensor. From the energy-momentum tensor 
in Eq. (8.3) we have p = p3 = —p and pı = 0. Therefore o + pi + p2 + p3 = 
—p < Oinstead of the positive value obtained for matter without pressure. Therefore 
the gravitational potential is repulsive instead of being attractive. 

Since the metric is Minkowski in the (tm, xm, ym, zm) coordinates, geodesics are 
given by 


xt (ty) = xË +u” (tm — to) (8.22) 


where, xj is the position of the particle at time fy = fo, and u” is the (constant) 
velocity vector. 


8.4 Gravity: thick planar wall 


Here we consider the gravitational field of a thick domain wall i.e. taking both the 
scalar field and Einstein equations into account. 
The Einstein equations are 


Guy = 82 GT yy 
- sa| a,oa - su| (ud? - vo (8.23) 
where we have used T, from Eq. (8.1). The scalar field equation is 
V V“ + V(b) =0 (8.24) 


where V, is the covariant derivative. 

These equations have been solved in [181] for the case when 167Gr* «1 
where 7 is the vacuum expectation value of the field ¢. The line element outside 
the thick wall is still given by Eq. (8.14) and there are no qualitative new effects. 
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The case when 16r Gn? > 1, however, does lead to new effects as first discussed 
in [170, 101, 102] and as summarized in the next section. 


8.5 Topological inflation 


If 167 Gn? > 1, the gravitational forces within the wall are stronger than the forces 
associated with the self-interaction of the scalar field. This can be seen by the 
following heuristic argument [170]. 

Consider the Z2 model with the quartic potential 


À 
V$) = ria -y (8.25) 


The thickness of the domain wall can be estimated by equating gradient and potential 
energies, which also agrees with the Bogomolnyi equation (see Eq. (1.31)), in the 
case when gravitational effects are ignored. The field @ gets an expectation value 
n and so, in the interior of the domain wall, 

2 


Lovo = 2 ~ VO) = Èr (8.26) 
2 282 4" 
and the thickness, ô, is 
5 é (8.27) 
àn ` 


This is an estimate of the length scale on which the scalar field interactions are 
working. 

Next, the length scale associated with gravitational effects is found from the 
Friedman-Robertson-Walker equation, which relates the space-time expansion rate, 


H, to the energy density 
8G 
H? ~ = p (8.28) 


which, when used inside the wall with p ~ An*/2, gives 


4 3 1 
SA ae Gap P 


Hence scalar field forces dominate over gravitational forces inside the domain wall 
if H7! > ô, or the order of magnitude condition, 


16xGn? < 1 (8.30) 


Therefore when 167 Gn? < 1 we can expect that gravitational effects are small 
in the interior of the domain wall. If, however, H7! < ô, the field is approximately 


134 Gravity and cosmology of domain walls 


smooth over a region where gravitational effects are strong. The field inside the 
domain wall has potential energy ~ 474 and this is what drives the gravitational 
effects. Therefore, we expect that the space-time inside the domain wall inflates 
in the direction normal to the wall, in addition to the inflation parallel to the wall 
that we have already seen in the thin-wall case (Eq. (8.14)). Furthermore, the field 
inside the wall is stuck on top of the potential owing to the topology that led to 
the existence of the wall. So the inflation goes on forever for topological reasons. 
Hence this inflation is called “topological inflation.” 

This picture has been confirmed by numerical solution of the coupled scalar 
field and Einstein equations in [131, 31, 94] with the conclusion that topological 
inflation inside the Z2 domain wall occurs for n > 0.33mp where mp is the Planck 
mass defined by G = 1/m3. 

Domain walls that are undergoing topological inflation cannot however form in 
the usual way during a cosmological phase transition as we discuss in Section 8.9 
below. 


8.6 Spherical domain wall 


The metric of a thin spherical domain wall has been discussed in the thin-wall limit 
in [80]. Inside the wall the metric is flat using Birkhoff’s theorem (e.g. see [177]) 


ds? = dT? — dr? — r?(d6* + sin? dg’), r < R(t) (8.31) 
where R(t) is the radius of the spherical wall and 
T =(1+ R%)!? (8.32) 


with overdots denoting derivatives with respect to the proper time of an observer 
moving with the domain wall. The proper time is related to the time coordinate t 


via the relation 
2GM\. IGM.  .,\ 17 
(1 = a ii = (1 e R?) (8.33) 


Outside the sphere, the metric is Schwarzschild with mass parameter M 


2GM 2GM\"! 
ds? = (1 — EM — (1 > cM) dr? — r?° (d0? + sin? 0do?®), r > R(t) 
r r 
(8.34) 
The mass is related to the maximum radius of the spherical wall, Rm, by 


M = 4ro R? (1 — 27 Go Rm) (8.35) 
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Figure 8.1 Energy in a large volume enclosing a collapsing spherical domain 
wall in the sine-Gordon and Z% models [183] as a function of time. The energy is 
roughly conserved until the radius becomes comparable to the wall thickness and 
then decreases sharply. The step-like features in the sine-Gordon models occur 
because the sphere bounces several times before annihilating. The spherical wall 
in the Z2 model annihilates without bouncing. [Figure reprinted from [183].] 


provided Rm < 1/4r Go. If Ra > 1/4r Go, it means that the spherical domain 
wall is a black hole even at the maximum value of its radius and the analysis breaks 
down. 


8.7 Scalar and gravitational radiation from domain walls 


A collapsing spherical domain wall emits scalar radiation and loses energy. It may 
be possible to extend the formalism in Section 3.5 to calculate this energy loss. 
However, such an analysis is not currently available. Instead the energy emission 
rate has been found numerically and is shown in Fig. 8.1 for spherical walls in the 
sine-Gordon and Z models [183]. 

A collapsing spherical domain wall does not emit gravitational radiation since 
the spherical symmetry implies a vanishing quadrupole moment of the energy- 
momentum distribution. However, colliding domain walls can lead to gravitational 
[157] and scalar radiation [175]. A dimensional analysis based on the quadrupole 
approximation for the gravitational power emitted when two relativistic spherical 
walls collide gives [157] 


_ GM 
R2 


P, (8.36) 
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where Mg ~ 47o R? is the mass of the bubble and R is the radius upon collision. 
Numerical analyses of bubble collisions (during first-order phase transitions) found 
that the quadrupole approximation overestimates the power radiated in gravitational 
radiation by about a factor of 50 [91]. 


8.8 Collapse into black holes 


If the radius, R(t) of a collapsing spherical domain wall remains larger than the 
Schwarzschild radius, Rs = 2GM, where M is the mass of the domain wall, then 
the domain wall does not become a black hole. As the wall collapses, it emits 
scalar radiation and, if this is rapid enough, M decreases sufficiently rapidly so 
that Rs < R at all times. Whether this happens can be checked explicitly by nu- 
merical evolution of the scalar field plus Einstein equations. We expect that if the 
Schwarzschild radius of the spherical domain wall is smaller than the width of the 
wall, black holes are not formed since rapid wall annihilation and radiation precede 
collapse to within the Schwarzschild radius. 

The converse case of black hole formation in the case when the scalar radiation 
is not too rapid is harder to demonstrate convincingly. The reason is that the time 
evolution of the fields gets slower as the black hole event horizon is about to form. By 
simply evolving the fields, it is impossible to see the formation of the event horizon 
and hence conclude that the domain wall collapses to form a black hole. However, it 
is hard to imagine any other outcome, especially since the scalar radiation rate only 
becomes significant once the spherical domain wall collapses to a size comparable 
to the thickness of the wall. 

The collapse of a slightly perturbed spherical domain wall has been studied 
numerically in [182] with the result that the amplitude of perturbations stays con- 
stant during the collapse. This means that the ratio of the perturbation amplitude 
to the radius grows during collapse as 1/R(t) and the shape of the wall deviates 
increasingly from being spherical. 


8.9 Cosmological domain walls: formation 


The formation of domain walls in a phase transition in flat non-expanding space- 
time has been discussed in Chapter 6. Since the universe is expanding and cooling, 
cosmic phase transitions can occur, just as in the laboratory, and domain walls can 
also form. If the phase transition proceeds quickly on cosmological time scales, the 
structure of these domain wall networks is similar to those formed in the laboratory 
and described in Chapter 6. The network is dominated by one infinite domain wall 
with very complicated topology. However, if the phase transition occurs slowly 
on cosmological time scales, the expansion can prevent the phase transition from 
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completion. For example, in a first-order phase transition, if the bubble nucleation 
rate is very slow, the bubbles will not be able to percolate because the expansion 
increases the separation of the bubbles that have already nucleated. These consid- 
erations are important for inflationary cosmology but here we will assume that the 
phase transition completes since otherwise domain walls would not be formed. 

In a model with 1627 Gn? > 1 (n is the vacuum expectation value of the scalar 
field), domain wall formation requires some new considerations [22, 166]. The 
reason is that the energy density inside the domain walls is larger than that outside. 
Hence if such inflating domain walls (see Section 8.5) were to form, the space-time 
expansion rate within them would be greater than that of the ambient cosmological 
expansion rate in which they were created. It is possible to show that a faster 
expanding region within a horizon of a slower expanding region can be created 
only if the null energy condition! is violated. The formation of defects proceeds 
according to the classical dynamics of a scalar field with energy-momentum tensor 
given by Eq. (8.1). Contracting the energy-momentum tensor twice with a null 
vector, N“, and using g,,N"“N” = 0 gives 


NENT, = (N3 p) = 0 (8.37) 


Hence the null energy condition is satisfied during defect formation and there is an 
obstruction to the formation of topologically inflating domain walls. The exception 
is if the faster expanding region has an extent that is larger than the cosmological 
horizon. In this situation, the domain wall is fatter than the horizon during the 
phase transition. Then the particle interaction rate is also slower than the Hubble 
expansion rate and the particles are not in a thermal state unless they were set up 
in that state as an initial condition at the Big Bang. The domain wall network that 
is produced will depend on the initial state of the particles. 


8.10 Cosmological domain walls: evolution 


If we assume that there is a dense network of walls within our cosmological horizon 
and that the network does not lose a significant amount of energy to radiation, we 
can work out the expansion rate of the universe and the scaling of the density of 
walls. 

The energy-momentum of the scalar field that forms the domain walls is given 
in Eq. (8.1). If we denote an average over a large volume by (-) we have 


1 
(Tuv) = (0,60,) na (emf 50? md væl] (8.38) 


1 The null energy condition is N” N "Tav > 0 where N” is any null vector and T,» is the energy-momentum 
tensor. For fluids with energy density p and isotropic pressure p, the null energy condition is p + p > 0. 
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We assume that g,,, is a background metric and only dependent on time. Also the 
field distribution is assumed to be isotropic so that 


Cap) = yp) = ((826)") (8.39) 
and 
(3;09;¢) =0, ix j (8.40) 
Define 
1 
(p°) = 3 (xp) + (8y6)° + (820V) (8.41) 


If we further assume that the field is dominantly in the form of domain walls that 
satisfy Eq. (1.31) to a good approximation, we get 


12 2 


(Q) = 3V) (8.42) 
which leads to 
ie R (8.43) 
XX. = 6 (0 a tt . 


For slowly varying fields this leads to the effective equation of state p = —2/3 
where p is the (isotropic) pressure and p the energy density [186]. If we assume 
that the time dependence of ¢ is only due to a boost of the domain walls, we can 
use @ = vy dx = vy /2V(¢) and 0,¢ = ydx@, where X = y(x — vt) and y is 
the Lorentz factor (see Eq. (1.10)). This leads to 


iy 2 
(Tix) = (v p= =) T) (8.44) 


Following Appendix F and treating the wall network as a fluid with equation 
of state p = —2p/3, we can write down the solutions for the scale factor and the 
scaling of the energy density in walls. If the initial conditions are such that the wall 
density is pọ when the scale factor is ao, the solution is 


ao t i 
Pwalls(a) = po—, a(t) = ao( =) (8.45) 
a to 


Note that this derivation ignores processes by which the wall network could lose 
energy into scalar and gravitational radiation. In addition, the walls interact with 
surrounding matter and experience friction. These effects make the problem of 
understanding the evolution of the wall network much more challenging. We now 
describe some numerical [125, 36, 97, 59] and analytical [77, 78] efforts to under- 
stand the evolution of the network. 
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8.11 Evolution: numerical results 


There are two numerical schemes for evolving a network of domain walls. The 
first is to use the zero thickness approximation for walls. In this approximation, it 
is hard to treat the collision of walls and the loss of energy from the network into 
radiation. The second approach is to solve the field theory equations of motion. In 
this approach, all the degrees of freedom of the system are retained. In fact, a lot 
of degrees of freedom that are evolved are inessential to the domain wall network 
and this additional baggage slows down the simulations. In an expanding universe 
the problem is even more severe because the overall length scales grow larger with 
time while the domain wall thickness remains the same. Thus the simulation needs 
to handle very disparate length scales. 

In [125, 36, 97, 59], the authors get around these problems by solving the field 
theory equations of motion but by letting the domain walls expand with the universe. 

More specifically, consider the Z2 model in an expanding space-time with metric 


uv = a(t) Nav (8.46) 


where n,» = diag(1, —1, —1, —1) and t is the conformal time. The equation of 
motion is 


ao + 2" ang — Ve +A? —a’n’)o =0 (8.47) 


In the approach pioneered in [125] the a?n? in the last term is replaced by a constant, 
effectively decreasing the vacuum expectation value, 7, with Hubble expansion. 
Since the width of the domain wall is proportional to 1/7, this amounts to letting 
the thickness of the walls grow in proportion to the scale factor. 

The result of this numerical study shows that the areal density, A, i.e. area of 
walls in a given region divided by the volume of the region, scales inversely as the 
first power of conformal time 


A= A (=y -pEi (8.48) 


where the subscript 0 refers to some initial time. This result holds in Minkowski 
space-time (a œ 1°), radiation-dominated (a « t!/*), and matter-dominated 
(a « t!/3) cosmologies. 

The domain wall network has also been studied by a combination of numerical 
and analytical techniques that use scaling arguments [13, 14]. 


8.12 Evolution: analytical work 


An analytic technique to study the evolution of non-relativistic interfaces in the 
condensed matter context was developed in [116] (also see [24, 65, 111]). The 
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technique has been extended to relativistic systems in [77, 78] and we now sum- 
marize the main features of this analysis. 

The starting point is to define a fictitious scalar field u(x“) such that it vanishes 
on the domain wall network 


Ux (oe )) =0, a=0,1,2 (8.49) 


where the domain wall network is located at X“(o“) and of denote world-volume 
coordinates. While u(x“) could have been taken to be the scalar field in the original 
field theory (say for the Z2 model), this is not suitable since u is later assumed to be 
a random field with a Gaussian distribution. The next step is to derive an equation 
of motion for u. 

We define the domain wall world-volume metric as in Eq. (7.16) 


hab = Suv(X)0qX" 0p X” (8.50) 


where g,,, is the ambient space-time metric and the indices a, b refer to world- 
volume coordinates. Two derivatives of Eq. (8.49) lead to 
1 

VIA] 
As long as the thin-wall limit is valid and, in particular, walls do not intersect, 
X“ satisfies the Nambu-Goto equation. When walls do intersect, the Nambu-Goto 
formalism breaks down. The formalism can continue to be valid provided we impose 
additional boundary conditions by hand at the intersection point. Depending on 
the boundary conditions that one imposes at the intersection point, the Nambu- 
Goto equation can describe intercommuting walls or walls that pass through each 
other. In the present formalism, the boundary conditions automatically arise from 
the evolution of the u field. The dynamics of the u field are such that they always 
describe walls that intercommute [77]. Using the Nambu-Goto equations of motion, 
Eq. (7.21), then leads to the equation of motion for the fictitious field u 


ða(VIh|h 3p X")ð u + herd, X" pX” d,,d,u = 0 (8.51) 


(3u) g — 3”“uð”u](ð ðu — rE 0) =0 (8.52) 


where T£, is the Christoffel symbol defined in Eq. (7.22). 

To solve Eq. (8.52) we must find a way to handle the non-linear terms. The key 
point now is that the domain wall network contains a random distribution of walls 
and hence u is a statistical field. One approach to treat the non-linear terms is to 
use the mean field approximation. In this approach non-linear terms are replaced 
by averages of non-linear terms multiplied by a single power of u. For example 


u? —> (u?)u (8.53) 


Further, the distribution of u is assumed to be Gaussian. 


8.13 Cosmological constraints 141 


After defining the correlators that enter the mean field theory version of 
Eq. (8.52), the field u satisfies the equation of motion 


u(t) 


3u + —— au — vV-u =0 (8.54) 


where, as in the previous section, t is the conformal time, The functions u and 
v? are defined in terms of the assumed two-point correlation functions of u (for 
details see [77, 78]). Once the solution for u is obtained from the linear differential 
equation, Eq. (8.54), the average areal density, A, and other quantities may be 


calculated. The results agree with the scaling in Eq. (8.48). 


8.13 Cosmological constraints 


The cosmological constraint on domain walls is remarkably robust, being almost 
independent of the field theory, details of the phase transition, and cosmology [186]. 
At any time after the domain wall forming phase transition, the vacua in different 
cosmological horizons are uncorrelated. This means that there is at least one domain 
wall per horizon. The minimum area of a horizon size domain wall is ~ H~? where 
H`! is the horizon size. Therefore the domain wall energy density averaged over 
a horizon volume is Pwais; ~ o H. Comparing this to the critical density of the 
universe,” we get 


Ge pe CU a OS eee (8.55) 

Pc H 
where t is the cosmic time. (We have taken H ~ 1/t which is true in a Friedman- 
Robertson- Walker cosmology in which a(t) « t* with 0 < œ < 1.) Hence, as time 
proceeds, there comes an epoch when the domain walls are the dominant form of 
energy in the universe. This happens at time t, given by 

1 

by Ge (8.56) 

Now o ~ 7} (e.g. Eq. (1.20)) up to factors of coupling constants which we assume 


are order unity. We also know particle physics fairly well up to an energy scale of 
about 100 GeV (approximately the electroweak scale) and have not seen any scalar 
fields yet. So the minimum value of ø is about (100 GeV)?. Walls of this tension 
would have started dominating the universe at (see Appendix A for numerical 
values) 


2 
| ~ 2 ~ 108s (8.57) 


min 


2 The critical density of the universe is defined as pe = 3H?/82G, where H(t) = a/a is the Hubble expansion 
rate defined in terms of the scale factor a(t) and its time derivative, å. 
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or approximately 10 years after the Big Bang. Once the walls dominate, the universal 
expansion becomes a « t? (Eq. (F.5)). This is unacceptable for several reasons. For 
example, since the domain wall dominated universe accelerates (ä > 0), density 
perturbations that are larger than the horizon keep getting stretched and stay larger 
than the horizon. This means that super-horizon density perturbations can never re- 
enter the horizon, which is an essential condition for them to start growing to form 
the galaxies, clusters, and large-scale structures that we currently observe. Even 
the growth of sub-horizon density perturbations is suppressed owing to cosmic 
acceleration. 

A second constraint on a network of cosmic domain walls acting as a fluid 
with equation of state p = —2p/3 comes from the measured expansion rate of the 
universe using supernovae data [127, 118]. These surveys find that the equation 
of state parameter, w = p/p, for our universe is less than about —0.8 [11, 128]. 
However, a universe dominated by a network of static (“frustrated”) domain walls 
[25] would have w ~ —0.67. 

Another possibility that has been considered is that perhaps there are some 
features that are missing in the standard model of particle physics, and that there 
indeed are very light domain walls in the universe [76]. Such light walls, if light 
enough, would be benign and could potentially play a role in cosmology. If we 
require that the domain walls not dominate the universe until the present time 
(~ 10!” s), Eq. (8.57), gives n < 100 MeV. Other cosmological constraints, such 
as arising from the isotropy of the cosmic microwave background can be used to 
put similar or somewhat stronger bounds on n [144, 158]. 


8.14 Constraints on and implications for particle physics 


Let us summarize the picture that has emerged in this chapter. 


If a field theory has discrete symmetries that are spontaneously broken in the 
ground state, it must contain domain wall solutions. 

If high-energy particle physics is described by such a field theory and the discrete 
symmetry gets spontaneously broken in the early universe, cosmic domain walls 
are produced. 

If the standard model is complete at energies below 100 GeV, then there can be no 
domain walls in the universe and no spontaneously broken discrete symmetries 
in particle physics. 


A closer examination of this sequence of arguments reveals a few loopholes that 
allow for spontaneously broken discrete symmetries in particle physics. First, there 
is the possibility that the discrete symmetry was broken right from the moment of 
the Big Bang. Then the whole universe could have been in one of the many discrete 
vacua at its very creation and no domain walls would be formed even though the 
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Figure 8.2 Sketch of the potential in the model of Eq. (8.58). Domain walls arise 
owing to a 27 change in the angular field variable and the location of the field 
inside the wall is marked by D. Such domain walls can terminate on strings, and 
the field within the string is located at ¢ = 0. 


underlying particle-physics theory could have broken discrete symmetries. This 
kind of scenario has been studied for magnetic monopoles in [50]. A related possi- 
bility is that, if the universe went through a period of superluminal expansion (“‘in- 
flation’”),* then correlations extend to scales that are vastly larger than our current 
horizon and our region of the universe is very likely to be free of any domain walls 
[97]. In another variant, domain walls are formed but subsequently inflated away. 

All the above loopholes only apply to very high-energy domain walls where 
quantum gravity and/or inflation effects are relevant. If a particle-physics model 
has spontaneously broken discrete symmetries at lower energy scales (but still larger 
than ~ 100 MeV) no loopholes are known and the model is ruled out based on the 
“cosmological domain wall catastrophe.” However, there is still the possibility that 
metastable or biased domain walls (see Section 6.8) can exist for some time in the 
universe. We now describe these two possibilities. 


8.15 Metastable domain walls 
In certain field theories, it is possible for domain walls to get punctured. To see how 
this can happen, consider the potential for a complex scalar field @ 
À an 
V@)= (ler =y = an? + 9°) (8.58) 


where we assume 0 < œ << à. The shape of this potential is shown in Fig. 8.2. The 
first term is minimized when |¢| = 7 and, restricting ¢ to the submanifold |¢| = n, 


3 Inflation occurs when the universe is dominated by a field that has an equation of state with —p < p < —p/3. 
Then the expansion rate of the universe is superluminal and volumes that are larger than the horizon can get 
correlated. 
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Figure 8.3 Cross-section of a wall that terminates on a string is shown on the 
left, and a wall with a puncture bordered by a string is shown on the right. 


the second term is minimized when ¢ + ¢* = +2n. Another way of writing the 
potential is by setting ¢ = y exp(ix) and then yy, x are real fields. Then 


À 
Vb) = SV? — Wy — By cos? x (8.59) 
The extrema of V are at Y = 0 and at 


3a + V9a? + 6422 
y=n , X =n 


81 (8.60) 


where n is an integer. The true vacua occur when n is an even integer. For ex- 
ample, domain wall solutions exist with the boundary condition x(x = —oo) = 0, 
x(x = +00) = 27. 

Now consider a domain wall in the model in 3 + 1 dimensions. Such a domain 
wall can terminate as shown in Fig. 8.3 since the path from x = 0 to x = 27 can 
be contracted by lifting it over the top (Y = 0) of the potential. While we have 
not described cosmic strings here, in models such as Eq. (8.58), we can have finite 
sections of open walls that are bordered by strings. Walls can also get punctured 
by holes that are bounded by strings. For further discussion of walls bounded by 
strings, we refer the reader to [88, 168, 171]. 

The evolution of a network of walls that can have punctures is very different 
from that of stable walls because a puncture can grow and eat up the wall. This 
provides a very efficient way for the wall network to lose energy and so the network 
never dominates the universe [168]. 

Another scheme that allows for the universe to have a finite period of time with 
domain walls is if a discrete symmetry is broken and then restored (see Section 6.1). 
Walls would be formed at the first phase transition and then they would dissolve at 
the second phase transition when the symmetry is restored. However, this scheme 
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would imply an unbroken discrete symmetry in the low-energy particle physics. 
We do not know of such a discrete symmetry although the possibility cannot be 
excluded. 

Finally, domain walls could have existed for some time in the early universe 
if there is an approximate discrete symmetry in the high-energy particle-physics 
model [60]. We have already seen an example of an approximate discrete symmetry 
inthe SU (5) Grand Unification model discussed in Section 2.1. If the cubic coupling 
in the potential in Eq. (2.5) is small, it can be ignored and the resulting model has 
SU(5) x Z2 symmetry, with all the domain wall solutions discussed in Section 2.2. 
A simpler example is that of the A#* together with a small cubic term. The potential 
is 


m? 2 3, À 4 
venye +ymọ Tg’ (8.61) 


Now the model still has two local minima but they are not exactly degenerate if y 
is very small (see Fig. 6.11). At the phase transition, a network of domain walls is 
formed and the typical separation and curvature scale of the domain walls is given 
by the correlation length &). With time the curvature scale grows and is denoted by 
R(t). So the force per unit area on the wall owing to tension is ~ o / R(t) where ø is 
the energy density of the wall. There is also a pressure difference pushing the wall 
toward the vacuum with the lowest energy. This pressure is given by the energy 
difference between the vacua and hence is proportional to y 


p~ ymi? (8.62) 


where 77 is the vacuum expectation value of the field. Therefore the tension is much 
larger than the pressure, and the dynamics of the wall network are unaffected by 


the pressure difference coming from the cubic term as long as 
o 1 
AE (8.63) 
ymn yn 


Once the network has evolved to a point where this condition is not met, the pressure 


R(t) < 


becomes important and drives the domain wall network such that the whole system 
reaches the true vacuum. From the area scaling law in Eq. (8.48), it follows that 
R(t) grows linearly with conformal time.* Therefore, in a radiation dominated 
universe, 


z +\12 
R = Ro— = r( 2) (8.64) 
To to 


4 The scaling law holds at late times after friction becomes unimportant. At earlier times, R grows as a different 
power of conformal time [87]. 
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where Ro can be taken to be the correlation length at the time of the phase transition 
fo. Inserting this relation in Eq. (8.63) and using Ro ~ 1/7, we get that the walls 
survive for a duration 


t 
balls ~ —5 (8.65) 
Y 


If y is small, the walls can survive for many Hubble expansions. In fact, if the walls 
survive for a long time, they might start dominating the density of the universe 
before they disappear. 

Even if domain walls are present in the universe for a relatively short time, they 
can still have important implications for cosmology. As the wall network evolves, 
the ambient matter interacts with the walls. Magnetic monopoles can get trapped on 
domain walls, leading to faster annihilation. This is the sweeping scenario discussed 
in [52]. In addition, the eventual collapse of domain walls can lead to black hole 
formation. These issues have received some attention but have yet to be studied in 
detail. 


8.16 Open questions 


jik 


. What happens when a black hole collides with a domain wall? Does it get stuck on the 
wall? Or does it pass through? For a discussion from the gravitational point of view, see 
[29, 148]. 

2. Develop an analytical formulation (perhaps along the lines of Section 3.5) to calculate 

the scalar radiation rate from collapsing domain walls. 

3. Discuss the cosmology of superconducting domain walls. 

4. What is the outcome of the SU (5) Grand Unified phase transition when the cubic coupling 

is small? Are domain walls formed? How does the network evolve? 


9 
Kinks in the laboratory 


In this chapter we discuss two laboratory systems where kinks are known to exist. 
The first system is trans polyacetylene which has a broken Z2 symmetry as in the 
i.¢* model. The second system is a Josephson junction transmission line, which is 
a laboratory realization of the sine-Gordon system. Helium-3 is another laboratory 
system that contains a wide variety of topological defects and the reader is referred to 
[174] for a discussion. In the third section of this chapter we describe Scott Russell’s 
solitons in water. These solitons are not topological like the others discussed in this 
book but we include the discussion anyway since the reader’s curiosity may have 
been aroused by the story in the Preface. 


9.1 Polyacetylene 


Polyacetylene consists of a linear chain of CH bonds. A sequence of x units is 
written as (CH),. In the ground state of polyacetylene, the carbon atom forms three 
o bonds, one of them is to the H in the CH unit, one to the unit on the left and one 
to the right. In addition, there is one more electron orbital that can cause bonding. 
This is called the x electron, and the x bond can form to the left or to the right. 
Then there are two possible sequences — first when the double (o and zr) bond is 
to the carbon on the right and the single to the left, the second when the double 
bond is to the left and the single to the right. These two possibilities are illustrated 
in Fig. 9.1 [149] in the trans configuration of polyacetylene.! 

The average bond lengtha ~ 1.22 A but the CH units are displaced so as to make 
double bonds (slightly) shorter than the single bonds. The physical displacements un 
along the horizontal axis in the two structures are depicted in Fig. 9.1. Qualitatively, 
the essential point is that the x electrons have to choose to either form the double 


1 In the cis configuration, there are also two states related by the left-right transformation but they are not degenerate 
in energy. 
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Figure 9.1 Structure of the two degenerate ground states of trans polyacetylene. 
The upper structure is denoted by A and the lower by B. Double bonds are denoted 
by heavy lines. 
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Figure 9.2 If the B state occurs on the left side of a chain and the A on the 
right, there is a kink in between where the simple alternating structure cannot be 
maintained. The kink is the region where the alternate single-double bonds do not 
exist. 


bond to the left or to the right. Hence there is a Z2 symmetry which is broken in the 
(“dimerized”’) ground state. Kinks form if different ground states are chosen at dif- 
ferent locations (Fig. 9.2). The center of the kink is located at the CH unit where the 
x electron wavefunction is equally shared between the CH units to the left and right. 

The Hamiltonian of the system depends on the displacement variables, u„ and 


on the locations of the z electrons 
K M 
H=- Ponneau Hed +D Slag wah + Sah OD 


where 


tn+1,n = fo — Æ (Un+1 — Un) (9.2) 
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is the hopping integral to leading order in displacements. The operators ch and 
Cn,s are creation and annihilation operators for electrons of spin s on the nth CH 
group. The parameter K is the effective spring constant of the o bonds and M is 
the mass of the CH group. 

To connect with the discussion of Chapter 1, the displacement variable 


dn = (Dun (9.3) 


can be viewed as a scalar field defined on a lattice interacting with a fermion (the z 
electron). The last two terms in Eq. (9.1) correspond to gradient and time derivative 
terms of a continuum field ¢(x) that corresponds to the discrete variables, ¢,. The 
first term describes interactions between ¢ and the electrons. The effective interac- 
tion for the @ field, after integrating out the fermionic variables, must respect the 
Zz symmetry, and hence corresponds to a ¢* interaction to lowest order. Therefore 
a non-relativistic version of the Z2 model of Eq. (1.2) captures some of the gross 
features of polyacetylene. 

The properties of kinks in polyacetylene (Fig. 9.2) have been studied in [150] 
using the Hamiltonian in Eq. (9.1) with the result that the kink width is approxi- 
mately 14 lattice spacings and the mass is approximately six electron masses [152] 
in good agreement with experiments [75]. 

The quantum properties of polyacetylene kinks have also been studied. In Sec- 
tion 5.3 we discussed how kinks can carry fractional quantum numbers [83]. Poly- 
acetylene kinks also carry fractional quantum numbers and electric charge, namely 
“half a bond” or +(2e)/2 charge since each bond consists of two electrons (one 
from each atom at either end of a bond) [150]. Indeed, in a chain where two single 
bonds are followed by a double bond (instead of the alternating single and double 
bonds in trans polyacetylene) the fractional charge can be shown to be one-third of 
a bond [151]. Reference [68] generalizes these ideas much further and shows that 
solitons may even carry irrational charges. 


9.2 Josephson junction transmission line 


We follow [134] in deriving the sine-Gordon equation for the Josephson transmis- 
sion line. 

Let us recall the basics of a transmission line, schematically shown in Fig. 9.3 
[54]. A potential difference is applied to the ends of two elements of a transmission 
line e.g. the two cables of a coaxial cable. The potential, V, and current, Z, in 
each of the wires are functions of the location on the transmission line, namely 
the x coordinate, and also of time. There is also a potential difference between the 
wires, and the current in the two wires can be different, but this is not shown in 
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Figure 9.3 Schematics of a transmission line with inductance Lo and capacitance 
Co per unit length. The symbols marked J represent a second coupling between the 
two transmission components. This coupling is absent in an ordinary transmission 
line but represents the tunneling current in the Josephson transmission line. 


the figure. We will only be considering the potential and current distribution along 
a single wire. Let Lo denote the inductance per unit length of the line, and Co 
the capacitance per unit length. Then Faraday’s law of induction tells us that the 
induced e.m.f. between points x + dx and x is proportional to the rate of change 
of the current in the segment within those points 


al 
Vit,x + dx) = V(r, x) = (Lode) = (9.4) 
or 
av al 
Bisa eee exc 9.5 
ax ET (23) 


Charge accumulates on the segment from x to x + dx in time dt owing to the 
different entering and exiting currents. The charge on the segment is also given by 
the capacitance times the potential. Hence 


av 
I(t, x + dx) = 1, x) = —(Codx) = (9.6) 
or 
ðI av 
pcan gies 9.7 
ax ar AT) 


Equations (9.5) and (9.7) can be combined to obtain wave equations for the current 
and the potential. 

A Josephson junction transmission line differs from the ordinary transmission 
line described above in that the two “wires” are superconductors and they are 
separated by a thin insulator. This set-up is shown in Fig. 9.4. Current can tunnel 
through the insulator and jump from one wire to the other. Hence the charge on a 
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Figure 9.4 A Josephson junction transmission line is constructed by separating 
two superconducting plates by a thin layer of an insulating material. 


segment also changes owing to the Josephson current and Eq. (9.7) gets modified 


ol aV 
— = —Co— — jx, t 9.8 
ae 0% hœ, t) (9.8) 
where jy is the Josephson current per unit length. 

The charge carriers (Cooper pairs) in either superconductor are described by 


macroscopic wavefunctions 


yı = Joie, th = ype” (9.9) 
where pı and pz are the charge carrier number densities in the superconductors. 
The tunneling Josephson current per unit area is [55] 


Ji = josing (9.10) 
where jo is the maximum Josephson current and is proportional to ./p; 62, and 
$=¢ġ -h (9.11) 
The Schrödinger equations for Yı and y2 imply 
dp q 
aV 9.12 
or Rh A 


where V is the potential difference across the junction and q = 2e is the electric 
charge of a Cooper pair. 
From Eqs. (9.5) and (9.12) we obtain 


| ene (9.13) 
Now we can insert this expression for J, and V as found from Eq. (9.12), in Eq. (9.8) 


to get 


ao 1 aod jq 
2 ing =0 9.14 
E LGI oa’ oe) 
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Rescaling t and x by the Josephson time and length scales 


RCA V2 h \ 12 
7 =) v- AE ( ) (0.15) 
qJo q JoLo 


gives the sine-Gordon equation as derived from Eq. (1.51) witha = 1 = £. 


9.3 Solitons in shallow water 


The solitons discussed in this book have all had a topological origin. In contrast, 
the solitons first discovered by Scott Russell in a water channel, and mentioned in 
the Preface, have their origin in the non-linearities of hydrodynamics and do not 
have a topological origin. 

The first step to show the existence of the water solitons is to derive the Korteweg- 
deVries (KdV) equation for waves of long wavelength moving in one direction in 
shallow water. We do not give this derivation here and instead refer the reader to, 
for example, Section 13.11 of [180]. 

The KdV equation is 


3 
ə əx 3x? 


where u is related to the height of the fluid surface and ô is a parameter. The soliton 
solution is [185] 


Jw- 8 3 
“tu “4p 0 (9.16) 


2| X — Xo 

u = Ugo + (uo — uw )sech = (9.17) 

where X = x — vt, Xo is a constant, v is the velocity of the soliton, 

iTA 
E A (9.18) 
12 
The velocity of the soliton is given by 

peung a (9.19) 


3 


in terms of the arbitrary constants uo and us. Note that the amplitude of the soliton 
and the velocity are related. 


9.4 Concluding remarks 


There are a number of situations where solitons have been discussed in the par- 
ticle physics literature. Most of these discussions, such as of domain walls in the 
SU(5) x Z2 model in Chapter 2, have been in the framework of Grand Unified 
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Theories. The attention has mostly focused on magnetic monopoles and strings 
because monopoles seem inevitable in this class of theories and strings are less 
constrained by cosmology. Similar topological structures also exist in the standard 
model of electroweak interactions but the monopoles are confined and the strings 
are unstable [2]. Domain wall and string solutions also exist in QCD in various 
external conditions, for example in high density matter such as might be present 
in the interiors of neutron stars [142]. Unlike solitons in the laboratory, however, 
solitons in particle physics and cosmology have not yet been discovered experi- 
mentally. Given the very similar underpinnings of laboratory and particle physics 
systems, there is hope that this situation will soon change. 


9.5 Open questions 


1. Is there a condensed matter system with spontaneously broken permutation symmetry? 
Discuss the domain walls in that system and whether a lattice can exist. Can the walls 
be observed experimentally? 

2. If there are QCD domain walls in neutron stars, how might they be observed from Earth? 


Appendix A 


Units, numbers and conventions 


We will work in natural units in which h = c = 1. In these units, all dimensionful 
quantities have dimensions of mass to some power. One way to convert from mass (g) to 
length (cm) and time (s), is to remember the values for the Planck mass, time, and length: 
mp = 1.2 x 10!° GeV, tp = 5.4 x 10° s, lp = 1.6 x 107? cm. Also, mptp = 1 = mplp 
in natural units. It is also useful to remember mp = 2.2 x 1075 g and, when dealing with 
magnetic fields, the conversion: 1 Gauss = 1.95 x 107% GeV?. In addition, for 
cosmological estimates it is convenient to know that 1 pe = 3.1 x 10!8 cm. 

The metric signature is taken to be (+, —, —, —). 
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Appendix B 
SU(N) generators 


SU(N) is the group of special (unit determinant), unitary, N x N complex matrices.' By 
considering the various constraints on the 2N? real components of the matrix owing to the 
special and unitary conditions, we can see that the matrix has N? — 1 independent degrees 
of freedom. Then, if g e SU(N), we can write 


g = exp(iagT“) (B.1) 


where a sum overa = 1,..., N? — 1 is implicit, w, are real constants, and T“ are the 
“generators” of the group. The T° satisfy the SU(N) Lie algebra and can be represented 
by matrices of various dimensions. In the N = 2 (SU (2)) case, the two-dimensional 
representation is in terms of Pauli spin matrices, T = of /2, or explicitly 


pa m e 9) oo 


The Lie algebra is 
[T", T?] = ie T°? (B.3) 


where €% is the totally antisymmetric tensor. One can also easily construct the higher 
dimensional representations. It is conventional to normalize the generators to satisfy 


1 
Tr(T'T®) = ad (B.4) 


where 5“ is the Kronecker delta. 

To get a set of generators for SU(N), it is simplest to build on the SU(2) generators in 
Eq. (B.2). First, one puts the Pauli spin matrices in the upper left-hand corner and obtains 
three SU(N) generators 


a l/o 0 
Paho 0 
Then one puts the off-diagonal Pauli spin matrices in the off-diagonal positions. Since 
there are N(N — 1)/2 off-diagonal positions of which two have already been filled by the 
a = 1, 2 generators, we can construct N(N — 1) — 2 more generators by filling each 


me a=1,2,3 (B.5) 


1 For a review of group theory in particle physics, see [62]. 
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remaining position by either 1 (as ino!) or by +i (as in o°). These look like 
On Or Beso tga.) Sax 0 0 


=e | æ 


N 
© 
— 
= 
Nie 
o 
= 
= 


where the subscripts j, k denote the position in the matrix. 

Finally we construct the diagonal generators. These are written by putting a series of 1s 
in say, n, successive diagonal positions, and then entering —n in the nn entry of the 
matrix. This scheme ensures that the generator is traceless and the resulting matrix is 


diag(1,...,1,,—n, 0,..., 0) (B.7) 


where 1, denotes 1 in the nn entry. The normalization is then fixed using the convention in 
Eq. (B.4) to get the generator 


1 
——diag(1,..., ln, —n, 0,...,0 B.8 
mae S } i 

In this way we construct N — 1 diagonal generators, one for each value of n. The third 
Pauli matrix is already included as the a = 3 generator. 

As a check, we find that the total number of generators constructed is 
3+(N(N — 1)— 2)+ (N — 2) = N? — Land this agrees with the degrees of freedom in 
SU(N). 

In the SU (5) Grand Unified model discussed in Chapter 2 an alternate set of diagonal 
generators is useful. 


1 
À3 = ztagl, —1, 0,0,0) 


1 
àg = —diag(1, 1, —2, 0, 0) 
aa S 


1 
T3 7 diag, 0,0, 1, —1) 


1 
Y = ——diag(?2, 2, 2, —3, —3 
2/15 a ) 


After the SU (5) symmetry is broken by the canonical vacuum expectation value of ® 
(Eq. (2.6)), A3 and Ag are generators of the unbroken SU (3), t3 of SU(2), and Y of U(1). 


Appendix C 


Solution to a common differential equation 


We have often encountered a differential equation of the type 


ey ; 2 2 
-Z + [e — vcosh 2u — v sinh 2u tanh x + vcosh* u sech x] y=0 


(C.1) 


where v, u are parameters and € is the eigenvalue. This differential equation has been 
solved in Section 12.3 of [113] where the Schrödinger problem has also been extensively 


studied. Here we reproduce the solution. 
The solution is given in terms of new parameters a and b 


1 1 1 1 
a= a vert — e — aN veW# — e = z+ T 3- 


2 
1 1 1 1 
b = -y ve — e + ~Vve 7 — € = >k} + =k- 
2 2 2 2 
Then, with 
Y = e™®sechtx F(x) 
the equation for F becomes 
F" — 2[a + b tanh x] F’ + [v cosh? u — b(b + 1)]sech?x F = 0 


where primes denote derivatives with respect to x. Defining 
1 
u = 5l! — tanh x] 


we get the hypergeometric equation 
dF dF 3 
u(1 -uz T [a +b+1-—2(b+ Duly + [v cosh" u — b(b + 1)]F = 0 
u u 
The general solution may be found in [71] 
F=AF\+ BF 
where A and B are constants of integration and 
F, = F(a, p; y; u) 
F =u" Fœ -y +1,ß -y +1;2-—y;u) 
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(C.2) 


(C.3) 


(C.4) 


(C.5) 


(C.6) 


(C.7) 


(C.8) 


(C.9) 
(C.10) 
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where 
i= \/ vcosh? ju + 
a= = — ,/ vcos = 
2 KI 
1 2 1 
p=b+5+ vcosh Big (C.11) 
y=at+b+l1 (C.12) 
and y is assumed to not be an integer. 
The general analysis can be taken further by considering the solution at x = +00. A 
solution that is regular at x — 00 (i.e. u = 0) is obtained by setting B = 0 in Eq. (C.8). 
Regularity at x = —oo (u = 1) is only obtained for certain values of €, and thus the 


energy levels are quantized. The details of the general analysis may be found in 
Section 12.3 of [113]. 

In this book, we have often encountered the special case with u = 0. Then, bound 
states are obtained for the following discrete values of b > 0 


by =o g— (45) (C.13) 


where n = 0, 1,2,..., N with N determined by by; < 0. The discrete eigenvalues of e€ 
follow from the definition in Eq. (C.3) 


1 1 
€n = (2n +1) v+- (+a) (C.14) 


Appendix D 


Useful operator identities 


Identity 1 
We wish to prove! 
eAtB — e^eBeC/2 


where C = [B, A] is assumed to commute with A and B. 
Let 


Sœ) = e(AtB)x 
where x is a parameter. Write 
S(x) =e“ U(x) 
where U is an unknown matrix-valued function. Then 


ds dU 
(A + B)S(x) = — = Ae“ U(x) + e4* — 
dx dx 


which leads to 


dU 
— =e “ Bet“ U(x) 
dx 
Now 
ay O'R [B, A" ]x” 
ee er tee ee er 
Also, 
[B, A"] = A”'[B, A] + [B, A” JA =--- = [B, A]n A"! 
provided C = [B, A] commutes with A. Therefore 
Be** = e® B + [B, A]xe” 
and 
dU 
—— = (B + [B, A]x)U(x) 
dx 


1 Several of the proofs in this Appendix were provided by Harsh Mathur, private communication (2005). 
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(D.1) 


(D.2) 


(D.3) 


(D.4) 


(D.5) 


(D.6) 


(D.7) 


(D.8) 


(D.9) 
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This equation can be solved to get 
U(x) = exp(Bx + [B, A]x?/2) (D.10) 


This solution satisfies the boundary condition U (0) = 1. 
Hence 


S(x) = e(AtB)x — e4% eBx lB. Al? /2 (D.11) 


With x = 1, we get the desired result. 


Identity 2 


Here we outline a proof of the identity 
:e4tB e ef ne: eP (D.12) 


where D = [A™, BT] is assumed to be a c-number. Also, a linear decomposition of A and 
B is assumed A = At + A~, B = B* + B- and the superscripts + refer to terms 
proportional to creation (+) and annihilation (—) operators. Expressions sandwiched 
between : : are normal ordered, and so annihilation operators are placed to the right of 
creation operators. 

The first step is the identity 


ef := ei et (D.13) 
This can be proved by explicit expansion of the exponentials. 
Then 
ref ne? me^ e^ e e" (D.14) 
and 
:e^tB = ef eb e^ eB" (D.15) 


since [A~, B7] = 0 = [A™, Bt]. 

Now we use the identity, Eq. (D.1) proved in the previous section to exchange the order 
of the middle two factors in Eq. (D.14) and together with Eq. (D.15) gives the identity in 
Eq. (D.12). 


Identity 3 
Here we wish to show 
A:e® :=:{A+[At, B]}e® : (D.16) 
i B. œ (4+) z n! —~\k ptyn—k 
A:e Tae J ia er? (Bt) (D.17) 
=) 1 ; n! k a4 k 
=A Series eae AEE aed 


+A~(B-)} (B+ (D.18) 


Useful operator identities 


Now use 
[A*, (B7)*] = k(B7)'[A*, B7] 
provided that [At, B~] is a c-number. Therefore 


(oe) 1 n 


! 
Are 2 ae co EAT + k(B7)E![At, B7] 


n=0 "* k=0 
+A (BY (Bt 
=: (A + [A}, B Je? : 
Similarly we can show 
:e? : A =: eP (A + [B*, AT]: 
Putting this together with Eq. (D.16) we get 
[A, : e? :] =: [A,e?]:+({At, B7] —[B*, AT]) : e? : 

In the case of interest for deriving Eq. (4.80) we have 


A*t =(A7), B* = -(B-) 


then 


[At, B7] —[Bt, A7] = [At, B7] + [At, BT] 
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(D.19) 


(D.20) 
(D.21) 


(D.22) 


(D.23) 


(D.24) 


(D.25) 


With A = ¢(y) and: e? := y(x), [At, BT] is purely imaginary and the right-hand side 


vanishes. Then Eq. (D.23) gives the identity 
[A, : e? :] =: [A,e?]: 


(D.26) 


Appendix E 


Variation of the determinant 


If M is a matrix function and 6M is a small variation of M, we wish to find the variation 
of the determinant of M (we follow Section 4.7 of [177]). 
Consider 


8[In(DetM(x))] = In(Det(M + 6M)) — In(Det(M)) 


_ p [PaM + 8m) 
= af Det(M) | 


= In(DetM~'Det(M + 6M)) 
= In(Det{M~!(M + 6M)}) 
= In(Det{1 + M~'sM}) 
= In(1 + Tr{M~!5M}) + OEM 
= Tr{M~'5M} + O((6M)*) (E.1) 
Hence 
5[(DetM(x))] = Tr[M~!(x)5M(x)|DetM (x) (E.2) 


which is the desired result. 
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Appendix F 


Summary of cosmological equations 


Assuming an isotropic and homogeneous universe, the cosmological line element is the 
Friedman-Robertson- Walker line element, and can be written as: 


dr 
X 442 2 
ds“ = dt - | 


+ r7(d6? + sin? 06?) (F.1) 
The function a(t) is known as the scale factor. k is a parameter that is —1 for a hyperbolic 
or negatively curved universe, 0 for a flat universe, and +1 for a positively curved 
universe. 

The equations of motion for a are derived from Einstein’s equations assuming that the 
universe is filled with one or more fluids! with total energy density p and pressure p. 


Then: 
ey) 
a 8G k 
HR? = = F.2 
(2) z ee (F.2) 
Ms 4r G 
ä = Sg +3p)a (F.3) 


and energy-momentum conservation gives: 


d 3 2 
aa 4 ) = —3 pa (F4) 


In addition, we need the equation of state for the fluid to connect p and p. Some examples 
of equations of state are: p = —p (cosmological constant), p = 0 (dust), p = p/3 
(radiation), and p = —2p/3 (slowly evolving wall network). Note that o may contain 
contributions from a large variety of forms of matter and then the corresponding energy 
densities and pressures must be added together. 


1 The fluid approximation means that the relaxation time of the various components — which in our case may 
be plasma, gas, stars, galaxies, walls, fundamental particles — is much shorter than the characteristic time for 
changes in the scale factor. 
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Summary of cosmological equations 


Assuming a single dominant component of energy density in a flat universe (k = 0), 
Eqs. (F.2) and (F.4) can be solved to obtain: 


Ht 0 


p=-prmaxe”, pxa 

1 

p=0>aat?, px 

a 

1 

=e 1/2 = 

p= 7z>7acaxt er 

2 1 
p=-3P> aat, P (F.5) 


p=wp ra xP Bey. px q3ury (F.6) 
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